线性筛的原理

线性筛

初始时，假设全部都是素数，当找到一个素数时，显然这个素数乘上另外一个数之后都是合数
把这些合数都筛掉，即算法名字的由来。但仔细分析能发现，这种方法会造成重复筛除合数，影响效率。

例如 30=5*3*2，它被2*15筛了一次，又被3*10筛了一次，所以就有了如下算法：

int Mark[MAXSIZE];  
int prime[MAXSIZE];  
  
//判断是否是一个素数  Mark 标记数组 index 素数个数  
int Prime(){  
    int index = 0;  
    for(int i = 2; i < MAXSIZE; i++){  
        //如果未标记则得到一个素数  
        if(Mark[i] == 0) prime[++index] = i;  
        //标记目前得到的素数的i倍为非素数  
        for(int j = 1; j <= index && prime[j] * i < MAXSIZE; j++){  
            Mark[i * prime[j]] = 1;  
            if(i % prime[j] == 0) break;  
        }  
    }  
    return index;  
}

利用了每个合数必有一个最小素因子。每个合数仅被它的最小素因子筛去正好一次。所以为线性时间。
代码中体现在：
if(i%prime[j]==0)break;
prime数组中的素数是递增的,当 i 能整除 prime[j]，那么 i*prime[j+1] 这个合数肯定也可以被 prime[j] 筛掉，因为i中含有prime[j], prime[j] 比 prime[j+1] 小。

接下去的素数同理，所以不用筛下去了。
在满足i%prme[j]==0这个条件之前以及第一次满足该条件时,pr[j]必定是pr[j]*i的最小因子

[End]