码迷,mamicode.com
首页 > 其他好文 > 详细

线性求逆元推导

时间:2018-10-20 16:21:57      阅读:174      评论:0      收藏:0      [点我收藏+]

标签:const   acpi   介绍   int   scanf   code   max   线性   rac   

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


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

·首先,我们设\(t=P/i,k=P%i\);

·对于\(i*t+k≡0 \pmod{P}\),我们可以做出如下推导:

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

·从而得到:\(\frac{P}{i}*inv[P\%i]+inv[i]≡0 \pmod{P}\);

·最后得到\(inv[i]=(-\frac{P}{i}+P)*inv[P\%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]);
}

线性求逆元推导

标签:const   acpi   介绍   int   scanf   code   max   线性   rac   

原文地址:https://www.cnblogs.com/KatouKatou/p/9821642.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!