当我真正理解素数线性筛法

参考自:点击链接
主要代码:

const int MAXN = 10000010;  

bool com[MAXN];  
int primes, prime[MAXN/10]; 

void solve(int n)  
{  
    primes = 0;  
    memset(com,false,sizeof(com));  
    com[0] = com[1] = true;  
    for (int i = 2; i <= n; ++i)  
    {  
        if (!com[i])  
        {  
            prime[++primes] = i;  
        }  
        for (int j = 1; j <= primes && i*prime[j] <= n; ++j)  
        {  
            com[i*prime[j]] = true;  
            if (!(i % prime[j]))  
                break;  
        }  
    }  
}  
现在具体讲解一下证明:
最难理解的是:
if (!(i % prime[j]))  
                break;  
要从下面两个方面:
每个数至少被访问一次
每个数至多被访问一次
进行理解.
对于质数,一定会在i的循环中访问到,并确定为质数。
对于合数,一定可以分解为一个最小素因子和其他数的乘积。
比如:
合数 i = p(最小素因子)* a;
若 i%prime[k] ==0;
则 p * a * prime[k+1] 可以被后面的 a * prime[k+1] 再乘以素数 p 筛选出来,
    (显而p<prime[k+1]) 所以i%prime[k] == 0 时要停止。
  • 2
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值