素数判断

普通方法

复杂度：O(n*sqrt(n))

bool is_prime(int n){
    int m=(int)sqrt(n);
    for(int i=2;i<=m;++i){
        if(n%i==0)
            return flase;
    }
    return true;
}

6倍法

原理：当x大于等于5时可以发现x可以写成6n-1,6n,6n+1,6n+2,6n+3,6n+4的形式，其中除了6n+1 6n-1这两个数可能为素数外，其他的都不可能为素数。故当判断时可以分为以下三部：

1. 判断5以下的情况

2. 当x>=5时

   a)判断是否与6n(n为正整数)相邻，如果相邻，继续下一步，否则返回；

   b)依次用6n相邻的数去除这个数

bool is_prime(int n){
   if(n==2||n==3){
       return true;
   }
   if(n%6!=1&&n%6!=5)
   return false;
   int m=(int)sqrt(n);
   for(int i=5;i<=m;i+=6){
       if(n%i==0||n%(i+2)==0)
       return false;
   }
   return true;
}

埃拉托斯特尼筛

vector<int> prime;
bool judge[range];
memset(judge,true,sizeof(judge)); 
judge[0]=judge[1]=flase;          //不是素数
for(int i=2;i<=sqrt(n);++i){      //
    if(judge[i]){
        prime.push_back(i);
        for(long long j=i*i;j<=maxn;j+=i){
            judge[j]=false;
        }
    }
}

优化：欧拉筛

在上面算法上加以改进，保证所有数只判断一次，不然可能被判断两次，如6=2*3；会被2和3都赋值

memset(judge,true,sizeof(judge))
judge[0]=judge[1]=flase;
for(int i=2;judge<=(int)sqrt(n);++i){
    if(judge[i]){
        prime.push_back(i);
        cnt++;
    }
    for(int j=0;j<cnt&&i*prime[j]<=range;++j){
        prime[i*pirme[j]]=flase;
        if(i%prime[j]==0){
            break;
        }
    }
}

仅作为个人纪录，如有错误，敬请指正。🙂