Recent experiments have observed large anisotropic collective flows in high multiplicity proton-lead collisions at the Large Hadron Collider (LHC), which indicates the possible formation of mini quark-gluon plasma (QGP) in small collision systems. However, no jet quenching has been confirmed in such small systems so far. To understand this intriguing result, the system size scan experiments have been proposed to bridge the gap between large and small systems. In this work, we perform a systematic study on both heavy and light flavor jet quenching in different collision systems at the LHC energies. Using our state-of-the-art jet quenching model, which combines the next-to-leading-order perturbative QCD framework, a linear Boltzmann transport model and the (3+1)-dimensional viscous hydrodynamics simulation, we provide a good description of nuclear modification factor $R_{rm AA}$ for charged hadrons and $D$ mesons in central and mid-central Pb+Pb and Xe+Xe collisions measured by CMS collaboration. We further predict the transverse momentum and centrality dependences of $R_{AA}$ for charged hadrons, $D$ and $B$ mesons in Pb+Pb, Xe+Xe, Ar+Ar and O+O collisions at the LHC energies. Our numerical results show a clear system size dependence for both light and heavy flavor hadron $R_{AA}$ across different collision systems. Sizable jet quenching effect is obtained for both heavy and light flavor hadrons in central O+O collisions at the LHC energies. Our study provides a significant bridge for jet quenching from large to small systems, and should be helpful for finding the smallest QGP droplet and the disappearance of QGP in relativistic nuclear collisions.