素数筛--

目的

快速的得到素数

一般的判断方法

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

▶ 时间复杂度为O(√𝑛）

埃式筛

①2是最小的素数，将表中所有是2的倍数删除

②3成为最小的素数，将表中所有是3的倍数删除

③以此类推

相关代码

#define maxn 10000000
bool vis[maxn];  //判断是否被标记过
int prime[maxn],x;  //记录素数

void isprime(int n) //埃氏筛
{
    for(int i=2;i<=n;i++)
    {
        if(!vis[i]) prime[x++]=i;  //把2，3，5等的素数记录下来
        for(int j=2;j*i<=n;j++)
        {
            vis[i*j]=true;  //标记，都不是素数了
        }
    }
}

▶ 时间复杂度为O(NloglogN)，非常接近线性

欧拉筛

由上面的过程可以看到，有的数字被删除了好几次

想要每个数只被其最大非自身因数筛选一次

8 被 4 筛去
10，15，25 被5筛去
12 被6筛去

const int N = 1e8+5;
int p[N];
int v[N];//存储 i 的最小质因子
int cnt;
void primes(int n)
{
	 for(int i = 2; i <= n; ++i){
		 if(!v[i]){//i 是质数
			 v[i] = i; p[++cnt] = i;
 }
 
 //给当前的数 i 乘上一个质因子
 for(int j = 1; j <= cnt; ++j){ 
		 //i 有比 p[j]更小的质因子,或着超出 n 的范围,停止循环
   		if(p[j] > v[i] || p[j] > n/i) break;
 				//p[j]是合数 i*p[j]的最小质因子
				 v[i*p[j]] = p[j];
      }
   }
} 

相关题目

① 两个四位数素数变换的最小次数