大佬部落格
個人比較菜會用板子就好了。
送上例題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;
}