简单数学模板

void get_prime(int n){//造素数
    check[1]=1;
    for(int i=2;i<=n;i++){
        if(!check[i])prime[tot++]=i;
        for(int j=0;j<tot;j++){
            if(i*prime[j]>n)break;
            check[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
}
void P_M(int n){//求1~n的逆元
    a[1]=1;
    for(int i=2;i<=n;i++)
        a[i]=1ll*(M-M/i)*a[M%i]%M;
}
ll gcd(ll a,ll b){//最大公因数 
    while(b)b^=a^=b^=a%=b;
    return a;
}
ll mul(ll a,ll b,ll p){//快速乘 
    ll res=0;
    while(b){
        if(b&1)res=(res+a)%p;
        a=(a+a)%p;
        b>>=1;
    }
    return res;
}
bool check_prime(ll x,ll p){
    ll i=ceil(log(p-1.0)/log(2.0)),d=1;
    for(;i>=0;i--){
        ll tmp=d;d=mul(d,d,p);
        if(d==1&&tmp!=1&&tmp!=p-1)return false;
        if((p-1)&(1<<i))d=mul(d,x,p);
    }
    return d==1;
}
bool miller_rabin(ll n){//判素数
    if(n==1)return false;
    for(int i=0;i<9;i++){
        if(n==prime[i])return true;
        if(n%prime[i]==0)return false;
    }
    for(int i=0;i<9;i++)
        if(prime[i]<n&&!check_prime(prime[i],n))return false;
    return true;
}
ll Pollard_Rho(ll n,ll c){
    ll x,y,d,i=1,k=2;
    x=mul(rand(),rand(),n-1)+1;
    y=x;
    while(1){
        i++;
        x=(mul(x,x,n)+c)%n;
        d=gcd(y-x,n);
        if(d>1&&d<n)return d;
        if(y==x)return n;
        if(i==k){
            y=x;
            k<<=1;
        }
    }
}
ll pri_fac[M],now[M];int tp,top;
void dec_prime(ll n){//分解质因数 
    //dfs,不想写了 
    //在此递归调用Pollard_Rho
    //将n分解 
}