一般单个筛选素数
O(n^2)
int prime(int a)
{
for(int i = 1;i <= a/i;i++)
{
if(a%i == 0) return i;
}
return a;
}
埃式筛:用于一连串的素数筛选
O(nloglogn)
int prime(int a)
{
bool prime[100001] = {false};
for(int i = 2;i <= n;i++)
{
if(!prime[i])
{
printf("%d\n",i);
for(int j = i*i;j <= n;j += i) prime[j] = true; //有重复,2*5 5*2 多的情况就是5*1/2/3/4 都会多余
// for(int j = 2*i;j <= n;j += i) prime[j] = true;
}
}
}
同样 欧拉筛(线性筛):也用于一连串素数筛选
O(n)
#include<stdio.h> //O(n) 用于找前n个质数
#include<string.h>
typedef long long ll;
const int maxn = 1001; //找1001前的质数,存在prime数组中
int prime[maxn];
bool judge[maxn];
int main()
{
int num = 0;
memset(judge,1,sizeof(judge));
for(int i = 2;i <= maxn;i++)
{
if(judge[i]) prime[num++] = i; //存质数
for(int j = 0;j < num;j++)
{
if(i*prime[j] > maxn) break; //超出范围就结束循环
judge[i*prime[j]] = 0; //去合数
if(i%prime[j] == 0) break; //非最小质数筛选就结束
}
}
for(int i = 0;i < num;i++) printf("%d\n",prime[i]); //所有的质数都保存在prime数组中,judge只是为了去除合数使用,所以可以用bool型
return 0;
}
还有min_25 筛法,有时间的同学可以了解一下…
素数的产生有一定规律
1.小于 x 的素数个数满足 x/lnx 个
2.素数一定能表示 6n + 1 的型式,及所有合数都是 1 - 6 的倍数
下面有一个图片特别清楚(理解着看一遍)