素数判定 2（codevs 1702）

题目描述 Description

一个数，他是素数么？

设他为P满足(P<=263-1)

输入描述 Input Description

P

输出描述 Output Description

Yes|No

样例输入 Sample Input

2

样例输出 Sample Output

Yes

数据范围及提示 Data Size & Hint

算法导论——数论那一节
注意Carmichael Number

#include<iostream>
#include<cstdio>
#define LL long long
using namespace std;
LL n;
int s[]={,,,,,,,,,};
LL slowmi(LL a,LL b)
{
    LL tot=;
    for(LL i=a;i;i>>=)
    {
        if(i&) tot=(tot+b)%n;
        b=(b+b)%n;
    }
    return tot%n;
}
LL mi(LL a,LL b)
{
    LL tot=;
    for(LL i=a;i;i>>=)
    {
        if(i&) tot=slowmi(tot,b)%n;
        b=slowmi(b,b)%n;
    }
    return tot%n;
}
bool check()
{
    if(n==) return true;
    else if(n<||!(n&))  return false;
    for(int i=;i<;i++)//费马小定理10次计算
    {
        if(n==s[i])return true;
        if(mi(n-,s[i])!=) return false;
    }
    return true;
}
int main()
{
    cin>>n;
    if(check())printf("Yes");
    else printf("No");
    return ;
}