线性筛素数法代码:
#include<iostream>
using namespace std;
const long N = 200000;
long prime[N] = {0},num_prime = 0;
int isNotPrime[N] = {1, 1};
int main()
{
for(long i = 2 ; i < N ; i ++)
{
if(! isNotPrime[i])
prime[num_prime ++]=i;
//关键处1
for(long j = 0 ; j < num_prime && i * prime[j] < N ; j ++)
{
isNotPrime[i * prime[j]] = 1;
if( !(i % prime[j] ) ) //关键处2
break;
}
}
return 0;
}
思想就是不断i乘上一个质数将其筛掉。
稍作修改,可以舍弃求模运算,并且可以只考虑奇数,并且可以同时求出质因子中的最小质数与最大质数。当然下面的maxprime数组对于求质数的过程没有用,可以删掉。
也可以在maxprime数组处稍作变动,同时求得一个数是不是唯一的质因子的幂次方。时间复杂度依然还是原先的O(n),并且没有求模运算。
const int MaxPrimeNum = 78499, MaxSize = 1001, MaxSize1 = 31251;
int MinPrime[MaxSize], MaxPrime[MaxSize]; //分别记录可以整除MaxSize的最小与最大prime的index //MaxPrime可以改成bit类型数组来记录是否是prime或者是否是Prime幂次方
int prime[MaxPrimeNum];
//自己修改,欧拉筛法求n以内的所有素数,存在prime数组里,index从1到num_prime-1 。
//并同时求出n以内的所有数的最小质因子与最大质因子,分别存在MaxPrime数组与MinPrime数组里
//复杂度O(n)
void GetPrimeSolution::euler_sieve(int n)
{
int i, j, len, num_prime = 1, tmpMaxIndex, tmp;
i = 2;
prime[len = tmpMaxIndex = MaxPrime[i] = MinPrime[i] = num_prime++] = i;
for (i = 3; i < n; i += 2)
{
if (!MinPrime[i])
{
prime[len = tmpMaxIndex = MaxPrime[i] = MinPrime[i] = num_prime++] = i;
}
else
{
tmpMaxIndex = MaxPrime[i];
len = MinPrime[i];
}
for (j = 2; j <= len && (tmp = i * prime[j]) < n; ++j)
{
MaxPrime[tmp] = tmpMaxIndex;
MinPrime[tmp] = j;
}
}
/* for (i = 0; i < num_prime; ++i)
{
printf("%d\t", prime[i]);
}
for (i = 0; i < n; ++i)
{
//printf("%d\t%d\t%d\n", i, MinPrime[i], MaxPrime[i]);
printf("%d\t%d\t%d\n", i, prime[MinPrime[i]], prime[MaxPrime[i]]);
}*/
}
const int MaxPrimeNum = 78499, MaxSize = 1001;
int MinPrime[MaxSize], MaxPrime[MaxSize]; //分别记录可以整除MaxSize的最小与最大prime的index //MaxPrime可以改成bit类型数组来记录是否是prime或者是否是Prime幂次方
int prime[MaxPrimeNum];
//自己修改,欧拉筛法求n以内的所有素数,存在prime数组里,index从1到num_prime-1 。
//并同时求出n以内的所有数的最小质因子与最大质因子,分别存在MaxPrime数组与MinPrime数组里
//复杂度O(n)
void euler_sieve(int n)
{
int i, j, len, num_prime = 1, tmp;
i = 2;
prime[len = MinPrime[i] = num_prime++] = i;
for (i = 3; i < n; i += 2)
{
if (!MinPrime[i])
{
prime[len = MinPrime[i] = num_prime++] = i;
}
else
{
len = MinPrime[i];
}
for (j = 2; j <= len && (tmp = i * prime[j]) < n; ++j)
{
MinPrime[tmp] = j;
}
}
for (i = 0; i < num_prime; ++i)
{
printf("%d\t", prime[i]);
}
/*
for (i = 0; i < n; ++i)
{
//printf("%d\t%d\t%d\n", i, MinPrime[i], MaxPrime[i]);
printf("%d\t%d\t%d\n", i, prime[MinPrime[i]], prime[MaxPrime[i]]);
}*/
}
int main()
{
euler_sieve(100);
return 0;
}