c语言阶乘取模,大素数阶乘取模模板

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> 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);

}