hdu5391Zball in Tina Town（素数鉴别的大猜想）

题目：点击打开链接

思路：balabala一大堆，其实本意就是考你威尔逊定理，关于威尔逊定理

如果p为素数，则(p-1)! mod p = p-1;

如果p为合数，则(p-1)! mod p = 0;当p等于4时除外，当p=4时，答案为2而不是0；

所以这个题就符合上边所说的定理，唯一的cheak就是p=4时。

解题：那么既然知道了是判断素数的问题，那么就好做了，但是判断素数的方法有这么多，选择哪一个呢，不妨让我们来试试。

1.素数表筛选法：

前边：今天算是见识了，疯狂输出，狂打表，结果疯狂RE，真心不知道为啥RE，最后终于知道为啥了，但是我真的不知道为什么会RE啊，好气啊，今天一直在纠结素数打表的问题，感觉爆炸了。

素数表的模板就在这一起整理了把

最后的结果是180ms左右

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
const int maxn = 1e6+10;//尽量开1e6，re不要怪我
bool isprime[maxn];
int prime[maxn];
//打印素数表的模板
void build()
{
    int i,j;
    memset(isprime,true,sizeof(isprime));
    isprime[1]=false;
    int nprime=0;
    for(i=2;i<maxn;i++){  //wtm为啥写sqrt（maxn）就RE
        if(isprime[i]){
            prime[nprime++]=i;
            for(j=i+i;j<maxn;j+=i) isprime[j]=false;//wtm为啥写i*i就RE！！！什么鬼
        }
    }
}

bool judgeprime(int x)
{
    for(int i=0;prime[i]*prime[i]<=x;i++)
        if(x%prime[i]==0)
        return false;

    return true;
}

int main()
{
    int T;
    build();
    scanf("%d",&T);
    while(T--){
        int x;scanf("%d",&x);
        if(judgeprime(x)) printf("%d\n",x-1);
        else{
            if(x==4) printf("2\n");
            else  printf("0\n");
        }
    }
    return 0;
}

2.miller_Rabin算法

前边：这个也是很奇怪，感觉算法很快，但是一直T，最后把乘法改进了下，结果速度变得奇快！！！这个乘法的快速计算的方法，没怎么明白。先收藏下来，不过miller_rabin倒是挺快的，但是最后的结果，没比打表快到那里去。速度是280ms

代码：

include<iostream>
#include<cstdio>
#include<ctime>
#include<cstdlib>
using namespace std;
typedef long long  ll;
const int Times=3;//MillerRabin检测次数

ll quickmul(ll x,ll y,ll Z)//大数相乘取模，非常神奇的一个快速乘法的计算方法
{
    ll tmp=x/(long double)Z*y+1e-3;
    return (x*y+Z-tmp*Z)%Z;
}

ll quickmod(ll a,ll b,ll c)
{
    ll ans=1;
    while(b){
        if(b&1) ans=(ans*a)%c;
        a=quickmul(a,a,c);
        b>>=1;
    }
    return ans;
}

bool millerRabin(ll x)//只利用费马小定理的测试函数 a^(x-1)%x==1则 x为素数
{
    if(x<=1) return false;
    if(x==2) return true;
    if(!(x&1)) return false;

    for(int i=0;i<Times;i++){
        ll a=rand()%(x-1)+1;//随机生成一个数 来进行费马小定理检测
        if(quickmod(a,x-1,x)!=1) return false;
    }
    return true;
}

int main()
{
    int T;
    while(scanf("%d",&T)!=EOF){
        while(T--){
            ll x;scanf("%lld",&x);
            if(millerRabin(x)) printf("%lld\n",x-1);
            else{
                if(x==4) printf("2\n");
                else printf("0\n");
            }
        }
    }
}