State switching processes play an important role in the control of nanodevices and properties of many quasi-1D objects in physics, chemistry and biology. The switching a metastable state into an energetically preferable state in extended systems often occurs through the nucleation and extension of domain of the new phase. In real systems new domains arise as a rule on defects or sample edges. Usually considered limit cases are the propagation of the edge domain wall over the whole system or multiple nucleation and coalescence of domains in the volume of the material. This paper deals with the more general situation of front propagation of the switching state from the sample edge simultaneously with the multiple nucleation of the new phase domains on randomly located defects in the volume. The distinction of the considered model from a large number of close ones is the account of the combination of two factors: (1) the study of the propagation of the state switching front from the sample boundary in a semi-infinite system, and (2) a detailed statistical description of the front runs in the fluctuating medium in terms of distribution functions, that is, beyond the framework of the mean-field approach. The corresponding statistically-kinetic problem of 'the relay' run lengths of the state switching front is exactly solved with the calculation of the distributions of the PAPER: Classical statistical mechanics, equilibrium and non-equilibrium