转载:线性求逆元推导

转载自:线性求逆元推导 - Katoumegumi - 博客园 (cnblogs.com)

本篇介绍线性求逆元的推导过程

  • 对于一个质数 \(P\),我们需要求出 \(1-N\)\(\mod P\) 意义下的逆元,如何使用线性的方法求其逆元呢?

  • 首先,我们设 \(t=\lfloor\frac{P}{i}\rfloor,k=P\mod i\);

  • 对于 \(i\times t+k\equiv0 \pmod{P}\),我们可以做出如下推导:

  • 等式两边同时除以 \(i\times k\),我们可以得到新式子 \(\frac{t}{k}+\frac{1}{i}\equiv0 \pmod{P}\);

  • 从而得到: \(\lfloor\frac{P}{i}\rfloor\times inv[P\mod i]+inv[i]≡0 \pmod{P}\);

  • 最后得到 \(inv[i]=(-\lfloor\frac{P}{i}\rfloor+P)\times inv[P\mod i]%P\);

\(code:\)

#include<stdio.h>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=(1e7*2)+2;
ll n,p,inv[maxn];
inline ll add(ll a,ll b){return a+b<p?a+b:a+b-p;}
inline ll mul(ll a,ll b){return a*b<p?a*b:a*b%p;}
int main()
{
scanf("%lld%lld",&n,&p);inv[1]=1;
for(int i=2;i<=n;i++) inv[i]=mul(add(-p/i,p),inv[p%i]);
}
posted @   ZnPdCo  阅读(30)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 使用C#创建一个MCP客户端
· ollama系列1:轻松3步本地部署deepseek,普通电脑可用
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· 按钮权限的设计及实现
点击右上角即可分享
微信分享提示