sagemath matlab,使用SageMath

A = 0xd9d8bd32fc078af28cb318a7ae07227accf335a1a02eddff421a50b7eceb9fd7

B = 0xfffffffeffffffffffffffffffffffff000000000000000000000001

M = 0xffffffff00000001000000000000000000000000ffffffffffffffffffffffff

x = 0xfa00a8d4d5755ae92322a317c3b1895d935ed95a36b4e4589328d7b7214a5d35

y = 0xaeeb0346ec45386861311d68619391181116843e51e91c7d80a17faa47792fe8

x1 = 0x739ca3c42d60579b8ab304fe511eb7501343d3a951a85ec67e74cfb88fdde40b

y1 = 0x3bcc7adbf46e9ce70453644e27b6947e7382c23c315f5ea7847f6cb4f10b62a2

P = (x, y)

Q = (x1, y1)

F = FiniteField(M)

E = EllipticCurve(F,[A,B])

P = E.point(P)

Q = E.point(Q)

print factor(E.order())

factors, exponents = zip(*factor(E.order()))

primes = [2, 4, 2137, 7649, 63127, 33636409, 55866443, 194282080159,269685708029]

dlogs = []

for fac in primes:

t = int(P.order()) / int(fac)

dlog =discrete_log(t*Q,t*P,operation="+")

dlogs += [dlog]

print("factor:"+str(fac)+", Discrete Log: "+str(dlog))

#The result after run

l = crt(dlogs,primes)

print(l)

print P*l == Q