POJ 3292 Semi-prime H-numbers(数筛)

题目大意：H-number指的是4*n+1的数（除4余1）。H-prime为没有H-number的因子的数,H-semi-prime为两个H-prime组成。求解1-h中的semi-prime个数。

思路：我们可以像素数筛一样来晒数，起初都是H-prime所以他们的乘积都是semi-prime。然后标记，如果有标记的数那么一定是H-composites,标记为不是H-semi-prime和H-prime。然后遍历（注意枚举的时候范围一定要大，否则因为可能是i最大的数不够大和j相乘得不到想要的数，然后控制停止）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<map>
#include<queue>
#include<cmath>
#include<algorithm>
#define MAX 1000
#define inf 0x3f3f3f3f
using namespace std;
int s[1000510],p[1000500];
int main(){
    int n,i,j;
    memset(s,0,sizeof(s));
    for(i=5;i<=1000100;i+=4){
        for(j=5;j<=100100;j+=4){
            if(i*j>1000010)break;
            if(!s[i]&&!s[j]){
                s[i*j]=1;
            }
            else{
                s[i*j]=-1;
            }
        }
    }
    int ans=0;
    for(i=1;i<=1000010;i++){
        if(s[i]==1){
            ans++;
        }
        p[i]=ans;
    }
    while(~scanf("%d",&n)&&n){
        printf("%d %d\n",n,p[n]);
    }
    return 0;
}