Miller-Rabin随机性素数测试算法

大佬博客
个人比较菜会用板子就好了。
送上例题hdu 2138
虽然暴力可以过但是还是用来学算法吧。

#include<cstdio>
#include<cstring>
#include<ctime>
#include<iostream>
#include<algorithm>
#define LL long long
using namespace std;
const int maxn=1e6+10;
int S;//规定随机多少次，次数越多正确率越高
LL quick_mod_multi(LL a,LL b,LL mod)//快速求a*b%mod
{
    LL t=0;
    a%=mod;
    while(b)
    {
        if(b&1) t=(t+a)%mod;
        a=(a<<1)%mod;
        b>>=1;
    }
    return t;
}
LL quick_mod_power(LL a,LL b,LL mod)//快速求a^b%mod
{
    LL res=1;
    a%=mod;
    while(b)
    {
        if(b&1)
        res=quick_mod_multi(res,a,mod);
        b>>=1;
        a=quick_mod_multi(a,a,mod);
    }
    return res;
}
bool miller_rabbin(LL n)//素性检验
{
    if(n==2) return true;
    if(n<2||!(n&1)) return false;
    int t=0;
    LL a,x,y,u=n-1;
    while((u&1)==0) t++,u>>=1;
    srand((unsigned)time(NULL));
    for(int i=0;i<S;i++)//随机多少次
    {
        a=rand()%(n-1)+1;
        x=quick_mod_power(a,u,n);
        for(int j=0;j<t;j++)
        {
            y=quick_mod_multi(x,x,n);
            if(y==1&&x!=1&&x!=n-1)
                return false;
            x=y;
        }
        if(x!=1) return false;
    }
    return true;
}
int main()
{
    S=3;
    int n;
    while(~scanf("%d",&n))
    {
        int sum=0;
        for(int i=0;i<n;i++)
        {
            LL a;
            scanf("%lld",&a);
            if(miller_rabbin(a)) sum++;
        }
        printf("%d\n",sum);
    }
}

这里有一个比暴力能快一点的算法很神奇刚发现素数还有这种规律。

#include<cstdio>
#include<cmath>
int pr[8]={4,2,4,2,4,6,2,6};//大于7的素数两个的差值是（这些数 or 相邻的加起来），反正我测得200以内是对的
int prime(int n)
{
    int i=7,j;
    if(n<2)
        return 0;
    if(n==2||n==3||n==5)
        return 1;
    if(!(n%2&&n%3&&n%5))
        return 0;
    for(;i<=sqrt(n);)
    {
        for(j=0;j<8;j++)
        {
            if(n%i==0)
             return 0;
            i+=pr[j];
        }
        if(n%i==0)
         return 0;
    }
    return 1;
}
int main()
{
    int i,n,m,s;
    while(scanf("%d",&n)!=EOF)
    {
        s=0;
        for(i=0;i<n;i++)
        {
            scanf("%d",&m);
            if(prime(m))
                s++;
        }
        printf("%d\n",s);
    }
    return 0;
}