We consider the problems whose mathematical model is determined by some Markov chain terminating with probability one; moreover, we have to estimate linear functionals of a solution to an integral equation of the second kind with the corresponding substochastic kernel and free term [1]. To construct weighted modifications of numerical statistical models, we supplement the coordinates of the phase space with auxiliary variables whose random values functionally define the transitions in the initial chain. Having implemented each auxiliary random variable, we multiply the weight by the ratio of the corresponding densities of the initial and numerically modeled distributions. We solve the minimization problem for the variances of estimators of linear functionals by choosing the modeled distribution of the first auxiliary random variable.