Vaccination is effective in preventing human papillomavirus (HPV) infection. It still remains
debatable whether males should be included in a vaccination program and unclear how to allocate
the vaccine in genders to achieve the maximum benefits. In this paper, we use a two-sex model
to assess HPV vaccination strategies and use the data from Guangxi Province in China as a case
study. Both mathematical analysis and numerical simulations show that the basic reproduction
number, an important indicator of the transmission potential of the infection, achieves its minimum
when the priority of vaccination is given to the gender with a smaller recruit rate. Given a fixed
amount of vaccine, splitting the vaccine evenly usually leads to a larger basic reproduction number
and a higher prevalence of infection. Vaccination becomes less effective in reducing the infection
once the vaccine amount exceeds the smaller recruit rate of the two genders. In the case study,
we estimate the basic reproduction number is 1.0333 for HPV 16/18 in people aged 15-55. The
minimal bivalent HPV vaccine needed for the disease prevalence to be below 0.05% is 24050 per
year, which should be given to females. However, with this vaccination strategy it would require
a very long time and a large amount of vaccine to achieve the goal. In contrast with allocating the
same vaccine amount every year, we find that a variable vaccination strategy with more vaccine
given in the beginning followed by less vaccine in later years can save time and total vaccine
amount. The variable vaccination strategy illustrated in this study can help to better distribute the
vaccine to reduce the HPV prevalence. Although this work is for HPV infection and the case study
is for a province in China, the model, analysis and conclusions may be applicable to other sexually
transmitted diseases in other regions or countries.