筛素数相关

筛素数

朴素筛法

原理表述

将每个数（无论素数还是和数）的倍数筛一遍

代码实现

#include<iostream>
using namespace std;
bool st[10010];
int primes[10010]，cnt = 0;
void get_primes(int n)
{
	st[1] = true;
	for(int i = 2 ; i <= n ; i++)
	{
		if(!st[i])   
			primes[++cnt] = i;
		for(int j = i+i ; j <= n ; j += i)	
			st[j] = true;
	}
}

int main()
{
	int n;
	cin >> n;
	get_primes(n);
	for(int i = 1 ; i <= cnt ; i++)
		cout << primes[i] << " "; 
	return 0;
 }

时间复杂度

O(nlogn)

埃氏筛法

原理表述

由于一个和数一定能分解为素数的乘积
那么将所有素数的倍数剔除，剩下的就是素数。

代码实现

#include<iostream>
using namespace std;
bool st[10010];
int primes[10010]，cnt = 0;
void get_primes(int n)
{
	st[1] = true;
	for(int i = 2 ; i <= n ; i++)
	{
		if(!st[i])
		{
			primes[++cnt] = i;	 
			for(int j = i+i ; j <= n ; j += i)	
				st[j] = true;
		}
	}
}

int main()
{
	int n;
	cin >> n;
	get_primes(n);
	for(int i = 1 ; i <= cnt ; i++)
		cout << primes[i] << " "; 
	return 0;
 } 

时间复杂度

O(nloglogn)

线性筛法 (欧拉筛法)

原理表述

在埃氏筛法的基础上,让每个合数只被它的最小质因子筛选一次,以达到不重复的目的

代码实现

#include<bits/stdc++.h>
using namespace std;
int primes[100000];
int st[100000],cnt=0;
void get_primes(int n)
{
	st[1]=true;
	for(int i=2;i<=n;i++)
	{
		if(!st[i])
		   primes[++cnt]=i;
		for(int j=1;j<=cnt&&i*primes[j]<=n;j++)
		{
			st[i*primes[j]] = 1;
			if(i%primes[j] == 0)  
				break;
		}
	}
}
int main()
{
	int n;
	cin>>n;
	get_primes(100000);
	for(int i=0;i<n;i++)
	{
		int x;
		cin>>x;
		if(!st[x])
		cout<<x<<" ";
	}
	return 0;
}

时间复杂度

O(n)