Miller Rabbin （判断大素数）

费马小定理：

a为整数，n是素数，且a，n互质，则有a^(n-1)≡1(mod n) ，即：a^(n-1)模n得1。

快速判定一个数是否为素数的方法：

如果存在一个整数a，使得a^(n-1)≡1(mod n) ，则称n为基于a的伪素数，当有多个满足关系的a时，则n为素数的概率趋向于1。所以取多个a测试一下即可。

模板代码如下：

# include<iostream>
# include<cstdio>
# include<cstring>
# include<cstdlib>
# include<algorithm>
using namespace std;
# define ll long long
ll mypow(ll a,ll b,ll m)
{
    if(b==0)
        return 1;
    if(b==1)
        return a%m;
    ll temp=mypow(a,b/2,m);
    temp*=temp;
    temp%=m;
    if(b&1)
        temp*=a;
    temp%=m;
    return temp;
}
bool Miller_Rabbin(ll x)
{
    if(x==2)
        return true;        ///2要直接判断
    for(int i=1;i<=50;++i){
        ll a=rand()%(x-2)+2;
        if(mypow(a,x-1,x)!=1)
            return false;
    }
    return true;
}
int main()
{
    ll n;
    while(scanf("%lld",&n)!=EOF)
    {
        if(Miller_Rabbin(n))
            cout<<"Yes"<<endl;
        else
            cout<<"No"<<endl;
    }
    return 0;
}