【数学 乘法逆元】luogu_3811 乘法逆元

本文介绍了一种求解模意义下乘法逆元的有效算法,通过递推公式加速计算过程,适用于大范围数值求解。文章详细解析了算法思路,并提供了完整的C++代码实现。

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

题意

给出n,pn,pn,p,求1∼n1\sim n1n中模ppp意义下的乘法逆元。

思路

假设已经求出[1,a−1][1,a-1][1,a1]的逆元,现在求aaa的逆元。

k=⌊pa⌋k=\left \lfloor \frac{p}{a} \right \rfloork=apr=pmod  ar=p\mod ar=pmoda,那么有
ka+r=pka+r=pka+r=p
所以
ka+r≡0ka+r\equiv 0ka+r0
(ka+r)a−1r−1≡0(ka+r)a^{-1}r^{-1}\equiv 0(ka+r)a1r10
kr−1+a−1≡0kr^{-1}+a^{-1}\equiv 0kr1+a10
a−1≡−kr−1a^{-1}\equiv -kr^{-1}a1kr1
因为−kr−1-kr^{-1}kr1是负数,所以我们让等式两边加上ppp
a−1+p≡−kr−1+pa^{-1}+p\equiv -kr^{-1}+pa1+pkr1+p
a−1≡p−⌊pa⌋r−1a^{-1}\equiv p-\left \lfloor \frac{p}{a} \right \rfloor r^{-1}a1papr1
由于已经得出r−1r^{-1}r1,所以可以直接推出a−1a^{-1}a1

代码

#include<cstdio>

int n, p;
int inv[3000001];

int main() {
	scanf("%d %d", &n, &p);
	inv[1] = 1;
	printf("%d\n", inv[1]);
	for (int i = 2; i <= n; i++) {
		inv[i] = p - (long long)(p / i) * inv[p % i] % p;
		printf("%d\n", inv[i]);
	}
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值