埃拉托色尼筛法
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll maxn=1e5;
ll n;
ll prime[maxn];
void eratos(){
for(ll i=0;i<n;i++){
prime[i]=true;
}
prime[0]=prime[1]=false;
for(ll i=2;i*i<=n;i++){
if(prime[i]){
ll j=i+i;
while(j<=n){
prime[j]=false;
j=j+i;
}
}
}
return;
}
int main(){
cin>>n;
eratos();
for(ll i=0;i<n;i++){
if(prime[i]==1){
cout<<i<<" ";
}
}
return 0;
}
欧拉筛法(比埃筛控制的区间更大)
这种方法比较好理解,初始时,假设全部都是素数,当找到一个素数时,显然这个素数乘上另外一个数之后都是合数
把这些合数都筛掉,即算法名字的由来。但仔细分析能发现,这种方法会造成重复筛除合数,影响效率。
比如10,在i=2的时候,k=215筛了一次;在i=5,k=56 的时候又筛了一次。所以,也就有了快速线性筛法。
利用了每个合数必有一个最小素因子。每个合数仅被它的最小素因子筛去正好一次。所以为线性时间。
代码中体现在:
if(i%prime[j]==0)break;
prime数组 中的素数是递增的,当 i 能整除 prime[j],那么 i*prime[j+1] 这个合数肯定被 prime[j] 乘以某个数筛掉。
因为i中含有prime[j], prime[j] 比 prime[j+1] 小。接下去的素数同理。所以不用筛下去了。
在满足i%prme[j]==0这个条件之前以及第一次满足改条件时,pr[j]必定是pr[j]*i的最小因子
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll maxn = 1e5;
ll prime[maxn];
ll vis[maxn];
ll index;
ll n;
void getprime()
{
for (int i = 2; i <= n; i++)
{
if (vis[i]==0)
prime[++index] = i;
for (int j = 1; j <= index && prime[j] * i <= n; j++)
{
vis[prime[j] * i] = 1;
if (i % prime[j] == 0)//if(i%prime[j]==0)break;
prime数组 中的素数是递增的,当 i 能整除 prime[j],那么 i*prime[j+1] 这个合数肯定被 prime[j] 乘以某个数筛掉。
因为i中含有prime[j], prime[j] 比 prime[j+1] 小。
break;
}
}
return;
}
int main()
{
cin >> n;
getprime();
for (ll i = 0; i < n; i++)
{
cout << prime[i] << " ";
}
return 0;
}