Based on inverse scattering transformation, a variable-coefficient fifth-order nonlinear Schrödinger equation is studied through the Riemann–Hilbert (RH) approach with zero boundary conditions at infinity, and its multi-soliton solutions with [Formula: see text] distinct arbitrary-order poles are successfully derived. By deriving the eigenfunction and scattering matrix, and revealing their properties, a RH problem (RHP) is constructed based on inverse scattering transformation. Via solving the RHP, the formulae of multi-soliton solutions are displayed when the reflection coefficient possesses [Formula: see text] distinct arbitrary-order poles. Finally, some appropriate parameters are selected to analyze the interaction of multi-soliton solutions graphically.