高效素数打表法

在网上看到的一种方法， 觉得很不错所以自己敲了一份记一下。。

比较省时间。。

#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
int prime[100010];
bool is_prime[10001000];
int main()
{
    int i, j, k;
    memset(is_prime, true, sizeof(is_prime));
    is_prime[1] = false;
    for(i = 2, k = 0; i <= 10001000; i++)
    {
        if(is_prime[i])
        {
            prime[k++] = i;
        }
        for(j = 0;i * prime[j] <= 10001000;j++)
        {
            is_prime[i * prime[j]] = false;
            if(i % prime[j] == 0)
                break;
        }
    }
    return 0;
}