c = 5382723168073828110696168558294206681757991149022777821127563301413483223874527233300721180839298617076705685041174247415826157096583055069337393987892262764211225227035880754417457056723909135525244957935906902665679777101130111392780237502928656225705262431431953003520093932924375902111280077255205118217436744112064069429678632923259898627997145803892753989255615273140300021040654505901442787810653626524305706316663169341797205752938755590056568986738227803487467274114398257187962140796551136220532809687606867385639367743705527511680719955380746377631156468689844150878381460560990755652899449340045313521804
d = 19275778946037899718035455438175509175723911466127462154506916564101519923603308900331427601983476886255849200332374081996442976307058597390881168155862238533018621944733299208108185814179466844504468163200369996564265921022888670062554504758512453217434777820468049494313818291727050400752551716550403647148197148884408264686846693842118387217753516963449753809860354047619256787869400297858568139700396567519469825398575103885487624463424429913017729585620877168171603444111464692841379661112075123399343270610272287865200880398193573260848268633461983435015031227070217852728240847398084414687146397303110709214913
e = 65537
from gmpy2 import *
from Crypto.Util.number import sieve_base
def solve_p_q():
# step 1
x = 0
k = ed - 1
# step 2
for i in range(len(sieve_base)):
g = sieve_base[i]
t = k
# step 3
if t % 2 == 0:
t = t // 2
x = powmod(g,t,n)
else:
continue
# step 4
y = gcd(x-1,n)
if 1<y<n:
p = y
q = n // p
return p,q
p,q = solve_p_q()
禾之舟 
