素数筛

素数判断：

bool is_prime(int a)
{
	for(int i=2;i*i<=a;i++)	
		if(a%i==0)		
			return false;
	return true;		
}

素数筛

两部分构成：埃氏筛+欧拉筛

• 埃氏筛：在给定范围内，找出n的倍数的所有的数并筛去，时间复杂度位O(nloglogn)。

void ai_prime(int m)
{
	for(int i=2;i<=m;i++)
	{
		if(!visit[i]){
			cnt++;
			prime[cnt]=i;
		}
		for(int j=2;j*i<=n;j++)
		{
			visit[j*i]=1;
		}	
	}
}

• 欧拉筛：让每个合数减少重复筛选的次数，时间复杂度O(n)。

   eg：2与4的公倍数会被重复筛选
  

void ou_prime(int m)
{
	for(int i=2;i<=m;i++)
	{
		if(!visit[i]){
			cnt++;
			prime[cnt]=i;
		}
		for(int j=1;j<=cnt;j++)//从已有的素数中筛选（判断） 
		{
			if(i*prime[j]>m)	break;//超出给定范围 
			visit[i*prime[j]]=1;//乘以相应的倍数 
			if(i%prime[j]==0)	break;//过滤点重复的筛选 
		}
	}	
}

有空看看：
数论部分——Matrix67
米勒罗宾素数检测法——smiling~

https://blog.csdn.net/dy416524/article/details/86431057