We present a method based on mean field theory to cope with the mixed nonlinear integer programming, especially with optimal power flow problems involving both continuous and discrete variables, in a more exact manner. That is, we first formulate OPF as a mixed integer programming, and then derive its mean field equations as well as the annealing algorithm, by taking advantage of the characteristics of the original problems. Numerical simulations have verified effectiveness of this approach for small power systems.