【数论】线性求[1,p-1]所有数逆元的方法

以前求逆元只会费马小定理和exgcd,看到别人都用递推求自己不会,今天学习了一下。
我们要在线性时间内求出11,21,(p1)1(modp)
p

111(modp)111(modp)

aa11(modp)    1<a<p

k=par=p  mod  a
p=ka+r       0<r<aka+r0(modp)(ka+r)a1r10(modp)kr1+a10(modp)a1kr1(modp)a1=pa(p  mod  a)1  mod  p=(ppa)(p  mod  a)1  mod  p

inv[1]=1;
for(int i=2;i<=n;i++)
    inv[i]=(p-p/a)*inv[p%a]%p;

同时,也可以据此来递归求出逆元,每次时间复杂度为O(log2n)

int Get_inv(int n){
    if(n==1)
        return 1;
    return (p-p/n)*(Get_inv(p%n))%p;
}

p
一组例子,n=7,p=18

posted @ 2016-05-26 17:53  outer_form  阅读(331)  评论(0编辑  收藏  举报