素数筛

借用网上某位大神的板子

#include 
//#include 
#include 
#include
using namespace std;

int prim[1000005]={2,3,5};
//素数是分为基本素数{2,3}、阳素数{6N+1，N>=1}形式的、阴素数{6N-1，N>=1}形式的
//为了代码的好写，在这里这样写的 :
//数除了{2,3,5}为素数，其他的数可以写成6N，6N+1,6N+2,6N+3,6N+4,6N+5  N>=1 可以表示全部的数
//6N，6N+2,6N+4都为偶数，不是素数，6N+3 == 3(2N+1) 不是素数，那么就只筛6N+1和6N+5就可以了
int main()
{
    int i, j, flag;
    int gcd = 2;
    int k = 3;

    for(i=7;i<=10000000;i+=gcd)
    {
        gcd = 6-gcd;//6N+1和6N+5的变换。
        flag = 1;
        for(j=0;prim[j]*prim[j]<=i;j++)//判断一个数是否为素数，只需判断这个数在2--sqrt(x)范围即可。
            if(i%prim[j] == 0){//因为一个开根号比乘法是要慢的，所以用乘法速度更快。
                flag = 0;
                break;
            }
        if(flag)//若这个数没有被其他素数整除 说明为素数。
            prim[k++] = i;
    }
        for(i=0;i<1000;i++)//输出前一千个素数
        printf("%d  ",prim[i]);
        return 0;
}