拿到题目，有两个文件，里面分别有n,e,c 看了一下两份文件的e相同，n不同 题目的名字是RSA_gcd，自然想到了求 n 1 n_1 n1, n 2 n_2 n2的最大公约数 求出的最大公约数若为素数的话，那么它就是p，题目自然就解出来了
import gmpy2 def get_p(n1, n2): p = gmpy2.gcd(n1, n2) assert gmpy2.is_prime(p) == True return p def int_str(data): flag = '' while data: flag += chr(data & 0xff) data >>= 8 return flag[::-1] if __name__ == "__main__": n1 = 23220619839642624127208804329329079289273497927351564011985292026254914394833691542552890810511751239656361686073628273309390314881604580204429708461587512500636158161303419916259271078173864800267063540526943181173708108324471815782985626723198144643256432774984884880698594364583949485749575467318173034467846143380574145455195152793742611717169602237969286580028662721065495380192815175057945420182742366791661416822623915523868590710387635935179876275147056396018527260488459333051132720558953142984038635223793992651637708150494964785475065404568844039983381403909341302098773533325080910057845573898984314246089 n2 = 22642739016943309717184794898017950186520467348317322177556419830195164079827782890660385734113396507640392461790899249329899658620250506845740531699023854206947331021605746078358967885852989786535093914459120629747240179425838485974008209140597947135295304382318570454491064938082423309363452665886141604328435366646426917928023608108470382196753292656828513681562077468846105122812084765257799070754405638149508107463233633350462138751758913036373169668828888213323429656344812014480962916088695910177763839393954730732312224100718431146133548897031060554005592930347226526561939922660855047026581292571487960929911 c1 = 9700614748413503291260966231863562117502096284616216707445276355274869086619796527618473213422509996843430296526594113572675840559345077344419098900818709577642324900405582499683604786981144099878021784567540654040833912063141709913653416394888766281465200682852378794478801329251224801006820925858507273130504236563822120838520746270280731121442839412258397191963036040553539697846535038841541209050503061001070909725806574230090246041891486506980939294245537252610944799573920844235221096956391095716111629998594075762507345430945523492775915790828078000453705320783486744734994213028476446922815870053311973844961 c2 = 20513108670823938405207629835395350087127287494963553421797351726233221750526355985253069487753150978011340115173042210284965521215128799369083065796356395285905154260709263197195828765397189267866348946188652752076472172155755940282615212228370367042435203584159326078238921502151083768908742480756781277358357734545694917591921150127540286087770229112383605858821811640935475859936319249757754722093551370392083736485637225052738864742947137890363135709796410008845576985297696922681043649916650599349320818901512835007050425460872675857974069927846620905981374869166202896905600343223640296138423898703137236463508 e = 65537 p = get_p(n1, n2) q1 = n1 // p q2 = n2 // p d1 = gmpy2.invert(e, (p - 1) * (q1 - 1)) d2 = gmpy2.invert(e, (p - 1) * (q2 - 1)) m1 = gmpy2.powmod(c1, d1, n1) m2 = gmpy2.powmod(c2, d2, n2) print(int_str(m1) + int_str(m2))