求逆元的2种模板

(1) exgcd 求逆元

//x->a关于b的逆元
//y->b关于a的逆元
typedef long long LL
void ex_gcd(LL a, LL b, LL &x, LL &y, LL &d){
      if (!b) {d = a, x = 1, y = 0;}
      else{
          ex_gcd(b, a % b, y, x, d);
          y -= x * (a / b);
      }
  }
LL inv(LL t, LL p){//如果不存在,返回-1 
     LL d, x, y;
     ex_gcd(t, p, x, y, d);
     return d == 1 ? (x % p + p) % p : -1;
 }

(2)快速幂+费马小定理求逆元

// 快速幂求逆元
LL pow_mod(LL a, LL b, LL p){//a的b次方求余p 
     LL ret = 1;
     while(b){
         if(b & 1) ret = (ret * a) % p;
         a = (a * a) % p;
         b >>= 1;
     }
     return ret;
 }
LL Fermat(LL a, LL p){//费马求a关于b的逆元 
         return pow_mod(a, p-2, p);
 }

猜你喜欢

转载自blog.csdn.net/bestsort/article/details/81132982