HDU 6608

Posted johnran

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了HDU 6608相关的知识,希望对你有一定的参考价值。

HDU 6608

题意:给以你一个质数,求小于它的最大质数的阶乘。

分析:Miller-Rabin快速判断素性,找到这个最大素数。此外,加上威尔逊定理推式子就好了。威尔逊定理描述的内容是对于一个正素数p:

? ((k-1)!modk)(= k-1)

有了这个定理,我们就可以很容易得到小于k最大素数阶乘结果了:

? (q!=frac{1}{(k-2)(k-3)...(q+1)} mod k)

#include<cstdio>

using namespace std;

typedef long long LL;

LL mul(LL x, LL y, LL mod){
    LL res=0;
    while(y){
        if(y&1)
            res=(res+x)%mod;
        x=(x+x)%mod;
        y>>=1;
    }
    return res;
}

LL powmod(LL x, LL y, LL mod){
    LL res=1;
    while(y){
        if(y&1)
            res=mul(res,x,mod);
        x=mul(x,x,mod);
        y>>=1;
    }
    return res;
}

int ft[12]={2,3,5,7,11,13,17,19,23,29,31,37};

bool check(LL a, LL n){
    LL x=n-1;

    int t=0;
    while((x&1)==0){
        x>>=1;
        ++t;
    }

    x=powmod(a,x,n);
    LL y;
    for(int i=1;i<=t;++i){
        y=mul(x,x,n);
        if(y==1&&x!=1&&x!=n-1)return true;
        x=y;
    }

    if(y!=1)return true;
    return false;
}

bool Miller_Rabin(LL x){
    if(x==2)return true;
    if(!(x&1)||x==1)return false;

    for(int i=0;i<12;++i){
        if(x<=ft[i])break;
        if(check(ft[i],x)) return false;
    }

    return true;
}

int main(){
    int T;scanf("%d", &T);
    while(T--){
        LL n;scanf("%lld",&n);
        LL q=n-1;
        while(!Miller_Rabin(q))--q;
        LL ans=1;
        for(LL i=q+1;i<n-1;++i)ans=mul(ans,powmod(i,n-2,n),n);
        printf("%lld
",ans);
    }

    return 0;
}

以上是关于HDU 6608的主要内容,如果未能解决你的问题,请参考以下文章

[hdu-6608] Fansblog 威尔逊定理 质数的密度分布 2019 多校 3

hdu多校第三场 1006 (hdu6608) Fansblog Miller-Rabin素性检测

h3c sr6608路由器IP限速

HDU4057 Rescue the Rabbit(AC自动机+状压DP)

HDU3247 Resource Archiver(AC自动机+BFS+DP)

2019 Multi-University Training Contest 3 - 1006 - Fansblog - 打表 - 暴力