This is a brief review of recent progress in constructing solutions to the matrix model Virasoro equations. These equations are parameterized by a degree n polynomial W n (x), and the general solution is labeled by an arbitrary function of n − 1 coefficients of the polynomial. We also discuss in this general framework a special class of (multi-cut) solutions recently studied in the context of N = 1 supersymmetric gauge theories.Introduction. It was realized in the beginning of nineties that matrix models partition functions typically satisfy an infinite set of Virasoro-like equations [1,2]. These were nothing but Ward identities (Schwinger-Dyson equations) which mainly fixed matrix model partition functions (because of the topological nature of matrix models [3], the Ward identities were restrictive enough). Moreover, it turned out that one of the most technically effective ways to deal with matrix models was to solve these Virasoro equations (they are also sometimes called loop equations) [4,5,6,7,8].At early times of matrix models one usually dealt with Virasoro equations describing relatively simple "phases" so that the equations had unambiguous solutions. An interest to more complicated phases of matrix models has revived after G.Bonnet, F.David and B.Eynard [9] proposed to deal with the multi-cut (or multi-support) solutions (known for long, [10,11,6,12]) in a new way: releasing the tunneling constraint [11]. Their approach was later applied by R.Dijkgraaf and C.Vafa [13] to description of low energy superpotentials in N = 1 SUSY gauge theories, [14].More concretely, the authors of [14] considered the N = 2 SUSY gauge (Seiberg-Witten) theory in special points where some BPS states become massless. Therefore, these states can condense in the vacuum which breaks half of the supersymmetries (leading to N = 1 SUSY) and gives rise to a non-zero superpotential. Values of this superpotential in minima are related to the prepotential of a Seiberg-Witten-like theory. In turn, R. Dijkgraaf and C.Vafa associated [13] the prepotential with logarithm of a partition function of the Hermitean one-matrix model in the planar limit of multi-cut solutions (it was later proved in [15]).In fact, actual definition of the multi-cut partition functions is a separate problem. For instance, one could simply define them as (arbitrary) solutions to the corresponding Virasoro equations (Dmodule point of view). Then, one may ask what is special about the concrete Dijkgraaf-Vafa (DV) §