NOJ 1453 筛法

2016.10.10

【题目描述】
1453.筛法(Sieve Method)

时限：1000ms 内存限制：10000K 总时限：3000ms

描述
用筛法求［a，b］中的素数。

输入
2个正整数：a b。
a、b均在1000以内，且a小于等于b。
输出
[a b]区间内的所有素数，每个单独一行。

输入样例
2 5

输出样例
2
3
5

【解题思路】
判素数，常规题。
最普通的做法就是枚举小于sqrt(n)的每一个数，看他们是否能被n整除。
一般的筛法的原理是判断当前数是素数，把这个数的倍数都排除掉
线性筛法在此基础上继续优化
①如果 i 都是是素数的话，一个大的素数 i 乘以不大于 i 的素数，这样筛除的数跟之前的是不会重复的。筛出的数都是 N=p1*p2的形式, p1，p2之间不相等
②如果 i 是合数，此时 i 可以表示成递增素数相乘 i=p1*p2*…*pn, pi都是素数（2<=i<=n）， pi<=pj ( i<=j )p1是最小的系数。当p1==prime[j] 的时候，筛除就终止了，也就是说，只能筛出不大于p1的质数*i。

可以证明：
一个数会不会被重复筛除。
合数肯定会被干掉。

【代码实现】
一般的筛法：

# include <stdio.h>

int main(void)
{
    int a, b;
    int i, j;
    int p[1005];

    scanf("%d%d", &a, &b);

    for (i = 2; i <= b; ++i)
        p[i] = 1;
    p[1] = 0;

    for (i = 2; i <= b; ++i)
        if(p[i])
            for (j = i*i; j <= b; j += i) //注意i*i哦
                p[j] = 0;

    for (i = a; i <= b; ++i)
        if (p[i])
            printf("%d\n", i);

    return 0;
}

线性筛法：

# include <stdio.h>

int main(void)
{
    int i, j, num, n = 1000;
    int np[1005] = {1,1};
    int p[1005];

    num = 0;
    for (i = 2; i <= n; ++i)
    {
        if (!np[i])
            p[++num] = i;//如果当前数是素数，加到素数表里 

        for (j = 1; j <= num && i * p[j] <= n; ++j) 
        {
            np[i * p[j]] = 1;//i*素数表中的素数筛去合数 
            if (!(i % p[j]))//当p[j]为i的因子的时候退出循环，说      
                            //明只乘以小于i的质数 
                break;
        }
    }

    return 0;
}

【心得体会】