python - 大数的模取幂

Question

我目前正在处理大量密码学。

a = 3087646334 

p = 1606938044258990275541962092341162602522202993782792835301611

b = int((p-1)/2).

欧拉准则说

a^b (mod p) = 1 or a^b (mod p) = p-1.

SAGE 给出正确答案，python 给出错误答案。这是为什么？

预期输出 - 1 or p-1

SAGE 输出 - 1

Python 输出 - 1047939464127281862631850078334726680804120494559424004614779

代码在下面。
智者

a = 3087646334 
p = 1606938044258990275541962092341162602522202993782792835301611
b = int((p-1)/2)

power_mod(a,b,p)

Python

a = 3087646334
p = 1606938044258990275541962092341162602522202993782792835301611
b = int((p-1)/2)

pow(a,b,p)

Answer 1

因为您使用的是 float 除法，这是不准确的。使用 b = (p-1) // 2反而。