筛法求素数

判断一个数是不是素数，最简单的方法

#include<stdio.h>
int is_sushu(int x)//是素数就返回1，否则返回0
{
    if(x==1)//1既不是素数也不是合数
        return 0;
    else if(x==2)//2是素数
        return 1;
    else
    {
        for(int i=2; i<=sqrt(x); i++)//sqrt降低时间复杂度
            if(x%i==0)//如果x能被i整除，说明x不是素数
                return 0;
        return 1;//是素数返回1
    }
}

但如果要判断的不是一个数，而是成千上万个数，上述方法显然太慢了。
于是引入了筛法的概念。
筛法时间复杂度为O(x*ln（lnx)）

用筛法求素数的基本思想是：把从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

代码实现

#include<stdio.h>
#include<string.h>
#include<math.h>
#define N 10000//N的大小可以根据需要变化
int a[N];//利用数组的下标记录要判断的数字
void is_sushu()
{
    memset(a,0,sizeof(a));//对数组a进行初始化为0，不是素数的标记为1，剩下为0的就是素数了
    a[1]=1;//1既不是素数也不是合数，先标记为0
    for(int i=2; i<=sqrt(N); i++)
    {
        if(a[i]==0)//如果i是素数
        {
            for(int j=2; j*i<=N; j++) //循环标记的范围是i*j<N
            {
                a[i*j]=1;//如果i是素数，那么i*j肯定不是素数
            }
        }
    }//最终所有非素数都标记为1，素数都标记为0
}

还有一种算法，对前面的筛法进行了优化，是湖南大学huicpc39同学设计的（可能是学来的，也可能是自创的，反正很腻害）。
在数量级更大的情况下就可以发现一般筛法和优化后的筛法有明显区别。
很可惜，源码让我给丢了，等找到再贴上来