一、若a与n互素,那么可以用扩展欧几里德和欧拉函数求出a对于n的逆元。
ax≡1(mod n),x为a对于n的逆元。
直接用欧几里德求逆元x:
//a与n互素,逆元才有解 #include <iostream> #include <cstdio> using namespace std; typedef long long LL; void exgcd(LL a,LL b,LL &d,LL &x,LL &y) { if(!b) { d=a;x=1;y=0; return ; } exgcd(b,a%b,d,y,x); y-=x*(a/b); } LL inv(LL a,LL n) { LL d,x,y; exgcd(a,n,d,x,y); return d==1?(x+n)%n:-1; } int main() { cout<<inv(2,1e9+7)<<'\n'; return 0; }
欧拉定理:对于任何两个互质的正整数a,n(n>=2),有a^(φ(n))≡1(mod
n)。
那么只要a与n互质,逆元x=a^(φ(n)-1)%n,可以用快速幂求出。
用欧拉定理求解逆元x:
LL inv(LL a,LL euler,LL n) { LL ret=1; --euler; while(euler) { if(euler&1) ret=ret*a%n; euler>>=1; a=a*a%n; } return ret; }
特殊地,若n是素数的话可以用费马小定理求解。
费马小定理:当n是素数时,a^(n-1)≡1(mod
n)。(费马小定理是欧拉定理的特殊情况)
那么逆元x=a^(n-2)%n。也可以用快速幂直接求出。
二、若a与n不互素,可以用通用的方法求出。(此方法求出的应该不能再叫逆元了)
引入问题:求p=a/b%n。(b与n不互素)
已知:除以一个数b等价与乘以它的乘法逆元inv(b)。
但是此时b与n不互素,因此不能用扩展欧几里德和欧拉定理求解。
可以直接用公式:p=a%(b*n)/b;公式不难证明。此方法是在poj1845的讨论区看到的........
另外附上大神的文章:ACdreamer-逆元详解
时间: 2024-10-12 17:32:56