poj 3292 Semi-prime H-numbers 筛素数（简单题）

题意：
H-numbers是4n+1形式的数
H-prime是由H-number构成的素数（即因子只有1和一个H-number）
H-semi是由两个H-prime的积构成
要求找出n（n<1e6）以内有多少个H-semi数？
分析：
简单题，模仿一下Euler筛法，得到H-prime表，然后O(n)去求H-semi数即可。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=1e6+9;
bool np[N];
int p[N],sum[N],semi[N];
void init()
{
    memset(np,0,sizeof(np));
    memset(semi,0,sizeof(semi));
    int cnt=0;
    for(int i=5;i<N;i+=4){
        if(!np[i])p[cnt++]=i;
        for(int j=0;j<cnt;j++){
            if(i*p[j]>=N)break;
            np[i*p[j]]=1;
            if(i%p[j]==0)break;
        }
    }
    for(int i=0;i<cnt;i++)
    for(int j=i;j<cnt;j++){
        if((ll)p[i]*p[j]>=N)break;
        semi[p[i]*p[j]]=1;
    }
    sum[0]=0;
    for(int i=1;i<N;i++)sum[i]=sum[i-1]+semi[i];
}

int main()
{
    //freopen("f.txt","r",stdin);
    int n;
    init();
    while(~scanf("%d",&n)&&n){
        printf("%d %d\n",n,sum[n]);
    }
    return 0;
}