P3811-[模板]乘法逆元【线性求逆元】
正題
評(píng)測(cè)記錄:https://www.luogu.org/recordnew/lists?uid=52918&pid=P3811
題目大意
求1~n1\sim n1~n中modpmod\ \ pmod??p的逆元。
解題思路
首先對(duì)于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)
左右兩邊同時(shí)乘上一個(gè)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)
反正本來(lái)就要%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
因?yàn)?span id="ozvdkddzhkzd" class="katex--inline">iii比(p%i)(p\%i)(p%i)小,按照遞推的順序我們?cè)谇蟪?span id="ozvdkddzhkzd" class="katex--inline">i?1i^{-1}i?1之前就已經(jīng)求出(p%i)?1(p\%i)^{-1}(p%i)?1了。但是我們還要保證不是負(fù)數(shù),所以我們可以直接計(jì)算
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]); }總結(jié)
以上是生活随笔為你收集整理的P3811-[模板]乘法逆元【线性求逆元】的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 祝福妈妈的话简短唯美 有关祝福妈妈的句子
- 下一篇: nssl1141,jzoj3470-最短