乘法逆元

最新推荐文章于 2023-08-09 13:52:48 发布

FSYo

最新推荐文章于 2023-08-09 13:52:48 发布

阅读量5.6k

点赞数 5

CC 4.0 BY-SA版权

分类专栏：乘法逆元文章标签：乘法逆元

本文链接：https://blog.youkuaiyun.com/sslz_fsy/article/details/81986892

乘法逆元专栏收录该内容

3 篇文章

订阅专栏

本文探讨了模运算的基本规则及其在大数运算中的应用，重点介绍了乘法逆元的概念及其实现方法，包括扩展欧几里得算法、欧拉定理和递推公式，为解决数学和计算机科学中的模运算问题提供了有效途径。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

若 $\tiny ax\equiv 1(mod p)$ ,gcd(a,b)=1

则称x为a模p的乘法逆元

我们先来看看有什么用

当输出结果很大时,要模一个mod在输出

(a+b)%mod=a%mod+b%mod

(a-b)%mod=a%mod-b%mod

(a*b)%mod=a%mod*b%mod

(a/b)%mod ... 呃

乘法逆元派上用场了,设b模p的乘法逆元为inv

(a/b)%p=(a*inv)%p=a%p*inv%p

为什么呢

因为b*inv $\tiny \equiv$ 1(mod p)

a/bbinv = a*inv $\tiny \equiv$ a/b(mod p)

举个例子

11x $\tiny \equiv$ 1 (mod 7) x=9

(110/11)%7=(110%7)*(9%7)=5*2%7=3

好了,想想怎么求吧

扩展欧几里得

ax $\tiny \equiv$ 1(mod p) 设ax=yp+1,b=-p 则ax+by=1

void exgcd(int a,int b)
{
	if(b==0) g=a,x=1,y=0;
	else{
		exgcd(b,a%b);
		int x1=y,y1=x-a/b*y;
		x=x1,y=y1;
		printf("%d*(%d)+%d*(%d)=%d\n",a,x,b,y,g);
	}
}
inv[a]=(x+p)%p;

复杂度 O(nlogn)

欧拉定理

欧拉函数

a ^ $\tiny \varphi(p)$ $\tiny \equiv$ 1(mod p)

a * a ^ $\tiny \varphi(p)-1$ $\tiny \equiv$ 1 (mod p)

所以a模p的乘法逆元是a ^ $\tiny \varphi(p)-1$

复杂度欧拉函数O(n),快速幂O(logn),O(nlogn)

void init(){//求欧拉函数(phi[i])
	for(int i=2;i<=n;i++){
		if(!p[i]){//是质数 
			p[i]=i,prime[++cnt]=i;//p->最小质因子
			phi[i]=i-1; //phi(prime)=prime-1
		}
		for(int j=1;j<=cnt;j++){
			if(prime[j]>p[i]||prime[j]*i>n) break;//不是最小质因子,或大于n
			p[prime[j]*i]=prime[j];
			phi[prime[j]*i]=phi[i]*(i%prime[j]?prime[j]-1:prime[j]);//欧拉函数的性质
		}
	}
}

递推

for(int i=2;i<=p-1;i++)
    inv[i]=(p-p/i)*inv[p%i]%p;

乘法逆元

若,gcd(a,b)=1

因为b*inv1(mod p)

a/b*b*inv = a*inv a/b(mod p)

扩展欧几里得

欧拉定理

递推

若 $\tiny ax\equiv 1(mod p)$ ,gcd(a,b)=1

因为b*inv $\tiny \equiv$ 1(mod p)

a/bbinv = a*inv $\tiny \equiv$ a/b(mod p)