求素数筛法总结

当使用一般算法 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;               
            }
        }
    }
} 
  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值