miller——rabin

突然发现自己在线性筛素数中有这个，忘了好久；

#include<iostream> 
#include<cstdio>
using namespace std;
long long num[10]={0,2,3,5,7,11,13,17,19};
long long mul(long long a,long long b,long long p)
{
    long long ans=1;
    a=a%p;
    while(b)
    {
        if(b&1)
            ans=(ans*a)%p;
        b>>=1;
        a=(a*a)%p;
    }
    return ans;
}
bool test(long long a)
{
    if(a==2)
        return true;
    if(a%2==0||a==1)
        return false;
    for(int i=1;i<=8;i++) 
        if(a==num[i])
            return true;
    long long t=0,temp=a-1,now;
    while((temp&1)==0)
    {
        temp>>=1;
        t+=1;
    }
    for(int i=1;i<=8;i++) 
    {
        now=mul(num[i],temp,a);
        long long nxt=now;
        for(int i=1;i<=t;i++)
        {
            nxt=(now*now)%a;
            if(nxt==1&&now!=1&&now!=a-1)
                return false;
            now=nxt;
        }
        if(now!=1)
            return false;
    }
    return true;
}
int main()
{
    long long n,m;
    scanf("%lld%lld",&n,&m);
    long long pass;
    for(int i=1;i<=m;i++)
    {
        scanf("%lld",&pass);
        if(test(pass))
            printf("Yes\n");
        else
            printf("No\n");
    }
}