线性筛质数

问题：生成小于等于max的所有质数

传统判断质数的方法是，利用for循环检测该数是否含有除1以外的因子，实现代码如下：

bool isZhishu(int n)
{
	int i;
	if(n==0||n==1)
	{
		return false;
	}
	for(i=2;i*i<=n;++i)
	{
		if(n%i==0)
		{
			return false;
		}
	}
	return true;
}

该方法在判断单个数是否为质数时可行，但是如果用来筛选较大范围内的质数时，由于时间复杂度较大，非常容易运行超时，因此下面采用线性筛筛选质数，下面是具体实现代码：

void filter_prime(int max)
{
	int*prime,*is_prime;//is_prime[i]若为0，则i为质数；is_prime若为1，则i为合数
	is_prime=(int*)malloc((max+1)*sizeof(int));
    //因为is_prime[i]表示的是i的状态，所以构造时需构造max+1
	prime=(int*)malloc(max*sizeof(int));//prime用来储存筛选出来的质数
	for(i=0;i<max+1;++i)
	{
		is_prime[i]=0;//初始化is_prime数组为0
	}
	for(i=2;i<=max;++i)
	{
		if(is_prime[i]==0)
		{
			prime[now++]=i;
            //如果i是质数，存储起来
		}
		for(j=0;j<now&&i*prime[j]<=max;++j)
		{
			is_prime[i*prime[j]]=1;//将对应的合数状态改为1，即筛出去
			if(i%prime[j]==0)
            //这是线性筛最关键的一步，区别于其他筛法，这是实现线性的核心
			{
				break;
			}
		}
	}
}

总结来说线性筛法的核心在于**把每个合数用最小质因子筛去，而实现这一步的关键就在于if(i%prime[j]==0){break;}**。

就比如12这个合数分解为2x2x3，最小质因子是2，代码执行最终12会被这样筛去is_prime[6*prime[0]]=1，即12=6x2，它不会被12=4x3筛掉，因为判断到4的时候，由于4%prime[0]即2==0时，break了。

而至于为什么要break？为什么可以break？

第一个问题，因为每个合数只能被最小质因子筛去，不加break会造成一个合数被筛多次，增加时间复杂度。

第二个问题，因为若 i%prime[j]==0，说明 i 存在质因子prime[j]，也就是说之后 i*prime[j+1]，prime[j+2]等等产生的合数会有质因子prime[j]，之后的循环一定会被prime[j]筛掉，不用担心会漏筛某个合数。