求一个数模p的逆元我们可以用费马小定理加快速幂求解或者用扩展欧几里得来求解,他们的复杂度都是o(log(n))的,都是很优秀的求法。但是当我们要求前n个数的逆元时就需要o(nlog(n))的复杂度了,如果我们可以通过第i个数前面的数的逆元来推出第i个数的逆元的话那么我们就可以线性的以o(n)的时间复杂度来求出前n个数的逆元了。证明如下。
设t=p/i,k=p%i
可知t*i+k≡0(mod p)
等式两边同时乘上inv(i)与inv(k)可得
t*inv(k)+inv(i)≡0(mod p)
所以
inv(i)=(p-p/i)*inv(p%i)%p
这个代码也很简单实际内容就一行公式。
#include<iostream>
#include<algorithm>
using namespace std;
const int N=20000999;
long long inv[N];
int main(){
long long p;
int n;
inv[1]=1;
scanf("%d%lld",&n,&p);
printf("%lld\n",inv[1]);
for(int i=2;i<=n;i++){
inv[i]=(p-p/i)*inv[p%i]%p;
printf("%lld\n",inv[i]);
}
return 0;
}