【模板】模意义下的乘法逆元
题目描述
给定 n,pn,pn,p 求 1∼n1\sim n1∼n 中所有整数在模 ppp 意义下的乘法逆元。
这里 aaa 模 ppp 的乘法逆元定义为 ax≡1(modp)ax\equiv1\pmod pax≡1(modp) 的解。
输入格式
一行两个正整数 n,pn,pn,p。
输出格式
输出 nnn 行,第 iii 行表示 iii 在模 ppp 下的乘法逆元。
样例 #1
样例输入 #1
10 13
样例输出 #1
1
7
9
10
8
11
2
5
3
4
提示
1≤n≤3×1061 \leq n \leq 3 \times 10 ^ 61≤n≤3×106,n<p<20000528n < p < 20000528n<p<20000528。
输入保证 ppp 为质数。
思路
递推式的起点:1−1≡1(modp)1^{-1}\equiv 1 \pmod{p}1−1≡1(modp)(i=1i=1i=1时)。
对于 i≥2i\geq2i≥2,设 p=k∗i+tp=k*i+tp=k∗i+t,则 k∗i+t≡0(modp)k*i+t\equiv 0 \pmod{p}k∗i+t≡0(modp),
然后同时乘上 i−1,t−1i^{-1},t^{-1}i−1,t−1,则可以得到 k∗t−1+i−1≡0(modp)k*t^{-1}+i^{-1}\equiv 0 \pmod{p}k∗t−1+i−1≡0(modp),
移项可得 i−1≡−k∗t−1(modp)i^{-1}\equiv -k*t^{-1} \pmod{p}i−1≡−k∗t−1(modp),
由于 k=⌊pi⌋k=\lfloor\frac{p}i\rfloork=⌊ip⌋,t=p%it=p\%it=p%i,
所以 i−1≡−⌊pi⌋∗(p%i)−1(modp)i^{-1}\equiv -\lfloor\frac{p}i\rfloor*(p\%i)^{-1} \pmod{p}i−1≡−⌊ip⌋∗(p%i)−1(modp),
即 i−1≡p−⌊pi⌋∗(p%i)−1(modp)i^{-1}\equiv p-\lfloor\frac{p}i\rfloor*(p\%i)^{-1} \pmod{p}i−1≡p−⌊ip⌋∗(p%i)−1(modp)。
代码
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 3e6 + 6;
ll inv[maxn];
int main()
{
ios::sync_with_stdio(0);
cin.tie(0);
cout.tie(0);
ll n, p;
cin >> n >> p;
inv[1] = 1; // 初始化inv[1]
cout << inv[1] << '\n';
for (ll i = 2; i <= n; i++)
{
inv[i] = p - p / i * inv[p % i] % p; // 递推求乘法逆元
cout << inv[i] << '\n';
}
return 0;
}
文章讲解了如何通过递推方式求解给定范围内的整数在质数p下的乘法逆元问题,给出了C++代码实例。

770

被折叠的 条评论
为什么被折叠?



