素数线筛-欧拉筛法为什么i%primes[j] == 0时要break

最新推荐文章于 2022-11-12 18:34:16 发布

绫零依

最新推荐文章于 2022-11-12 18:34:16 发布

阅读量423

点赞数 2

分类专栏：算法文章标签：算法素数筛

本文链接：https://blog.csdn.net/Qiao712/article/details/112795572

版权

算法专栏收录该内容

1 篇文章 0 订阅

订阅专栏

//素数线筛-欧拉筛法
int primes[100000];	//用来储存已获得的素数（当栈用）
int num = 0;		//已得到的素数的个数  （栈顶）
bool getPrime(int n){
    for(int i = 0; i<100000; i++) isNotPrime[i] = 0;
    isNotPrime[1] = true;

    for(int i = 2; i<=n; i++){
        if(!isNotPrime[i]) primes[num++] = i;
        for(int j = 0; j<num; j++){     //用已知的素数乘当前的数i，筛掉后面的
            if(i * primes[j] >= n) break;
            isNotPrime[i * primes[j]] = true;
            if(i % primes[j] == 0) break; 	/*why*/
        }
    }
}

目的就是为了保证每个合数不被重复筛掉，这是核心

原因：

很明显，当i % primes[j] == 0时，primes[j]是i的最小的素因数。
若此时不break的话，在后面的循环（内层j的循环）中筛掉的 x = i*primes[k](k>j) 一定含有这个素因数primes[j]。而拿除掉这个最小的素因数primes[j]得到的x/primes[j] 一定是比i大的 （∵primes[k]/primes[j]>1）
所以在后面的循环(外层i的循环）中x一定又会被某个更大的i筛掉（这个i就是 x/primes[j])
所以要在i % primes[j] == 0时，break，保证后面不会再筛一次