In this talk, I will present recent developments in the thermal bootstrap with applications to the O(N) model, based on 2312.13030, 2411.00978, and ongoing work. I will describe how analytical and numerical methods can be combined to formulate a bootstrap approach for computing thermal OPE coefficients (i.e., thermal one-point functions). This approach uses zero-temperature conformal data as input, with the Kubo-Martin-Schwinger (KMS) condition serving as the key consistency constraint. Special emphasis will be given to applications in the 3D Ising, O(2), and O(3) models. Additionally, I will discuss preliminary results on the perturbative thermal bootstrap.