In this article, we have reported the transverse momentum ($p_T$) spectra and nuclear modification factors of (multi-) strange hadrons produced in Au+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV using the Monte Carlo HYDJET++ model. For the present study, we have tuned few of the input parameters, in both soft as well as hard parts, of the model to appropriately describe the production of (multi-) strange hadrons. We have presented the $p_T$-spectra of (multi-) strange hadrons in different centrality intervals. On comparing $p_T$-spectra with the experimental data, we observe that the model does not achieve thermal equilibrium for (multi-) strange hadrons in peripheral collisions and at high $p_T$. The $R_{AA}$ results obtained by the model reproduce the experimental data well while the $R_{CP}$ results show a larger suppression at high $p_T$ than the experimental data. Further, HYDJET++ results are also compared with different theoretical models and discussed, wherever possible.