pow(a, p - 2, p) pow(a, -1, p) できるpythondef mod_inverse_ee(a, m):
"""
Solve ax mod m = 1 with extended euclidean.
x = a^-1.
"""
x, y, g = extended_euclidean(a, m)
assert g == 1
return x % m>nのmodulo pでの逆元はフェルマーの小定理より
> pow(n, p-2, p)
> で求められます。計算量はO(logp)です。二項係数とかでよく使います。