具体方法的介绍可以看下面详细代码,有注释
//选出素数方法
#include<iostream>
#include<cstring>
#include<map>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn = 1010;
//方法1:试除法
bool s_prime(int a)
{
if (a <= 1) return false;
for (int i = 2; i*i < a; ++i) {
if (a%i == 0) return false;
}
return true;
}
//方法2:埃氏筛法,时间复杂度O(nloglog2(n)),缺点,比如10000会被重复筛除很多次(2,5,10...)
bool num[maxn];
int E_sieve(int n) //筛除0~n之间的非素数
{
for (int i = 0; i <= n; ++i) num[i] = false;
for (int i = 2; i*i <= n; ++i) {
if (!num[i]) {
for (int j = 2*i; j <= n; j = j + i) {
num[j] = true;
}
//改进,//当i为素数时,能够被i整除的数肯定不是素数,这里j=2*i可以改进为j=i*i,因为当i=5时,2*5,3*5,4*5在前面i为2,3,4时已经算过了
}
int k = 0;
int prime[maxn];
for (int i = 2; i <= n; ++i) {
if (!num[i]) prime[k++] = i; //筛出素数
}
}
}
//方法3:欧拉筛:时间复杂度O(n)
//对于每一个数,要么他是素数,要么它是素数的倍数
//所以用最小的质因子来筛,又因为主循环从2~n,所以每个数都会被筛到,而且只会被筛一次
int prime[maxn], i, j, countnum = 1;
bool num[maxn];
int oula(int n)
{
memset(num, true, sizeof(num));
for (i = 2; i <= n; ++i) { //i作为本身被筛掉或者作为因子来筛选其他数
if (num[i]) prime[countnum++] = i;
for (int j = 1; j < countnum; ++j) {
if (i*prime[j] > n) break; //超出所需求解的上限,退出
num[i*prime[j]] = false; //排除非素数
if (i%prime[j] == 0) break; //这里的i就是表示倍数,6*3=2*9,因为6是2的倍数,如果6*2的时候没有退出,那么当i到9的时候,18又会被筛掉一次
}
}
}
int main()
{
return 0;
}