We mapped the kinetic temperature structure of two massive star-forming regions, N113 and N159W, in the Large Magellanic Cloud (LMC). We have used $sim$1hbox{$,.!!^{primeprime}$}6,($sim$0.4,pc) resolution measurements of the para-H$_2$CO,$J_{rm K_ aK_c}$,=,3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$ transitions near 218.5,GHz to constrain RADEX non-LTE models of the physical conditions. The gas kinetic temperatures derived from the para-H$_2$CO line ratios 3$_{22}$--2$_{21}$/3$_{03}$--2$_{02}$ and 3$_{21}$--2$_{20}$/3$_{03}$--2$_{02}$ range from 28 to 105,K in N113 and 29 to 68,K in N159W. Distributions of the dense gas traced by para-H$_2$CO agree with those of the 1.3,mm dust and emph{Spitzer},8.0,$mu$m emission, but do not significantly correlate with the H$alpha$ emission. The high kinetic temperatures ($T_{rm kin}$,$gtrsim$,50,K) of the dense gas traced by para-H$_2$CO appear to be correlated with the embedded infrared sources inside the clouds and/or YSOs in the N113 and N159W regions. The lower temperatures ($T_{rm kin}$,$<$,50,K) are measured at the outskirts of the H$_2$CO-bearing distributions of both N113 and N159W. It seems that the kinetic temperatures of the dense gas traced by para-H$_2$CO are weakly affected by the external sources of the H$alpha$ emission. The non-thermal velocity dispersions of para-H$_2$CO are well correlated with the gas kinetic temperatures in the N113 region, implying that the higher kinetic temperature traced by para-H$_2$CO is related to turbulence on a $sim$0.4,pc scale. The dense gas heating appears to be dominated by internal star formation activity, radiation, and/or turbulence. It seems that the mechanism heating the dense gas of the star-forming regions in the LMC is consistent with that in Galactic massive star-forming regions located in the Galactic plane.