inline void getprime(int n) { for(register int i=1; i<=n; i++) prime[i] = 1 ; prime[1] = 0 ; for(register int i=2; i<=n; i++) { if(!prime[i]) continue; for(register int j=2; j<=n/i; j++) prime[i*j] = 0 ; } }
这个的算法是 O(n log log n)html
两只log 两只log 跑的快~算法
接近线筛 线筛太麻烦了。。主要是这个算法比较简便 线筛不会重复筛选spa
而这个有一部分是要重复筛选的code
下面来看真的线筛。。htm
“最小质因数 × 最大因数(非本身) = 这个合数”blog
的途径删掉。因为每一个数只被筛一次,时间复杂度为 O(n)。get
void GetPrime(int n)//筛到n
{ memset(isPrime, 1, sizeof(isPrime)); //以“每一个数都是素数”为初始状态,逐个删去
isPrime[1] = 0;//1不是素数
for(int i = 2; i <= n; i++) { if(isPrime[i])//没筛掉
Prime[++cnt] = i; //i成为下一个素数
for(int j = 1; j <= cnt && i*Prime[j] <= n/*不超上限*/; j++) { //从Prime[1],即最小质数2开始,逐个枚举已知的质数,并指望Prime[j]是(i*Prime[j])的最小质因数 //固然,i确定比Prime[j]大,由于Prime[j]是在i以前得出的
isPrime[ i*Prime[j] ] = 0; if(i % Prime[j] == 0)//i中也含有Prime[j]这个因子
break; //重要步骤。见原理
} } }