11.24 作业 Problem D: 筛选素数

Description

在数学上，素数的分布没有任何已知规律，因此检测一个数是否素数，只能用比它小的素数来检测整除性质。如果要求出一定范围内的素数表直接检测的代价就太高了。一般采用筛选法的思想：
把从1开始的、某一范围内的正整数从小到大顺序排列，1不是素数，首先把它筛掉。剩下的数中选择最小的数是素数，然后去掉它的倍数。依次类推，直到筛子为空时结束。
如有：
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
1不是素数，去掉。剩下的数中2最小，是素数，去掉2的倍数，余下的数是：
3 5 7 9 11 13 15 17 19 21 23 25 27 29
剩下的数中3最小，是素数，去掉3的倍数……如此下去直到所有的数都被筛完，求出的素数为：
2 3 5 7 11 13 17 19 23 29
用筛选的办法求出素数表之后，再进行查询。

求区间m～n之间的素数，其中m,n为整数，且0<=m<=n<=500000。

Input

输入多行，到EOF结束。每行为一个测试用例，两个整数m和n，满足0<=m<=n<=500000。总测试用例数不超过1000个。

Output

输出为m～n之间的所有素数，每个数一行，两个测试用例之间用一个空行分隔。

如果m～n之间没有素数，则输出一个空行。

思路：

与上次的3n+1相似将该删掉的数储存为-1 防止读取注意输出格式

注：有个研究了半天的素数判断函数效率比较高本题在oj上运行的速度为200毫秒

代码：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#define MAX_LEN 500000

int sushu(int n)//应该是我所知的最快的判断素数的函数
{
    int sqrtn=(int)sqrt((double)n);
    int i;
    if(n==2)
        return 1;
    if(n%2==0||n==1)
        return 0;
    else
    {
    for(i=3; i<=sqrtn; i+=2)
    {
        if(n%i==0)
        return 0;
    }
    }
    return 1;
}
int main()
{
    int x,y,count;
    while(scanf("%d%d",&x,&y) != EOF)
    {
        if(x > y)//比较x y大小
        {
            int temp = y;
            y = x;
            x = temp;
        }
        int i,str[MAX_LEN];
        memset(str,0,sizeof(str));
        for(i = x; i <= y; i++)
        {
            if(str[i] < 0)//跳过删掉的数
                continue;
            if(sushu(i) == 1)
            {
                int flag;
                for(flag = i+1; flag <= y; flag++)//将后面所有的素数倍数删掉
                {
                    if(flag % i == 0)
                        str[flag] = -1;
                }
                str[i] = 1;//将该素数储存
            }
        }
        count = 0;
        for(i = x; i <= y; i++)
            if(str[i] == 1)
            {
                printf("%d\n",i);
                count++;
            }
        if(count == 0)
            printf("\n");

        printf("\n");//注意输出格式
    }

    return 0;
}