POJ - 3292 Semi-prime H-numbers (埃氏筛思想）

题意：H数为4n + 1这样的数，H数的乘法是封闭的，H素数只能表示为1和它本身的乘积，除了H素数，其余的都是H合数，半素数是恰好能表示为两个H素数的乘积，给出一个数X，输出这个数以及小于等于这个数的半素数的个数

思路：埃氏筛思想

#include <stdio.h>
const int maxn = 1000005;
int a[maxn],cnt[maxn];
void init()
{
    int i,j,ans;
    for(i = 5; i < maxn; i += 4) {
        for(j = 5; i * j < maxn; j += 4) {
            if(a[i] == 0 && a[j] == 0) a[i * j] = 1;
            else a[i * j] = -1;
        }
    }
    ans = 0;
    for(i = 1; i < maxn; i++) {
        if(a[i] == 1)
            ans++;
        cnt[i] = ans;
    }
}
int main(void)
{
    int n;
    init();
    while(scanf("%d",&n) != EOF && n)
        printf("%d %d\n",n,cnt[n]);
    return 0;
}