拿到题,发现给的e不常规,p1和p2相等,有两个不同n,两个不同c和两个不同e。给定两个密文的情况下,通常需要找到两者之间存在的关系,“合并”密文求解才能得到正确的明文。先计算phi1和phi2,然后分别计算e和phi的gcd,发现都等于14,所以这里问题考察e和phi不互素情况下的求解。
1 gcd(e,phi)=b
2 e = a * b
3 ed ≡ 1 mod phi
4 me ≡ c mod n
5 abd ≡ 1 mod phi 由2和3得到,这里可以推出abd=1
6 mab ≡ c mod n 由2和4得到,这里可以变换成c ≡ mab mod n
存在性质:
由以上可以推出cbd ≡ mabbd mod n 即得cbd ≡ mb mod n
所以我们先要求出bd,然后就能算出mb。求bd我们可以根据5式子,所以要算出a。由题目分析我们知道b=14,所以a=e/b。代码如下:
p = p1
phi1 = (p - 1) * (q1 - 1)
phi2 = (p - 1) * (q2 - 1)
b = gmpy2.gcd(e1, phi1)
a1 = e1 // b
a2 = e2 // b
bd1 = gmpy2.invert(a1, phi1)
bd2 = gmpy2.invert(a2, phi2)
由以上我们分别求出了bd,然后就可以分别算出mb。到了这一步还是无法算出m的,我第一次看别人题解时就想,直接调用gmpy2.iroot(mb,14)不就算出来m了吗?试了,结果不行。原因?欢迎你告诉我。
无法算出mb所以我们要再找新的解题思路。于是找到了同余式组,如下:
c1 ≡ mb1 mod n1 变换得到 c1 ≡ mb1 mod (p1 * q1)
c2 ≡ mb2 mod n2 同理 c2 ≡ mb2 mod (p2 * q2)
根据同余性质:
a ≡ b mod (mn) 变换成 a ≡ b mod m且a ≡ b mod n
所以得到下列式子:
c1 ≡ mb1 mod p1
c1 ≡ mb1 mod q1
c2 ≡ mb2 mod p2
c2 ≡ mb2 mod q2
因为p1和p2相等,根据如下性质
我们可以得到c3 ≡ mb1 * mb2 mod p
代码写法如下:
mb1 = pow(c1, bd1, n1)
mb2 = pow(c2, bd2, n2)
c3 = mb1 * mb2 % p
c2 = mb2 % q2
c1 = mb1 % q1
接下来要引用中国同余定理,所以先要有所了解,百度百科解释如下:
简单来说,就是求一个整数满足一组同余式组。所以这里把c1,c2和c3‘连接’起来,求一个特解c,即c满足如下:
c3 = c % p
c2 = c % q2
c1 = c % q1
所以求特解c的代码如下:
res = solve_crt([c1, c2, c3], [q1, q2, p])
于是我们重新得到了c ≡ me mod n的形式:
res ≡ m14 mod q1*q2
根据上面式子我们就可以求出真正的m14,也许你会想:之前不是都求出了mb了吗?一开始我也挺疑惑的,思考发现之前的mb都是分开的,只有求出特解c之后求的
mb才是真正的m14。根据以上的求解的思路,可以得到如下式子:
因为n = q1 * q2,所以我们只要对7和phi求逆元就可以求出2d,于是就可以求出m2,接着对m2开根就可以得到m。代码如下:
n = q1 * q2
f = (q1 - 1) * (q2 - 1)
m = res % n
d2 = gmpy2.invert(7, f)
m = pow(m, d2, n)
msg = gmpy2.iroot(m, 2)[0]
print(long_to_bytes(msg))
from Cryptodome.Util.number import *
import gmpy2
from libnum import *
e1 = 14606334023791426
p1 = 121009772735460235364940622989433807619211926015494087453674747614331295040063679722422298286549493698150690694965106103822315378461970129912436074962111424616439032849788953648286506433464358834178903821069564798378666159882090757625817745990230736982709059859613843100974349380542982235135982530318438330859
q1 = 130968576816900149996914427770826228884925960001279609559095138835900329492765336419489982304805369724685145941218640504262821549441728192761733409684831633194346504685627189375724517070780334885673563409259345291959439026700006694655545512308390416859315892447092639503318475587220630455745460309886030186593
c1 = 11402389955595766056824801105373550411371729054679429421548608725777586555536302409478824585455648944737304660137306241012321255955693234304201530700362069004620531537922710568821152217381257446478619320278993539785699090234418603086426252498046106436360959622415398647198014716351359752734123844386459925553497427680448633869522591650121047156082228109421246662020164222925272078687550896012363926358633323439494967417041681357707006545728719651494384317497942177993032739778398001952201667284323691607312819796036779374423837576479275454953999865750584684592993292347483309178232523897058253412878901324740104919248
n1 = p1 * q1
e2 = 13813369129257838
p2 = 121009772735460235364940622989433807619211926015494087453674747614331295040063679722422298286549493698150690694965106103822315378461970129912436074962111424616439032849788953648286506433464358834178903821069564798378666159882090757625817745990230736982709059859613843100974349380542982235135982530318438330859
q2 = 94582257784130735233174402362819395926641026753071039760251190444144495369829487705195913337502962816079184062352678128843179586054535283861793827497892600954650126991213176547276006780610945133603745974181504975165082485845571788686928859549252522952174376071500707863379238688200493621993937563296490615649
c2 = 7984888899827615209197324489527982755561403577403539988687419233579203660429542197972867526015619223510964699107198708420785278262082902359114040327940253582108364104049849773108799812000586446829979564395322118616382603675257162995702363051699403525169767736410365076696890117813211614468971386159587698853722658492385717150691206731593509168262529568464496911821756352254486299361607604338523750318977620039669792468240086472218586697386948479265417452517073901655900118259488507311321060895347770921790483894095085039802955700146474474606794444308825840221205073230671387989412399673375520605000270180367035526919
n2 = p2 * q2
p = p1
phi1 = (p - 1) * (q1 - 1)
phi2 = (p - 1) * (q2 - 1)
b = gmpy2.gcd(e1, phi1)
a1 = e1 // b
a2 = e2 // b
bd1 = gmpy2.invert(a1, phi1)
bd2 = gmpy2.invert(a2, phi2)
mb1 = pow(c1, bd1, n1)
mb2 = pow(c2, bd2, n2)
c3 = mb1 * mb2 % p
c2 = mb2 % q2
c1 = mb1 % q1
res = solve_crt([c1, c2, c3], [q1, q2, p]) #crt中国同余定理,求出特解作为
print(res)
n = q1 * q2
f = (q1 - 1) * (q2 - 1)
m = res % n
d2 = gmpy2.invert(7, f)
m = pow(m, d2, n)
msg = gmpy2.iroot(m, 2)[0]
print(long_to_bytes(msg))
得到flag :
中间的\xcf\x86为utf-8编码,转换得到字符:φ
最终flag为:flag{gcd_e&φ_isn't_1}
这题真的真的真的花费了我好多时间去理解,到目前我依然存在一些不理解的地方。当然我也从中学到了一些东西,比如CRT,e和phi不互素时的求解思路等。一开始我是不会做这道题的,于是我去寻找wp,发现每个wp中总是说都不够透彻,所以看了wp总觉得会了但又没完全会,于是写了这篇博客分享我的理解,随便整理思路过程。如果你觉得有所收获,可以跟我分享一下你对这题的理解(也许我理解在某个点不对呢?),这也许能帮助到我。
数论-模运算与同余的性质_AcmAlgorithm的博客-CSDN博客_模同余关系性质的证明
手机扫一扫
移动阅读更方便
你可能感兴趣的文章