Three dimensional conformal field theories(CFTs) describe important critical physics in the real world. In the past few years we have made a lot of progress in 3D CFTs using bootstrap methods. In particular, using numeric bootstrap, critical exponents of the 3D Ising model, Super-Ising model, O(2) model, O(3) model have been determined to state-of-the-art precision. On the other hand, using analytic bootstrap, the information of large spin operators in the leading twist can be computed based on the numeric data of lower lying operators. The two methods are sensitive to different regions of the spectrum and complement each other. In this talk, I will first give a general review of recent numerical bootstrap results. Then I will discuss how to combine the numerics with the analytics. As a result, the bootstrap prediction can be greatly improved.