The question can be answered fairy simply. Main sequence stars are fusing hydrogen in their core via a mechanism that closely regulates the core temperature, so to within about a factor of 2, they all have the same core temperature-- and that's your question. But this does not require that they should have the same surface temperature. The logic is, the surface temperature is set by the luminosity and the radius of the star, via the formula L = kR2Ts4, for k a constant, which solves for the surface temperature Ts = (L/kR2)1/4. So this shows that you need to know R and L. The L is generally set by radiative diffusion in the interior of the star, which depends on mass essentially because the mass is the "stuff" the light has to diffuse through, and this leads to a relation like L is proportional to M raised to the power 3.5 or so, depending on some opacity details. The point is, M determines L. So now you only need R, and this depends on the history of contraction required to get the core to fusion temperature. That can be determined by using the fact that the average energy per particle necessary for fusion must be about equal to the potential energy per particle, set by M/R. So that sets R-- it is the R needed, given M, to get fusion in the core, so roughly R is proportional to M. Put it all together, and you find Ts is proportional to M to a power of about 0.4, roughly. That says a star 10 times more massive, or 1/10 as massive, than the Sun should have a surface T that is about 3 times higher, or 1/3 as high. That's about right.