P3811 【模板】乘法逆元

题目

这是一道模板题

题目描述

给定 $n,p求1\sim n$中所有整数在模 $p$ 意义下的乘法逆元。

输入格式

一行两个正整数 $n,p$。

输出格式

输出 $n$ 行，第 $i$ 行表示 $i$ 在模 $p$ 下的乘法逆元。

输入输出样例

输入

10 13

输出

1
7
9
10
8
11
2
5
3
4

说明/提示

$1\le n\le 3\times 10^6,n<p<20000528$

输入保证 $p$ 为质数。

题解

模板乘法逆元
可把刚打好的扩欧乘法逆元板子交上去结果发现 T 了，
看了眼题解才发现要用线性方法求逆元，因为题目要求1~n的逆元。

线性法要推导一个递推公式：
设 $p=k\ast i+r$
也就是说 $p/i$ 的商为 $k$，余数为 $r$。
$\therefore k=\lfloor \frac{p}{i}\rfloor ,r=pmodi$

得 $k\ast i+r\equiv 0(modp)$

左右同乘 $(i^{-1}\ast r^{-1})$，              \;\;\;\;\;\;\;// $x^{-1}$ 表示 x 的逆元

得：$k\ast r^{-1}+i^{-1}\equiv 0(modp)$

代入$k=\lfloor \frac{p}{i}\rfloor ,r=pmodi$

得： $\lfloor \frac{p}{i}\rfloor \ast (pmodi{)}^{-1}+i^{-1}\equiv 0(modp)$

$\therefore i^{-1}\equiv -\lfloor \frac{p}{i}\rfloor \ast (pmodi{)}^{-1}(modp)$

$\therefore i^{-1}\%p=-\lfloor \frac{p}{i}\rfloor \ast (p\%i{)}^{-1}\%p$

$\therefore i^{-1}=-\lfloor \frac{p}{i}\rfloor \ast (p\%i{)}^{-1}\%p$ ( i 的逆元必小于 p，所以左边模 p可省掉)

所以我们便有了个递推式：

$inv[i]=-p/i\ast inv[p\%i]\%p$

$\because p\%p等于0$

$\therefore inv[i]=p-p/i\ast inv[p\%i]\%p加上p答案不变$

这样就能保证 $inv[i]>0$

超级简洁的代码

#include<iostream>
using namespace std;
long long p,n,inv[10000001];
int main()
{
	scanf("%d%d",&n,&p);
	inv[1]=1; //初始化，1在所有意义下的逆元都为1
	for(int i=2; i<=n; i++) inv[i]=(-p/i)*inv[p%i]%p; //从2开始套式子递推
	for(int i=1; i<=n; i++) printf("%d\n",inv[i]); //最后输出
	return 0;
}