The complete phase diagram of Random Energy Model (REM) is obtained for complex temperatures using the method proposed by Derrida. We find the density of zeroes for statistical sum. Then the method is applied to Generalized Random Energy Model (GREM). This allows us to propose new analytical method for investigating zeroes of statistical sum for finite-dimensional systems.Random Energy Model (REM) [1] is one of the most popular contemporary models of statistical mechanics. Besides its direct applications in spin-glasses the model has a wide range of employments in very different areas of the modern theoretical physics and biophysics. Recently, several new, and more intimate applications have been founded in the theories of well-developed turbulence [2], and strings [3][4][5]. On the other hand, it is well-known that the model allows to consider the fundamental problems of information theory in the language of statistical mechanics. In particular, problems connected with transmission of information through a noisy channel can be formulated and effectively solved in the language of REM [6,7]. This astonishing range of different applications cannot be accidental one: It seems, that REM is a paradigm for a complex system physics, or at least contains all essential ingredients of that physics. In other words, the nature sings his sacral songs using REM's language.The thermodynamical structure of REM has simple structure but at the same contains all essential ingradients of phase transitions [1,11,12]: There is well defined critical temperature T c which can be obtained from clear and semi-intuitive physical reasons. For T < T c the system is frozen in the spin-glass phase with one level of replica symmetry breaking. Besides this rough thermodynamical phase transition there are some "mild", mesoscopic phase transitions in the domain (T , T c ),T > T c [1,12]. The temperatureT strongly depends from the concrete form of energy's distribution; it can be infinite for the gaussian distribution [1] but it has some finite value for REM generated by more physical dilute hamiltonians [12]. This rich and interesting spectra of phase transitions waits its complete and exhaustive investigation, and we hope that the present work will shed some light to this question also.It is well-known that structure and properties of phase transitions can be investigated through analytical properties of thermodynamical quantities in the complex plane of temperature or/and magnetic field [18] [17] [16]. Indeed, if a phase transition is assosiated with singular behavior of thermodynamic potential in the thermodynamic limit (in our case it is free energy or statistical sum), then by consideration of its analytical properties in the complex plane the important physical information can be obtained. The method was proposed by Yang and Lee [14] (for the case of complex magnetic fields), and Fisher [15] (for the case of complex temperatures). It has large variety of different applications in statistical physics; furthermore it is one one of the most ex...