[笔记] 线性筛质数（素数）

原理 ：保证n只会被n的最小质因子筛掉。
操作 ：
枚举i:(2~n)

1. 若当前i还没被筛掉，i就是质数。
2. 然后往后枚举primes[j]，筛掉primes[j] * i。
   因为：
   在找到i的最小质因子primes[j]前，
   由于是从小到大枚举的质数，说明此时primes[j]比i的最小质因子还小。
   此时primes[j]一定是primes[j] * i的最小质因子。
   可以保证primes[j] * i被它的最小质因数primes[j]筛掉。
3. i % primes[j] == 0时，不再往后搜primes[j+x]。
   因为：
   此时primes[j]是i的最小质因子。
   已知：primes[j] < iprimes[j+x]且primes[j]是i的最小质因子。
   可知：primes[j]才是primes[j+x] * i的最小质因子。
   primes[j+x] * i将会在后面枚举到i = primes[j+x] * i/primes[j]时，被primes[j]筛掉。
   可以保证primes[j+x] * i被它的最小质因数primes[j]筛掉，

补充

1. i的最小质因子primes[j]是第一个满足i % primes[j] == 0的primes[j]。

时间：数据范围在 10⁷以上的时候,线性筛法比埃氏筛法快一倍。

---

用代码说明：

//线性筛素数-注释版
//质数数组primes，计数器cnt记录下一个质数的位置
int primes[N], cnt;
//st数组记录筛除状态，0为没被筛掉，1为已被筛掉。
bool st[N];
void get_primes(int n ){
	//寻找不大于n的正整数中的质数
    for(int i = 2; i <= n; i ++){
    	//没被筛掉说明是质数,将这个新的质数加入primes里
        if(!st[i]) 
        	primes[cnt ++] = i;
        //从小到大枚举所有已知的质数 primes[j]
        for(int j = 0; primes[j] * i <= n; j ++){	
        	//当质数大于n / i的时候break;
       		//primes[j] * i <= n就是筛掉所有小于n的合数就可以了
       		//筛掉合数 i*primes[j]
            st[primes[j] * i] = 1;	
            //当这句话发生的时候,primes[j]一定是i的最小质因子
            //j再往后(设为j+x)就不能保证primes[j+x]是primes[j]*i的最小质因数了
            //因为i的质因数包含有primes[j],primes[j]是i的最小质因数
            //那么也将会是primes[j]*i的最小质因数
            //j最大会到i就结束。
            if(i % primes[j] == 0) 
            	break;
        }
    }
}

---

精简模板：

int primes[N], cnt;
bool st[N];
void get_primes(int n ){
    for(int i = 2; i <= n; i ++){
        if(!st[i]) 
        	primes[cnt ++] = i;
        for(int j = 0; primes[j] * i <= n; j ++){	
            st[primes[j] * i] = 1;	
            if(i % primes[j] == 0) 
            	break;
        }
    }
}