We derive the single-fluid evolution equations describing a mixture made of a gas phase and an arbitrary number of dust phases, generalizing the approach developed by Laibe & Price. A generalization for continuous dust distributions as well as analytic approximations for strong drag regimes is also provided. This formalism lays the foundation for numerical simulations of dust populations in a wide range of astrophysical systems while avoiding limitations associated with a multiple-fluid treatment. The usefulness of the formalism is illustrated on a series of analytical problems, namely the DUSTYBOX, DUSTYSHOCK and DUSTYWAVE problems as well as the radial drift of grains and the streaming instability in protoplanetary discs. We find phys...