The productions of massive gauge bosons, Z 0 and W + /W − , in heavy-ion reactions at the LHC, provide an excellent tool to study the cold nuclear matter effects in high-energy nuclear collisions. In this paper we investigate Z 0 and W + /W − productions in p+Pb and Pb+Pb at the LHC, at NLO and NNLO with DYNNLO incorporating the nuclear PDFs (nPDFs) parametrization sets EPS09 and DSSZ, within the framework of perturbative QCD. The numerical simulations of the transverse momentum spectra, rapidity dependence, and related nuclear modification factors for Z 0 and W particles, as well as the charge asymmetry for W boson, are provided and tested against the latest experimental data. It is found that the theoretical results with EPS09 and DSSZ nPDFs can give good descriptions of the recent data on Z 0 and W ± particles in p+Pb and Pb+Pb within the experimental error bars, though some differences between results with EPS09 and DSSZ can be observed, especially in the rapidity dependence of the Z 0 yield. Theoretical predictions for future measurements on Z 0 and W in p+Pb and Pb+Pb collisions at the LHC are also provided.