Top panels: Energy dependence comparison of the transverse momentum spectra of K*$^0$ (left) and $\phi$ (right) as a function of $p_{\rm T}$ for the NSD events, measured in the rapidity interval $-$0.5 $<$ $y$ $<$ 0 for p$-$Pb collisions at $\sqrt{s_{\rm NN}} =$ 5.02 and 8.16 TeV. Bottom panels: The ratio of $p_{\rm T}$ spectrum at $\sqrt{s_{\rm NN}} =$ 8.16 TeV to the $p_{\rm T}$ spectrum at $\sqrt{s_{\rm NN}} =$ 5.02 TeV. The ratio is compared with the predictions from EPOS-LHC, DPMJET and HIJING. The statistical and systematic uncertainties are shown as bars and boxes, respectively.