The development of optimal, or near optimal solution strategies for higher-order discretizations, including steady-state solutions methodologies, and implicit time integration strategies, remains one of the key deter-mining factors in devising higher-order methods which are not just competitive but superior to lower-order methods in overall accuracy and efficiency. The goal of this work is to investigate and develop a fast and robust algorithm for the solution of high-order accurate discontinuous Galerkin discretizations of non-linear systems of conservation laws on unstructured grids. Herein we extend our previous work to the three-dimensional steady-state Euler equations, by coupling the spectral p-multigrid approach with a more traditional agglom-eration h-multigrid method for hybrid meshes, in a full-multigrid iteration strategy. In this hp-multigrid ap-proach the coarse "grid" levels are constructed by reducing the order (p) of approximation of the discretization using hierarchical basis functions (p-multigrid), together with the traditional (h-multigrid) approach of con-structing coarser grids with fewer elements. The overall goal is the development of a solution algorithm which delivers convergence rates which are independent of "p" (the order of accuracy of the discretization) and inde-pendent of "h" (the degree of mesh resolution), while minimizing the cost of each iteration. The investigation of efficient smoothers to be used at each level of the multigrid algorithm is also pursued, and comparisons between different integration strategies are made as well. Current three-dimensional results demonstrate con-vergence rates which are independent of both order of accuracy (p) of the discretization and level of mesh resolution (h).