This paper proposed an efficient and adaptive frequency sampling algorithm for frequency response analysis using dynamic condensation-based reduced-order modeling. For the degree of freedom-based model reduction method, the reduced-order basis becomes a frequency-dependent matrix since the relationship between master and slave degrees of freedom stems from partial equations of a second-order dynamical system. Such frequency-dependency makes the analysis inefficient for investigating the frequency response of the system. Considering that the coverage of a local reduced-order basis at a single frequency varies depending on the frequency, a new frequency sampling algorithm was proposed with a strategy of constructing multiple local reduced-order models (ROMs) at sample frequencies. For adaptive sampling, the frequency range of a local ROM was evaluated, and a new sample was added if there was a gap between two adjacent ROMs. As a result, the accuracy of the local ROM can be estimated, and the efficiency in the online stage was greatly enhanced. The proposed method was verified by performing frequency response analysis with several numerical examples, including a large-scale structural and dynamic system.