乘法逆元整理

乘法逆元

1.费马小定理

如果 $ a,p $ 互质，那么 $ a^{p-1} ≡ 1(mod p) $
又由逆元方程知 $ a*x ≡ 1 (mod p) $ ，得到 $ ax ≡ a^{p-1}(mod p) $
所以逆元 $ x $ 为 $ a^{p - 2} mod p $，快速幂求解即可。

    #include<iostream>
    #include<cstdio>
    #include<cstring>
    #include<algorithm>
    #include<cmath>
    
    using namespace std;
    
    #define LL long long
    
    LL n,p;
    
    inline LL fast_pow(LL a,LL b,LL p) {
        LL ans = 1;
        a %= p;
        while(b) {
            if(b&1) ans = ans * a % p;
            a = a * a % p;
            b >>= 1;
        }
        ans %= p;
        return ans;
    }
    
    int main() {
        scanf("%lld%lld",&n,&p);
        for(int i = 1 ; i <= n ; i++) 
            printf("%lld \n",fast_pow(i,p-2,p));
        return 0;
    }

2.线性求逆元

这个不会证。
不过不会证问题也不大，可以直接当一个结论用。
$ inv[i] = p - \frac{p}{i} * inv[p % i] % p$
不过要注意 $ inv[1] = 1,inv[0] = tan90^o = 0 $
$ for $ 循环从2开始，时间复杂度 $ O(n) $ 。

  
  using namespace std;
  
  #define LL long long
  const int N = 3e6 + 10;
  
  LL n,p,inv[N];
  
  int main() {
      scanf("%lld%lld",&n,&p);
      inv[1] = 1;inv[0] = 0;
      for(int i = 2; i <= n ; i++) 
          inv[i] = p - (p / i) * inv[p % i] % p;
      for(int i = 1 ; i <= n ; i++)
          printf("%lld \n",inv[i]);
      return 0;
  }

转载于:https://www.cnblogs.com/Repulser/p/11207140.html