算法竞赛—乘法逆元

乘法逆元

定义：对于正整数$a$，存在至少一个$x$使得$ax\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$ ，则称所有满足以上条件的$x$均为$a$在模$p$意义下的乘法逆元，一般记作$a^{-1}$。

1. 用扩展欧几里得算法求逆元
   本质是求线性同余方程$ax\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$的解$x$。
   上式可变形为$ax+py=1(x,y\in N)$，由裴蜀定理可知，$gcd(a,p)=1$。
   扩展欧几里得算法即是基于以上等式诞生的。

int exgcd(int a,int b,int &x,int &y){
	if(b==0){
		x=1,y=0;
		return a;
	}
	int d=exgcd(b,a%b,y,x);
	y-=a/b*x;
	return d;
}

2. $p$为质数时，利用费马小定理简化求逆元
   费马小定理：当$p$为质数时，线性同余方程$ax\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$ 必存在解$x=a^{p-2}$。
   因此在此特殊条件下，求乘法逆元即为求$a^{p-2}\mathrm{m}\mathrm{o}\mathrm{d}p$，而这个高次幂可以简单地通过循环求出，或者用快速幂加速得出。故此代码省略。