线性求逆元

来自Miskcoo's Space
题目链接
要求线性内求出1-n中所有数的逆元\((n<=3000000)\)
这时候\(nlog(n)\)算法就有点悬了
所以需要一种\(O(n)\)的算法
这个方法是这样的
首先\(1^{-1} \equiv 1 \pmod p\)
然后我们设\(p = k\cdot i + r,~r < i,~1 < i < p\)
再将这个式子放到\(\mod p\)意义下就会得到
\[k\cdot i + r \equiv 0 \pmod p\]
两边同时乘上\(i^{-1}\cdot r^{-1}\)就会得到
\[k\cdot r^{-1}+i^{-1}\equiv 0\pmod p\]
\[i^{-1}\equiv -k\cdot r^{-1} \pmod p\]
\[i^{-1}\equiv -\left\lfloor\frac{p}{i}\right\rfloor\cdot \left(p\bmod i\right)^{-1} \pmod p\ \]

于是就可以从前面推出当前的逆元了,代码也就一行

A[i]=-(p/i)*A[p%i];
#include<iostream>
#include<cstdio>
#define N 3000005
using namespace std;

typedef long long ll;

ll _inv[N];  
void pre(int MOD){  
    _inv[0]=_inv[1]=1;  
    for(int i=2;i<N;i++)
        _inv[i]=((MOD-MOD/i)*_inv[MOD%i])%MOD;  
}    
int main(){  
    int n,p;
    cin>>n>>p;
    pre(p);
    for(int i=1;i<=n;++i)
        printf("%d\n",_inv[i]);
    return 0;  
}  

\[ k\cdot r^{-1} + i^{-1} &\equiv& 0 &\pmod p
i^{-1} &\equiv& -k\cdot r^{-1} &\pmod p
i^{-1} &\equiv& -\left\lfloor\frac{p}{i}\right\rfloor\cdot \left(p\bmod i\right)^{-1} &\pmod p \]

转载于:https://www.cnblogs.com/qdscwyy/p/7795368.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值