素数筛 codevs 1675 大质数 2

时间:2023-03-09 18:49:47
素数筛  codevs 1675 大质数 2

1675 大质数 2

 时间限制: 1 s
 空间限制: 1000 KB
 题目等级 : 钻石 Diamond
题解
 查看运行结果
题目描述 Description

小明因为没做作业而被数学老师罚站,之后数学老师要他回家把第n个质数找出来。

小明于是交给聪明的你。请你帮忙!【wikioi-1530】

…………………………以上为背景…………………………

老师怀疑小明仅仅是找到第n个质数,于是又叫小明把1到n以内(不包括n)的质数全部找出来。小明又找到了你……

输入描述 Input Description

一个正整数n。

(1<=n<=1000000)

输出描述 Output Description

n以内的质数,每个一行。

样例输入 Sample Input

233

样例输出 Sample Output

2
3
5
7
11
13
17
19
23
29
31
37
41
43
47
53
59
61
67
71
73
79
83
89
97
101
103
107
109
113
127
131
137
139
149
151
157
163
167
173
179
181
191
193
197
199
211
223
227
229          //(不含n=233)

Euler筛法(线性筛法):
 #include<iostream>
using namespace std;
#include<cstdio>
#define N 1001000
#include<bitset>
bitset<N>check;
int prim[N],tot=;
int main()
{
int n;
cin>>n;
check.set();
for(int i=;i<n;++i)
{
if(check[i])
{
printf("%d\n",i);
prim[tot++]=i;
}
for(int j=;j<tot;++j)
{
if(prim[j]*i>=n) break;
check[prim[j]*i]=false;
if(i%prim[j]==) break;
}
} return ;
}

埃氏筛法:

 /*/*这个程序在i==四万多的时候溢出了,而欧拉筛法因为有“if(i%prim[j]==0) break;
”,所以一般不会乘法溢出*/*/
#include<iostream>
using namespace std;
#include<cstdio>
#define N 1001000
#include<bitset>
bitset<N>check;
int prim[N],tot=;
int main()
{
int n;
cin>>n;
check.set();
for(int i=;i<n;++i)
{
if(check[i])
{
printf("%d\n",i);
for(int j=i*i;j<n;j+=i)
check[j]=false;
} } return ;
}