A mathematical model for the natural history of human papillomavirus (HPV) is designed and used to assess the impact of a hypothetical anti-HPV vaccine and Pap cytology screening on the transmission dynamics of HPV in a population. Rigorous qualitative analysis of the model reveals that it undergoes the phenomenon of backward bifurcation. It is shown that the backward bifurcation is caused by the imperfect nature of the HPV vaccine or the HPV-induced and cancer-induced mortality in females. For the case when the disease-induced and cancer-induced mortality is negligible, it is shown that the disease-free equilibrium (i.e., equilibrium in the absence of HPV and associated dysplasia) is globally-asymptotically stable if the associated reproduction number is less than unity. The model has a unique endemic equilibrium when the reproduction threshold exceeds unity. The unique endemic equilibrium is globally-asymptotically stable for a special case, where the associated HPV-induced and cancer-induced mortality is negligible. Numerical simulations of the model, using a reasonable set of parameter values, support the recent recommendations by some medical agencies and organizations in the USA to offer Pap screening on a 3-year basis (rather than annually). (rather than annually).