素数
素数最重要的时素数的求法,可以根据定义求还有它的优化等,其实都用不大上(复杂度太高),这里不啰嗦了直接上埃筛:
#include<cstdio>
#include<iostream>
using namespace std;
int N;
bool prime[10010];
void init(){
for(int i = 2; i < N; i ++) prime[i] = true;
for(int i = 2; i*1 < N; i ++){//判断改成i*i<N,解决sqrt精度问题
if(prime[i]){
for(int j = i*i; j < N; j += i){//从i*i开始就可以了,因为i*(小于i的数)之前已经算过了
prime[j] = false;
cout<<j<<" ";
}
cout<<endl;
}
}
}
int main(){
cin>>N;
init();
}
乍一看可能有点懵,细看还是很好理解的,其实就是把素数的倍数剪去,实在不理解课以把N开小点跟着走一遍就好了(额(⊙﹏⊙) 傻傻的我就是这么干的),后面的好多地方都用到这种思想。复杂度O(nloglogn)。
乃个其实埃筛就是让你缓冲一下,我们还有O(n)的做法,线性筛:
#include<cstdio>
const int N = 100000 + 5;
bool prime[N];//prime[i]表示i是不是质数
int p[N], tot;//p[N]用来存质数
void init(){
for(int i = 2; i < N; i ++) prime[i] = true;//初始化为质数
for(int i = 2; i < N; i++){
if(prime[i]) p[tot ++] = i;//把质数存起来
for(int j = 0; j < tot && i * p[j] < N; j++){
prime[i * p[j]] = false;
if(i % p[j] == 0) break;//保证每个合数被它最小的质因数筛去
}
}
}
int main(){
init();
}
细看会发现这里不止把质数的整数倍剪去还极大的减少了埃筛第二重循环的复杂度,在N特别大的时候尤为明显。
昨天再队友那搞来一个判断一个数是否为素数的很快的模板,感觉比较有用,补上:
int sushu(long long num)
{
if(num == 1) return 0;
if(num ==2|| num==3 ) return 1 ;
if(num %6!= 1&&num %6!= 5) return 0 ;
int tmp =sqrt( num);
for(int i= 5;i <=tmp; i+=6 )
if(num %i== 0||num %(i+ 2)==0 )return 0 ;
return 1 ;
}
再补充一个数论定理:
1~n内的质数个数:nln(n)(实际偏小,越趋近于+∞越准确,估计复杂度是够的)