In this paper we consider the very high order approximation of solutions of the Euler equations. We present a systematic generalization of the residual distribution method of (Abgrall, J.ComputPhys 2006) to very high order of accuracy, by extending the preliminary work discussed in (Abgrall, Larat, Ricchiuto, Tave, Computers and Fluids 2009) to systems and hybrid meshes. We present extensive numerical validation for the third and fourth order cases with Lagrange finite elements. In particular, we demonstrate that we both have a non-oscillatory behavior, even for very strong shocks and complex flow patterns, and the expected accuracy on smooth problems.