筛质数
程序员文章站
2022-04-01 15:41:47
...
题目链接
题意:给定一个正整数n,请你求出1~n中质数的个数。
输入格式
共一行,包含整数n。
输出格式
共一行,包含一个整数,表示1~n中质数的个数。
数据范围
1≤n≤1e6
输入样例:
8
输出样例:
4
思路1 (朴素做法):把1~n中2的倍数,3的倍数, 4的倍数……都删除,那么剩下的数就都是质数了。
代码实现: O(n * logn)
#include<iostream>
#include<cstdio>
using namespace std;
const int N = 1e6 + 5;
int n, x;
int primes[N], res;
bool st[N];
void get_primes(int n)
{
for(int i = 2; i <= n; i ++){
if(!st[i]){
res ++ ;
}
for(int j = i + i; j <= n; j += i)
st[j] = 1;
}
}
int main()
{
cin >> n;
get_primes(n);
cout << res << endl;
return 0;
}
埃氏筛法:我们会发现4的倍数也是2的倍数,所有我们只需要把所有质数的倍数删除就好。
代码实现: O(n * loglogn)
#include<iostream>
#include<cstdio>
using namespace std;
const int N = 1e6 + 5;
int n, x;
int primes[N], res;
bool st[N];
void get_primes(int n)
{
for(int i = 2; i <= n; i ++){
if(!st[i]){
res ++ ;
for(int j = i + i; j <= n; j += i)
st[j] = 1;
}
}
}
int main()
{
cin >> n;
get_primes(n);
cout << res << endl;
return 0;
}
当然这还可以优化一下,变成一个基本上线性的算法。
代码实现: O(n)
#include<iostream>
#include<cstdio>
using namespace std;
const int N = 1e6 + 5;
int n, x;
int primes[N], res;
bool st[N];
void get_primes(int n)
{
for(int i = 2; i <= n; i ++){
if(!st[i]) primes[res ++ ] = i;
//这里不用加j <= cnt的限制条件,因为当i是合数的时候,primes[j]遇到i的最小因子就会停下来,当i是质数的时候,primes[j] = n / i后也会停下来
for(int j = 0; primes[j] <= n / i; j ++){
st[primes[j] * i] = 1;
//primes[j]一定是i的最小质因子
if(i % primes[j] == 0) break;
//一定会筛掉所有合数,例如4会把它的最小质因子2筛掉
}
}
}
int main()
{
cin >> n;
get_primes(n);
cout << res << endl;
return 0;
}
上一篇: 51单片机攻略—外部中断、定时器
推荐阅读
-
C语言程序设计100例之(12):Eratosthenes筛法求质数
-
素数(质数)判断、打印素数表(Eratosthenes筛法)、质因子分解————附完整代码
-
质数怎么判断(线性筛查质数)
-
【leetcode系列】【算法】【简单】计数质数(埃拉托斯特尼筛法sieve of Eratosthenes)
-
埃式筛/厄拉多塞筛法/Sieve_of_Eratosthenes/计数质数
-
【代码超详解】HDU 2138 How many prime numbers? 有多少个质数?(质数埃筛打表 + 判定,31 ms)
-
C语言程序设计100例之(12):Eratosthenes筛法求质数
-
线性筛质数
-
筛质数
-
质数怎么判断(线性筛查质数)