简单认识乘法逆元
如果
,且
(a与p互质),则称a关于模p的乘法逆元为x。
用代码不太恰当的举个栗子:
if(a*x%p==1){ cout<<x<<endl; }
- 费马小定理
由费马小定理 (p为素数),稍作变形即是 ,a^(p-2)即是a的逆元。 - 乘法逆元的作用:a除以一个数模p,等于a乘这个数的乘法逆元模p。
例题:P3811【模板】乘法逆元(luogu)
题目:给定n,p求1~n中所有整数在模p意义下的乘法逆元。
- 递推计算连续的数的逆元
#include<cstdio>
#define ll long long
const int maxn=3e6+5;
ll inv[maxn]= {0,1};
int main()
{
int n,p;
scanf("%d%d",&n,&p);
inv[1]=1;
printf("1\n");
for(int i=2; i<=n; i++)
inv[i]=(p-p/i)*inv[p%i]%p, printf("%d\n", inv[i]);
return 0;
}
- 解题也可以用费马小定理+快速幂,不过这道题卡时间,会TL。
#include <cstdio>
#define ll long long
using namespace std;
ll n, p;
long long PowMod(long long a,ll b)
{
long long ret=1;
while(b)
{
if(b&1)ret=ret*a%p;
a=a*a%p;
b>>=1;
}
return ret;
}
int main()
{
scanf("%d%d", &n, &p);
for(ll i=1; i<=n; i++)
{
printf("%d\n",PowMod(i, p-2));
}
return 0;
}