大素数阶乘取模模板

hdu-6608 Fansblog

#include
using namespace std;
#define ll long long

ll p;
ll mul_mod(ll a,ll b,ll mod){    //快乘法 基于快速幂的二分思想 
    ll ans=0;                    //由于考虑到取模数很大 快速幂会溢出 
    while(b){                    //必须使用该方法 
        if(b&1)                  //非递归版 
            ans=(ans+a)%mod;
        a=a*2%mod;
        b>>=1;
    }
    return ans;
}

ll pow_mod(ll a,ll n,ll mod){            //快速幂 递归版 
    if(n>1){                            
        ll tmp=pow_mod(a,n>>1,mod)%mod;
        tmp=mul_mod(tmp,tmp,mod);
        if(n&1) tmp=mul_mod(tmp,a,mod);
        return tmp;
    }
    return a;
}

bool Miller_Rabbin(ll n,ll a){//米勒拉宾素数判断函数主体
    ll d=n-1,s=0,i;    
    while(!(d&1)){            // 先把(2^s)*d 算出来 
        d>>=1;
        s++;
    }
    ll t=pow_mod(a,d,n);    //a^d取一次余判断 
    if(t==1 || t==-1)        //一或负一则可以声明这可能是质数 
        return 1;
    for(i=0;itab[i] && !Miller_Rabbin(n,tab[i]))
            return 0;
    }
    return 1;
}

void ex_gcd(ll a,ll b,ll &d,ll &x,ll &y){
	if(!b){
		d=a;
		x=1;
		y=0;
	}
	else {
		ex_gcd(b,a%b,d,y,x);
		y-=x*(a/b);
	}
}

ll inv(ll a,ll p){
	ll d,x,y;
	ex_gcd(a,p,d,x,y);
	return d==1?(x+p)%p:-1;
}

int main(){
	int pa;
	cin>>pa;
	while(pa--){
		scanf("%lld",&p);
		
		for(ll q=p-2;;q-=2){
			if(is_prime(q)){
				ll ans=1ll;
				for(ll i=q+1;i> 25ll) % mod * (1ll << 25) % mod;
    ll R = a * (b & ((1ll << 25) - 1)) % mod;
    return (L + R) % mod;
}
inline ll qpow_mul(ll x,ll n,ll mod) {
    ll ans = 1,a = x%mod;
    while(n) {
        if(n&1) {
            ans=mul(ans,a,mod);
        }
        a = mul(a,a,mod);
        n>>=1;
    }
    return ans;
}
  
  
inline ll qmul(ll a,ll b,ll mod) {
    ll ans = 0,x = a;
    while(b) {
        if(b&1) {
            ans = (ans + x%mod)%mod;
        }
        x = x*2%mod;
        b>>=1;
    }
    return ans%mod;
}
  
inline ll qpow_qmul(ll x,ll n,ll mod) {
    ll ans = 1,a = x%mod;
    while(n) {
        if(n&1) {
            ans=qmul(ans,a,mod);
        }
        a=qmul(a,a,mod);
        n>>=1;
    }
    return ans;
}
  
inline ll qpow(ll x,ll n,ll mod) {
    ll ans = 1,a = x;
    while(n) {
        if(n&1) {
            ans = ans*a%mod;
        }
        a = a*a%mod;
        n>>=1;
    }
    return ans;
}
  
inline ll gcd(ll m,ll n) {
    return n==0?m:gcd(n,m%n);
}