素数筛法

注释很详细 直接上代码了 唯一分解定理是重点

#include<iostream>
#include<algorithm>
#include<cstring>
#define maxn 1000
#define ll long long
using namespace std;
int prim[maxn],tot;
int flag[maxn<<2];//flag要开比打表范围大一点 否则会有数组越界
int cnt=0;
void prime(){//maxn为打表范围 
	ll i,j;
	for( i=2;i<=maxn;i++){
	
		if(!flag[i]){
		prim[cnt++]=i;
		}//如果没有被标记 将素数存入prim数组 
		for(j=0;j<cnt;j++){
			if(i*prim[j]>maxn)break;//唯一分解定理 
			flag[i*prim[j]]=1;//当它不大于maxn范围  标记这个合数

			if(i%prim[j]==0){//如果当前i和素数有因子 退出（唯一分解定理  ex：20%2==0 如果不退出的话 会标记20*3  但是在后面的打表范围中30*2 会再次标记 造成时间损失） 
				break;
			}
		}
	}	
}
int main(){
prime();
for(int i=0;i<cnt;i++){
	cout<<prim[i]<<" ";
}		
	return 0;
}