P3811-[模板]乘法逆元【线性求逆元】
正題
評測記錄:https://www.luogu.org/recordnew/lists?uid=52918&pid=P3811
題目大意
求1~n1\sim n1~n中modpmod\ \ pmod??p的逆元。
解題思路
首先對于p,我們將其分解為ki+r(k=?pi?,r=p%r)ki+r(k=\lfloor \frac{p}{i}\rfloor,r=p\%r)ki+r(k=?ip??,r=p%r),然后有
ki+r≡0(modp)ki+r\equiv 0(mod\ \ p)ki+r≡0(mod??p)
左右兩邊同時乘上一個i?1?r?1i^{-1}*r^{-1}i?1?r?1
k?r?1+i?1≡0(modp)k*r^{-1}+i^{-1}\equiv 0(mod\ \ p)k?r?1+i?1≡0(mod??p)
i?1≡?k?r?1(modp)i^{-1}\equiv -k*r^{-1}(mod\ \ p)i?1≡?k?r?1(mod??p)
i?1≡??pi??(p%i)?1(modp)i^{-1}\equiv -\lfloor \frac{p}{i}\rfloor*(p\%i)^{-1}(mod\ \ p)i?1≡??ip???(p%i)?1(mod??p)
反正本來就要%p\%p%p
i?1=??pi??(p%i)?1%pi^{-1}= -\lfloor \frac{p}{i}\rfloor*(p\%i)^{-1}\%pi?1=??ip???(p%i)?1%p
因為iii比(p%i)(p\%i)(p%i)小,按照遞推的順序我們在求出i?1i^{-1}i?1之前就已經求出(p%i)?1(p\%i)^{-1}(p%i)?1了。但是我們還要保證不是負數,所以我們可以直接計算
i?1=p??pi??(p%i)?1%pi^{-1}= p-\lfloor \frac{p}{i}\rfloor*(p\%i)^{-1}\%pi?1=p??ip???(p%i)?1%p
code
#include<cstdio> using namespace std; int n,p; long long inv[3000010]; int main() {scanf("%d%d",&n,&p);inv[1]=1;printf("%d\n",inv[1]);for(int i=2;i<=n;i++)inv[i]=(long long)p-p/i*inv[p%i]%p,printf("%d\n",inv[i]); }總結
以上是生活随笔為你收集整理的P3811-[模板]乘法逆元【线性求逆元】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 祝福妈妈的话简短唯美 有关祝福妈妈的句子
- 下一篇: nssl1141,jzoj3470-最短