线性筛法（一）--素数筛法（二）

换一种筛子

上一篇我们讲到用素数的倍数一定不是素数这一个方法来构建筛子。不清楚的童鞋可以再看看。

但是这种方法还是有问题：比如筛30时，在 i=2 的时候， k = 2 * 15 筛了一次；在 i=5， k = 5 * 6 的时候又筛了一次。

今天我们换一种思路，并添加一个优化，从而解决这个问题，还是求1-n的素数。
核心思想：所有的合数可以表示为两个或多个素数的积。（实际上也相当于合数为素数乘一个比它小的合数-。-）

实现：我们遍历每个数，把它与之前找到的素数相乘，把相乘的结果标记为合数。
——实际上是把上一篇的“某个素数乘以每个数”改为了“每个数乘以某个素数”，但是由于我们做了一点优化，使前面提到的重复情况去除，所以算法得到了改进。

1，初始化（a数组标记，queue数组存素数与素数乘积）

a 0 1 2 3 4 5 6 7 - 0 0 0 0 0 0 0 0 queue 0 1 2 3 4 5 6 7 - 0 0 0 0 0 0 0 0

2，我们先找到2，存入数组queue；

a 0 1 2 3 4 5 6 7 - 0 0 0 0 0 0 0 0 queue 0 1 2 3 4 5 6 7 - 2 0 0 0 0 0 0 0

3，进行素数相乘 2x2，标记 4；

a 0 1 2 3 4 5 6 7 - 0 0 0 0 1 0 0 0 queue 0 1 2 3 4 5 6 7 - 2 0 0 0 0 0 0 0

4，找到3，并对3进行操作。。。；

a 0 1 2 3 4 5 6 7 8 9 10 - 0 0 0 0 1 0 1 0 0 1 0 queue 0 1 2 3 4 5 6 7 - 2 3 0 0 0 0 0 0

5，找到4，但不放入queue，并对 4 处理

a 0 1 2 3 4 5 6 7 8 9 10 11 12 - 0 0 0 0 1 0 1 0 1 1 0 0 1 queue 0 1 2 3 4 5 6 7 - 2 3 0 0 0 0 0 0

6，如此操作。。。

我们来看看程序：

#include <iostream>#include <stdio.h>#include <string.h>using namespace std;const int N = 100001;int main(){    int n;    int a[N];    int queue[N];    while(scanf("%d",&n) == 1)    {        memset(a,0,sizeof(a));        memset(queue,0,sizeof(queue));        int tail = 0;        for(int i = 2; i < n; i++)        {            if(a[i] == 0)                queue[tail++] = i;            for(int j = 0; j < tail && i * queue[j] < N; j++)            {                a[i * queue[j]] = 1;                if(!(i % queue[j]))//这是一个优化，去除15*2，6*5的重复                    break;            }        }        printf("2");        for(int i = 3; i < n; i++)        {            if(a[i] == 0)                printf(" %d",a[i]);        }    }    return 0;}

上一篇：线性筛法（一）–素数筛法（一）