乘法逆元

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

摘要生成于 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 \equiv1(mod p)

a/b*b*inv = 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 \equiv1(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;

 

评论 8
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值