大素数阶乘取模模板

hdu-6608 Fansblog

#include<bits/stdc++.h>
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;i<s;i++){                //不是的话继续乘上s个2 
        if(t==n-1)            //(n-1)*(n-1)%n=1 这一步是优化 
            return 1;
        t=mul_mod(t,t,n);    // 快乘 
    }
    return 0;
}

bool is_prime(ll n){
    ll i,tab[4]={3,4,7,11};//本来应该取[1,n]内任意整数 
    for(i=0;i<4;i++){                //但一般这几个数足以,不需要太多组测试 
        if(n==tab[i])
            return 1;        //小判断小优化~ 
        if(!n%tab[i])
            return 0;
        if(n>tab[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<p-1;i++){//利用威尔逊定理把 q的阶乘模p 化简为 q+1到p-2的连乘模p 
					ans=mul_mod(ans,i,p);
				}
				ans=mul_mod(1,inv(ans,p),p);
				printf("%lld\n",ans);
				break;
			}
		}
	}
	return 0;
}

inline ll mul(ll a,ll b,ll mod){
//改良版快速乘 趋近O(1)
//无法解决洛谷屠龙勇士的毒瘤溢出还是得用龟速乘
//只适用于小数据中的溢出  且两个均为正整数
    ll L = a * (b >> 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);
}