题目;
from Crypto.Util.number import *
from flag import flag
def nextPrime(n):
n += 2 if n & 1 else 1
while not isPrime(n):
n += 2
return n
p = getPrime(1024)
q = nextPrime(p)
n = p * q
e = 0x10001
d = inverse(e, (p-1) * (q-1))
c = pow(bytes_to_long(flag.encode()), e, n)
# d = 19275778946037899718035455438175509175723911466127462154506916564101519923603308900331427601983476886255849200332374081996442976307058597390881168155862238533018621944733299208108185814179466844504468163200369996564265921022888670062554504758512453217434777820468049494313818291727050400752551716550403647148197148884408264686846693842118387217753516963449753809860354047619256787869400297858568139700396567519469825398575103885487624463424429913017729585620877168171603444111464692841379661112075123399343270610272287865200880398193573260848268633461983435015031227070217852728240847398084414687146397303110709214913
# c = 5382723168073828110696168558294206681757991149022777821127563301413483223874527233300721180839298617076705685041174247415826157096583055069337393987892262764211225227035880754417457056723909135525244957935906902665679777101130111392780237502928656225705262431431953003520093932924375902111280077255205118217436744112064069429678632923259898627997145803892753989255615273140300021040654505901442787810653626524305706316663169341797205752938755590056568986738227803487467274114398257187962140796551136220532809687606867385639367743705527511680719955380746377631156468689844150878381460560990755652899449340045313521804
首先,我们审计题目,发现我们已知的条件有d,c和e。同时,我们知道p和q是两个相邻的质数。于是,根据以上条件,我们开始解题;
一、首先,我们必须知道,想要解出此题。我们必须知道n的值,而想要求得n的值,则必须要求p和q的值。(n=p*q)
二、rsa基础公式:
e*d = 1 mod phi_n (其中phi_n = (p-1)*(q-1)=1 + k*phi_n
通过对e*d的计算,我们知道其长度为2066位,而生成p的条件为 getPrime(1024),所以(p-1)*(q-1)应该为2048位,2066-2048=16,所以k值必在 pow(2,15)至pow(2,16)之间。于是,我们可以根据此条件暴力求解k值,从而求出(p-1)*(q-1)的值。
三、当我们求得phi_n的值后,我们已知p和q是两个相邻的质数,所以,我们可以使用sympy库对p q两值进行求解。思路为先对(p-1)*(q-1)开方,再求得大于开方所得数和小于开方所得数的质数。
p = sympy.prevprime(gmpy2.iroot((e*d-1)//i,2)[0])
q = sympy.nextprime(p)
#sympy.prevprime(x)是求大于x最近的质数,sympy.nextprime(x)是求小于x最近的质数。
写出脚本对题目求解:
from Crypto.Util.number import long_to_bytes
import gmpy2
import sympy
d = 19275778946037899718035455438175509175723911466127462154506916564101519923603308900331427601983476886255849200332374081996442976307058597390881168155862238533018621944733299208108185814179466844504468163200369996564265921022888670062554504758512453217434777820468049494313818291727050400752551716550403647148197148884408264686846693842118387217753516963449753809860354047619256787869400297858568139700396567519469825398575103885487624463424429913017729585620877168171603444111464692841379661112075123399343270610272287865200880398193573260848268633461983435015031227070217852728240847398084414687146397303110709214913
c = 5382723168073828110696168558294206681757991149022777821127563301413483223874527233300721180839298617076705685041174247415826157096583055069337393987892262764211225227035880754417457056723909135525244957935906902665679777101130111392780237502928656225705262431431953003520093932924375902111280077255205118217436744112064069429678632923259898627997145803892753989255615273140300021040654505901442787810653626524305706316663169341797205752938755590056568986738227803487467274114398257187962140796551136220532809687606867385639367743705527511680719955380746377631156468689844150878381460560990755652899449340045313521804
e=65537
p=0
q=0
for i in range(pow(2,15),pow(2,16)):
if (e*d-1)%i==0:
p = sympy.prevprime(gmpy2.iroot((e*d-1)//i,2)[0])
q = sympy.nextprime(p)
if (p-1)*(q-1)==(e*d-1)//i:
break
n=p*q
m=gmpy2.powmod(c,d,n)
m1=long_to_bytes(m)
print(m1)
解得:
b'NCTF{70u2_nn47h_14_v3ry_gOO0000000d}'