洛谷 P3811 [模板] 模意义下的乘法逆元

文章讲解了如何通过递推方式求解给定范围内的整数在质数p下的乘法逆元问题,给出了C++代码实例。

【模板】模意义下的乘法逆元

题目描述

给定 n,pn,pn,p1∼n1\sim n1n 中所有整数在模 ppp 意义下的乘法逆元。

这里 aaappp 的乘法逆元定义为 ax≡1(modp)ax\equiv1\pmod pax1(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 ^ 61n3×106n<p<20000528n < p < 20000528n<p<20000528

输入保证 ppp 为质数。

思路

递推式的起点:1−1≡1(modp)1^{-1}\equiv 1 \pmod{p}111(modp)i=1i=1i=1时)。
对于 i≥2i\geq2i2,设 p=k∗i+tp=k*i+tp=ki+t,则 k∗i+t≡0(modp)k*i+t\equiv 0 \pmod{p}ki+t0(modp)
然后同时乘上 i−1,t−1i^{-1},t^{-1}i1,t1,则可以得到 k∗t−1+i−1≡0(modp)k*t^{-1}+i^{-1}\equiv 0 \pmod{p}kt1+i10(modp)
移项可得 i−1≡−k∗t−1(modp)i^{-1}\equiv -k*t^{-1} \pmod{p}i1kt1(modp)
由于 k=⌊pi⌋k=\lfloor\frac{p}i\rfloork=ipt=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}i1ip(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}i1pip(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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值