代码：
p = 8637633767257008567099653486541091171320491509433615447539162437911244175885667806398411790524083553445158113502227745206205327690939504032994699902053229 q = 12640674973996472769176047937170883420927050821480010581593137135372473880595613737337630629752577346147039284030082593490776630572584959954205336880228469 dp = 6500795702216834621109042351193261530650043841056252930930949663358625016881832840728066026150264693076109354874099841380454881716097778307268116910582929 dq = 783472263673553449019532580386470672380574033551303889137911760438881683674556098098256795673512201963002175438762767516968043599582527539160811120550041 c = 24722305403887382073567316467649080662631552905960229399079107995602154418176056335800638887527614164073530437657085079676157350205351945222989351316076486573599576041978339872265925062764318536089007310270278526159678937431903862892400747915525118983959970607934142974736675784325993445942031372107342103852 from Crypto.Util.number import inverse import libnum #这里求的是p的逆元 I = inverse(p,q) m1 = pow(c,dp,p) m2= pow(c,dq,q) m = int((((m1-m2)*I)%p)*q+m2) print(libnum.n2s(m))
在python3.12 中可以执行以上代码，当然我们使用gmpy2也可以，执行效果一样，只是使用的库函数不同
代码解释
from Crypto.Util.number import inverse import libnum 导入模块 inverse用于计算逆元 libnum 用于下面输出函数print(libnum.n2s(m))将数字输出为字符串 I = inverse(p,q) #计算p的逆元 m1 = pow(c,dp,p) #计算m1 m2= pow(c,dq,q) #计算m2 m = int((((m1-m2)*I)%p)*q+m2) #使用int强制将输出内容转为整数int print(libnum.n2s(m)) #将整数转换为字符串 p = 8637633767257008567099653486541091171320491509433615447539162437911244175885667806398411790524083553445158113502227745206205327690939504032994699902053229 q = 12640674973996472769176047937170883420927050821480010581593137135372473880595613737337630629752577346147039284030082593490776630572584959954205336880228469 dp = 6500795702216834621109042351193261530650043841056252930930949663358625016881832840728066026150264693076109354874099841380454881716097778307268116910582929 dq = 783472263673553449019532580386470672380574033551303889137911760438881683674556098098256795673512201963002175438762767516968043599582527539160811120550041 c = 24722305403887382073567316467649080662631552905960229399079107995602154418176056335800638887527614164073530437657085079676157350205351945222989351316076486573599576041978339872265925062764318536089007310270278526159678937431903862892400747915525118983959970607934142974736675784325993445942031372107342103852 import gmpy2 import libnum \#这里求的是p的逆元 I = gmpy2.invert(p,q) m1 = pow(c,dp,p) m2= pow(c,dq,q) m = int((((m1-m2)*I)%p)*q+m2) print(libnum.n2s(m))
以上是在python3.10的执行效果
输出都是
