埃式筛法比较好理解,线性筛法涉及到的东西蛮多。
标题名字,就是链接。
试除法
bool isPrime(int n) {
if (n==1)return false;
for (int i = 2; i <= n / i; i++) {//不使用iXi<=n 避免超出int范围
if (n % i == 0) return false;
}
return true;
}
埃氏筛法
#include<iostream>
using namespace std;
// 通过筛法求素数 即排除掉倍数的部分
const int N = 200010;
int primes[N], cnt;
bool st[N];
void get_prime(int n) {
for (int i = 2; i <= n; i++) {
if (!st[i]) {//如果没有被筛掉,则是质数
primes[cnt++] = i;
for (int j = i * 2; j <= n; j+=i) {
//j每次都是i的倍数
st[j] = true;//将倍数的值删除掉,就是正确的了
}
}
}
}
int main() {
get_prime(200000);
int k;
while (cin >> k)cout << primes[k - 1] << endl;
//从0开始计数,所以是k-1;
return 0;
}
线性筛法
#include<iostream>
using namespace std;
// 通过筛法求素数 即排除掉倍数的部分
const int N = 200010;
int primes[N], cnt;
bool st[N];
void get_prime(int n) {
for (int i = 2; i <= n; i++) {
if (!st[i]) primes[cnt++] = i;//放入质数
for (int j = 0; primes[j] * i <= n; j++) {
st[primes[j] * i] = true;//质数的倍数,一定是合数,所以被筛掉
if (i % primes[j] == 0)break;//i中的最小质因子是primes[j]
}
}
}
int main() {
get_prime(200000);
int k;
while (cin >> k)cout << primes[k - 1] << endl;
//从0开始计数,所以是k-1;
return 0;
}