HDU多校第三场6608 Fansblog（米勒罗宾+威尔逊定理）

题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=6608

$题目说给你一个1e9-1e14的素数p，让你找到比p小的最大素数q，求q!\mathrm{m}\mathrm{o}\mathrm{d}p的值$

$我们知道相邻两个素数之间的距离不会很大，一般小于logn，所以我们可以从p-1开始枚举q，用米勒罗宾判定是否为素数即可找到q，根据威尔逊定理(p-1)！\equiv -1(modp),用-1乘上p-1,p-2,p-3\dots \dots q+1的逆元即可$

#include <bits/stdc++.h>
typedef long long ll;
using namespace std;
const int maxn = 2e5 + 500;
int t;
ll p,q,ans;
ll qmul(ll a, ll b, ll mod)
{
    ll res = 0;
    while(b)
    {
        if (b & 1)
            res = (res + a) % mod;
        a = (a << 1) % mod;
        b >>= 1;
    }
    return res;
}
ll qpow(ll x, ll n, ll mod)
{
    ll res = 1;
    while(n)
    {
        if (n & 1)
            res = qmul(res, x, mod);
        x = qmul(x, x, mod);
        n >>= 1;
    }
    return res;
}
bool Miller_Rabin(ll n){
    if(n==2)return true;
    if(n<2||!(n&1))return false;
    int cishu=2,ji=0;
    ll m=n-1;
    while(!(m&1)){
        ji++;
        m>>=1;
    }
    srand(100);
    while(cishu--){
        ll a=rand()%(n-1)+1;
        ll x=qpow(a,m,n),tmp=0;
        for(int i=0;i<ji;i++){
            tmp=qmul(x,x,n);
            if(tmp==1&&x!=1&&x!=n-1)return false;
            x=tmp;
        }
        if(tmp!=1)return false;
    }
    return true;
}
int main()
{
    scanf("%d", &t);
    while(t--)
    {
        scanf("%lld", &p);
        q = p - 1;
        while(!Miller_Rabin(q))
            --q;
        printf("%lld\n", q);
        ans = p - 1;
        for (ll i = p - 1; i > q; --i)
            ans = qmul(ans, qpow(i, p - 2, p), p);
        printf("%lld\n", ans);
    }
    return 0;
}