This study developed an evolutionary algorithm, called the telescopic evolutionary algorithm (TEA), for the inversion of the spectral ratio of the horizontal to vertical components of strong motions (eHVSR), based on the diffuse-field concept. We verified the TEA using synthetic eHVSR curves from the literature and examined its convergence stability with respect to the independent variables of the algorithm. Furthermore, we present the application of TEA in the inversion of eHVSR curves for the shear wave velocity of the subsurface structure down to the seismic bedrock at seven stations of the K-NET/KiK-net strong-motion networks in Japan. The corre-sponding geological maps of the regions from the Geological Survey of Japan as well as the site condition database of K-NET/KiK-net and the Japan Seismic Hazard Information Station (J-SHIS) were processed. Subse-quently, the TEA retrieved the Vs profiles of the subsurface structures at these stations to highlight the perfor-mance of TEA in various geological settings. Finally, important site proxies were calculated from the retrieved subsurface structures.