#include<cstring>
using namespace std;

bool vis[1000000];    //存儲某個數(下標)是不是素數，是則為true，否則為false
int prime[1000000];     //存儲前n個素數，註意不是不超過n的素數，是前n 個。


void init_prim() {
    memset(vis, true, sizeof(vis));
    vis[0] = vis[1] = false;     //0和1不是素數
    int num = 0;
    for (int i = 2; i <= n; ++i) {
        if (vis[i] == true) {
            num++;
            prime[num] = i;
        }
        for (int j = 1; ((j <= num) && (i * prime[j] <= n));  ++j) {
            vis[i * prime[j]] = false;
            if (i % prime[j] == 0)
                break; //點睛之筆
        }
    }
}
 

素數打表高效算法