p = 3660057339895840489386133099442699911046732928957592389841707990239494988668972633881890332850396642253648817739844121432749159024098337289268574006090698602263783482687565322890623 b = 1515231655397326550194746635613443276271228200149130229724363232017068662367771757907474495021697632810542820366098372870766155947779533427141016826904160784021630942035315049381147 x1 = 2157670468952062330453195482606118809236127827872293893648601570707609637499023981195730090033076249237356704253400517059411180554022652893726903447990650895219926989469443306189740 x2 = 1991876990606943816638852425122739062927245775025232944491452039354255349384430261036766896859410449488871048192397922549895939187691682643754284061389348874990018070631239671589727 k_2 = x2 + 2*x1 R.<a> = PolynomialRing(Zmod(p), implementation='NTL')#最好以后这么写 f = 4*k_2*(x1^3+a*x1+b)-(3*x1^2+a)^2 f = f.monic() x=f.roots() for each in x: m = int(each[0]) print(long_to_bytes(m))