大数乘法的逆元

 暴力法：  

 int i;
    for (i=1;;i++) { if (((long long int)(n)*i-an)%M==0) break; }

欧拉函数

long long inv( long long n )
{
    return pow( n, M - 2 )%M;
}

扩展欧几里得

//扩展欧几里德
void exp_gcd( LL a ,LL b,LL &x,LL &y) {
     if( b == 0 ) {
         x = 1;
         y = 0;
     }
     else {
          exp_gcd( b,a%b,x,y );
          LL t;
          t = x;
          x = y;
          y = t - a/b*y;
     }
}
//逆元
inline LL getNN(LL x) {
        LL now , t;
        exp_gcd( x, M,now,t );
        return (now%M+M)%M;
}

还有一个

int64 inv(int64 x) {  
    if(x == 1) return 1;  
    return  inv(M%x) * (M - M/x) % M;  
}  

void Inv( )
{    
    inv[1]=1;
    for(int i=2;i<=n;i++){
        inv[i]=M-(M/i)*inv[M%i]%M;
}