埃式筛法,线性筛求素数

1.最简单的算法O(√n)

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

}

2.埃式筛法算法O(nloglogn)

int prime[MAXN];
bool is_pri[MAXN];
int cnt=0;

void  isPrime(int n){
	for(int i=0;i<=n;i++) is_pri[i]=1;
	is_pri[0] = is_pri[1] = false;
	for(int i = 2;i <= n; i++){
		if(is_pri[i]){
			prime[++cnt] = i;
			for(int j=2*i;j<=n;j+=i) 当i满足质数时,筛选出i所有的倍数
				is_pri[j]=false; 
		}
	}
	
}

3.欧拉筛(线性筛)O(n)

时间复杂度最小的最优的筛法,相当于对埃式筛的进一步优化,因为埃式筛的一个合数会被多次筛选,比如20既会被2筛去,也会被5筛一次,所以这是冗余的操作

欧拉筛的理论基础是“任何一个合数都能写成两个质数相乘的形式”,由此可见,对于每一个合数,我们只需要用它的最小质因数把它筛去即可

int prime[MAXN];
bool is_pri[MAXN];
int cnt=0;

void  isPrime(int n){
	for(int i=0;i<=n;i++) is_pri[i]=1;
	is_pri[0] = is_pri[1] = false;
	for(int i = 2;i <= n; i++){
		if(is_pri[i]){
			prime[++cnt] = i;
		}
		for(int j=1;j<=cnt;j++){
			if(i*prime[j]>n)
                break;
			is_pri[i*prime[j]]=false;
			if(i%prime[j]==0)//这时此质数已经是这个合数的最小质因数
				break;
			
		}
		
	}
	
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
埃氏筛法(Eratosthenes筛法)和线性筛(Linear Sieve)都是用于求解素数的算。它们的主要区别在于筛选过程中的操作不同。 1. 埃氏筛法埃氏筛法是一种古老的求解素数的算,由古希腊数学家埃拉托斯特尼(Eratosthenes)提出。它的基本思想是从2开始,将所有2的倍数标记为非素数,然后找到下一个未被标记的数(即3),将所有3的倍数标记为非素数,依此类推。最后留下的未被标记的数就是素数。 C++实现埃氏筛法的代码如下: ```cpp #include <iostream> #include <vector> using namespace std; const int N = 1000000;vector<int> is_prime(N, true); void sieve_of_eratosthenes() { is_prime = is_prime = false; for (int i = 2; i * i < N; ++i) { if (is_prime[i]) { for (int j = i * i; j < N; j += i) { is_prime[j] = false; } } } } int main() { sieve_of_eratosthenes(); for (int i = 2; i < N; ++i) { if (is_prime[i]) { cout << i << " "; } } return 0; } ``` 2. 线性筛线性筛是一种改进的埃氏筛法,它将筛选过程从平方根优化到线性时间。线性筛的基本思想是对于每个素数p,筛选出所有小于等于p^2的合数。这样可以减少筛选的次数,提高效率。 C++实现线性筛的代码如下: ```cpp #include <iostream> #include <vector> using namespace std; const int N = 1000000; vector<int> is_prime(N, true); vector<int> primes; void linear_sieve() { is_prime = is_prime = false; for (int p = 2; p * p < N; ++p) { if (is_prime[p]) { for (int i = p * p; i < N; i += p) { is_prime[i] = false; } } } for (int p = 2; p < N; ++p) { if (is_prime[p]) { primes.push_back(p); } } } int main() { linear_sieve(); for (int i = 0; i < primes.size(); ++i) { cout << primes[i] << " "; } return 0; } ``` 这两种算都可以有效地求解素数,但线性筛相对于埃氏筛法在筛选次数上有很大优势。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值