A highly accurate equation of state (EOS) for chain molecules formed from spherical segments interacting through Mie potentials (i.e., a generalized Lennard-Jones form with variable repulsive and attractive exponents) is presented. The quality of the theoretical description of the vapour-liquid equilibria (coexistence densities and vapour pressures) and the second-derivative thermophysical properties (heat capacities, isobaric thermal expansivities, and speed of sound) are critically assessed by comparison with molecular simulation and with experimental data of representative real substances. Our new EOS represents a notable improvement with respect to previous versions of the statistical associating fluid theory for variable range interactions (SAFT-VR) of the generic Mie form. The approach makes rigorous use of the Barker and Henderson high-temperature perturbation expansion up to third order in the free energy of the monomer Mie system. The radial distribution function of the reference monomer fluid, which is a prerequisite for the representation of the properties of the fluid of Mie chains within a Wertheim first-order thermodynamic perturbation theory (TPT1), is calculated from a second-order expansion. The resulting SAFT-VR Mie EOS can now be applied to molecular fluids characterized by a broad range of interactions spanning from soft to very repulsive and short-ranged Mie potentials. A good representation of the corresponding molecular-simulation data is achieved for model monomer and chain fluids. When applied to the particular case of the ubiquitous Lennard-Jones potential, our rigorous description of the thermodynamic properties is of equivalent quality to that obtained with the empirical EOSs for LJ monomer (EOS of Johnson et al.) and LJ chain (soft-SAFT) fluids. A key feature of our reformulated SAFT-VR approach is the greatly enhanced accuracy in the near-critical region for chain molecules. This attribute, combined with the accurate modeling of second-derivative properties, allows for a much improved global representation of the thermodynamic properties and fluid-phase equilibria of pure fluids and their mixtures.