The immunophenoscore (IPS) is an important indicator for evaluating immunotherapy response. This work was designed to establish a prognostic model based on IPS-related genes in cervical cancer. Weighted correlation network analysis (WGCNA) was utilized to identify key modules related to IPS in cervical cancer data from The Cancer Genome Atlas (TCGA). The results show that the yellow module (158 genes) had a high correlation with both IPS_CTLA4_blocker and IPS_CTLA4_and PC1/PDL1/PDL2 blocker. Univariate cox regression analysis and LASSO regression analysis were performed based on 158 genes, and 9 characteristic genes were finally identified to construct the model. According to the differentially expressed genes, cervical cancer samples were divided into high-risk and low-risk groups and cluster 1.2.3. Higher risk scores associated with poorer prognosis. cluster2 and cluster3 were identified as cervical cancer subtypes with significant survival differences. cluster2 had higher immune cell infiltration levels and better prognosis, with greater sensitivity to Cyclopamine, Imatinib, MG-13, Paclitaxel, PHA-665752, Rapamycin, Sorafenib, Sunitinib, and VX-680. In contrast, cluster3 had higher TTN and PIK3CA mutations and greater sensitivity to AZ628, Dasatinib, Doxorubicin, HG-6-64-1, JQ12, Midostaurin, PF-562271, TAE684, and WH-4-023. In conclusion, we developed a feasible risk score model based on IPS-related genes for cervical cancer prognosis and identified potential drugs for different cervical cancer subtypes.