當運算時需要求模的時候,可以直接做的有+-*但是不滿足/
而逆元就是通過某種運算來達到求(a/b)%p的結果
當b*c≡1(mod p)時,
有(a/b)%p=(a/b*b*c)%p=(a*c)%p
這裏c就是b關於p的逆元
那麼如何求c呢
費馬小定理
當p爲素數時,有 a^p≡a(mod p)
故有a^p≡1(mod p)
因此a關於p的逆元就是a^(p-2)
這個冪可以使用快速冪來得出
擴展歐幾里得算法
對於不完全爲0的非負整數a,b,必存在整數x,y,滿足ax+by=gcd(a,b)
因此,求a*x≡1(mod p)
即求a*x+p*y=1=gcd(a,b)/gcd(a,b)
也就是先求出a*x+p*y=gcd(a,b)
然後x=x/gcd(a,b)即可
注意:x可能爲負的,需要x=x+p
typedef long long ll;
ll exgcd(ll a, ll b, ll &x, ll &y) //這是擴展歐幾里得
{
if(b == 0){
x = 1;
y = 0;
return a;
}
ll t = exgcd(b, a%b, x, y);
ll tmp = y;
y = x - a / b * y;
x = tmp;
return t;
}
k = exgcd(k ,p ,x ,i);
if(x < 0) x += p; //防止x爲負
x = x / k;