当使用一般算法 i->sqrt(n) 超时时, 就必须思考对整个程序的优化。
优化的地方在于降低重复率。即降低对2,4,6这种有共同因数的运算枚举,而只是对一个更小的循环周期范围(如6N算法,任何一个自然数都可以表述为6N,6N+1,6N+2,6N+3…6N+5,而显然当N>=1时,6N,6N+2,6N+3,6N+4一定不为质数,即缩小了需要遍历的范围)甚至对素数组来求解。
埃式筛法
算法思想: 从2开始向下遍历,每找到一个素数,则从数组中剔除以这个素数为因子的数(仍具有重复性的浪费,如6会在检查2和检查3时被重复剔除)
#define MAXSIZE 100000
bool prime[MAX];
void isprime()
{
prime[2] = 0;
for(int i = 2;i <= floor(sqrt(max));i++) //只需要运行到最大的数
的根号的值即可如果想顺便输出质数的数量,则需要修改条件
{
if(!prime[i] || i != 2) // 将2的倍数从数组中剔除
{
prime[i] = i%2 == 0 ? 1 : 0 ;
}
if(!prime[i])
{
for(int j = i;j <= max;j++)
{
prime[j] = 1;
}
}
}
}
此处仍可对函数进行部分优化。
将最后函数内的循环改为
for(int j = i*i,j <= max; j += 2*n) //进一步降低了重复率.
欧拉筛
与埃式筛法的区别:
埃式筛法有可能对一个数据筛去多次,
而欧拉筛只用每个合数的最小质因数去筛。
#define MAX 1000000
int primes[MAX/3]; //大概..防止堆栈溢出
bool judge[MAX] = {false};
void isprime()
{
int pr_num = 0;
for(int i = 2;i <= MAX;i++)
{
if(!judge[i])
{
primes[pr_num] = i;
pr_num++;
}
for(int j = 0;i * primes[j] <= MAX;j++)
{
judge[i * primes[j]] = true;
if(i % primes[j] == 0) //始终保证由最小质数筛选
{
break;
}
}
}
}