Simulating hypersonic flow around a space vehicle is challenging because of the multi -scale and nonequilibrium nature inherent in these flows. To effectively deal with such flows, a hybrid scheme combining the stochastic particle Bhatnagar-Gross-Krook (BGK) method with direct simulation Monte Carlo (DSMC) was developed recently, but only for monatomic gases (Fei et al. (2021) [29]). In this paper, the particle-particle hybrid method is extended to polyatomic gas flows. In the near continuum regime, employing the Ellipsoidal-Statistical BGK model proposed by Dauvois et al. (2021) [22] with discrete lev-els of vibrational energy, the stochastic particle BGK method for polyatomic gases is first established following the idea of the unified stochastic particle BGK (USP-BGK) scheme. It is proven to be of second-order accuracy in the fluid limit. After that, the USP-BGK scheme with rotational and vibrational energies is combined with DSMC to construct a hybrid scheme. The present hybrid scheme for polyatomic gases is validated with numer-ical tests of homogeneous relaxation, 1D shock structure and 2D hypersonic flows past a wedge and a cylinder. Compared to the other stochastic particle methods, the proposed hybrid scheme can achieve higher accuracy at a much lower computational cost. Therefore, it is a more efficient tool to study multiscale hypersonic flows.(c) 2022 Elsevier Inc. All rights reserved.