逆元 递推求逆元

其实有些题需要用到1-p模p的所有逆元，这里p为奇质数。那么如果用快速幂求时间复杂度为O(p log(p))，
如果对于一个1000000级别的素数，这样做的时间复杂度是很高了。实际上有的算法，有一个递推式如下

inv[i]=(M-M/i)*inv[M%i]%M (其中M为模数，要求为奇质数)

它的推导过程如下：

设t=M/i,k=M%i，那么
    t*i+k≡0（Mod M）
    -t*i≡k(Mod M)
对上式两边同时除 i×k，进一步得到
    -t*inv[k]≡inv[i](Mod M)
再把和替换掉，最终得到
    inv[i]=(M-M/i)*inv[M%i]%M

初始化inv[1]=1，这样就可以通过递推法求出模素数的所有逆元了。

另外1-p模的所有逆元值对应1-p中所有的数，比如p=7，那么对应的逆元是1 4 5 2 3 6。

代码：

    #include<iostream>  
    #include<cstdio>  
    #include<cstdlib>  
    #include<cmath>  
    #include<algorithm>  
    #include<cstring>  
    using namespace std;  
    int A[100001];  
    int p;  
    int main()  
    {  
        cin>>p;  
        A[1]=1;  
        for(int i=2;i<=10;i++)  
            {  
                A[i]=(p-(p/i))*A[p%i]%p;  
                printf("%d %d %d\n",i,A[i],(i*A[i])%p);  
            }  
    }  

---

顺便复习一下逆元：

方法1：扩展欧几里得。  ax=1(mod P), gcd(a,p)=1， 其中x为a的逆元，就是我们所求，ax=PY+1,   ax-Py=1,  所以用扩展欧几里得可以求出x。

方法2：费马小定理：  如果模P是素数的话，那么inv(a)=pow(a,p-2)%p; 等式右边用快速幂运算可以得出。

^_^