Recent works on three-planet mean motion resonances (MMRs) have highlighted their importance for understanding the details of the dynamics of planet formation and evolution. While the dynamics of two-planet MMRs are well understood and approximately described by a one-degree-of-freedom Hamiltonian, little is known of the exact dynamics of three-body resonances besides the cases of zeroth-order MMRs or when one of the bodies is a test particle. In this work, I propose the first general integrable model for first-order three-planet mean motion resonances. I show that one can generalize the strategy proposed in the two-planet case to obtain a one-degree-of-freedom Hamiltonian. The dynamics of these resonances are governed by the second fundamental model of resonance. The model is valid for any mass ratio between the planets and for every first-order resonance. I show the agreement of the analytical model with numerical simulations. As examples of application, I show how this model could improve our understanding of the capture into MMRs as well as their role in the stability of planetary systems.