Effective preconditioning of neutron diffusion problems is necessary for the development of efficient DSA schemes for neutron transport problems. This paper uses P-multigrid techniques to expand two preconditioners designed to solve the MIP diffusion neutron diffusion equation with a discontinuous Galerkin (DG-FEM) framework using first-order elements. These preconditioners are based on projecting the first-order DG-FEM formulation to either a linear continuous or a constant discontinuous FEM system. The P-multigrid expansion allows the preconditioners to be applied to problems discretised with second and higher-order elements. The preconditioning algorithms are defined in the form of both a V-cycle and W-cycle and applied to solve challenging neutron diffusion problems. In addition a hybrid preconditioner using P-multigrid and AMG without a constant or continuous coarsening is used. Their performance is measured against a computationally efficient standard algebraic multigrid preconditioner. The results obtained demonstrate that all preconditioners studied in this paper provide good convergence with the continuous method generally being the most computationally efficient. In terms of memory requirements the preconditioners studied significantly outperform the AMG.