C++算法——埃氏筛

5 篇文章 0 订阅
3 篇文章 0 订阅

C++判断素数:埃氏筛

思路

这个算法是利用打表的方法来计算的:
首先,我们要知道一个特性
就是一个质数的倍数,一定是一个合数

利用这个特性
我们可以写出以下代码

for (int i = 2; i * i <= n; i++)
{
	if (!prime[i] == true)
	{
		for (int j = i + i; j <= n; j++)
		{
			prime[j] = true;//这个数就是质数的倍数
		}
	}
}

这段代码就是埃氏筛的核心代码(就这些)
上面的判断是为了以防重复打表,下面的for是为了标记这个质数所有小于n的合数

代码

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1E9 + 7;

bool prime_n[MAXN];

int EhrlichSieve(bool prime[], int n)
{
	for (int i = 2; i * i <= n; i++)
	{
		if (!prime[i] == true)
		{
			for (int j = i + i; j <= n; j++)
			{
				prime[j] = true;
			}
		}
	}
	return 0;
}
int main()
{
	EhrlichSieve(prime_n, 1000000);
	int q;
	cin >> q;
	while (q--)
	{
		int x;
		cin >> x;
		if (prime_n[x] == false)
		{
			cout << x << "是素数\n";
		}
		else
		{
			cout << x << "不是素数\n";
		}
	}
	return 0;
}
//这里的"EhrlichSieve"使用的是百度机翻,不知道对不对
埃氏法(Sieve of Eratosthenes)是一种用于求解素数的经典算法。它的基本思想是从2开始,将每个素数的倍数标记为合数,直到遍历完所有小于等于n的数。最后剩下的未被标记的数即为素数。 以下是使用埃氏法求解1到n的素数的C++代码示例: ```cpp #include <iostream> #include <vector> std::vector<int> sieveOfEratosthenes(int n) { std::vector<bool> isPrime(n + 1, true); // 初始化标记数组,默认所有数都是素数 std::vector<int> primes; // 存储素数的容器 for (int p = 2; p * p <= n; p++) { if (isPrime[p]) { for (int i = p * p; i <= n; i += p) { isPrime[i] = false; // 将p的倍数标记为合数 } } } for (int p = 2; p <= n; p++) { if (isPrime[p]) { primes.push_back(p); // 将未被标记的数加入素数容器 } } return primes; } int main() { int n; std::cout << "请输入一个正整数n:"; std::cin >> n; std::vector<int> primes = sieveOfEratosthenes(n); std::cout << "1到" << n << "的素数有:"; for (int prime : primes) { std::cout << prime << " "; } std::cout << std::endl; return 0; } ``` 运行示例: ``` 请输入一个正整数n:20 1到20的素数有:2 3 5 7 11 13 17 19 ``` 在上述代码中,我们使用了一个布尔类型的标记数组`isPrime`来记录每个数是否为素数。初始时,将所有数都标记为素数。然后从2开始遍历,如果当前数是素数,则将其倍数标记为合数。最后遍历标记数组,将未被标记的数加入素数容器。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值