This study investigated the initial mass function (IMF) and star formation history of high-mass stars in the Small Magellanic Cloud (SMC) using a population synthesis technique. We used the photometric survey catalog of Lee (2013) as the observable quantities and compare them with those of synthetic populations based on Bayesian inference. For the IMF slope (γ) range of -1.1 to -3.5 with steps of 0.1, five types of star formation models were tested: 1) continuous; 2) single burst at 10 Myr; 3) single burst at 60 Myr; 4) double bursts at those epochs; and 5) a complex hybrid model. In this study, a total of 125 models were tested. Based on the model calculations, it was found that the continuous model could simulate the high-mass stars of the SMC and that its IMF slope was -1.6 which is slightly steeper than Salpeter's IMF, i.e., γ=-1.35.