poj 3292 Semi-prime H-numbers

http://poj.org/problem?id=3292

如果一个数 i 是 H-素数，那么5i+4i*x一定是H数但不是H-素数，
因为（5i+4i*x）mod 4=5i mod 4=(5 mod 4)(i mod 4)=1*1=1,且5i+4i*x=i(4x+5).

做法：先找出所有的H-素数，在两两相乘得出H-合成数

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define LL long long
#define M 1000001
using namespace std;
bool ish_p[M+10],ish_m[M+10];
int ans[M+10],hp[M+10],cnt=0;
int main()
{
    for(int i=5;i<=M;i+=4)
    {
        if(ish_p[i]) continue;//不是H-素数
        hp[++cnt]=i;
        for(int j=i*5;j<=M;j+=i*4)
         ish_p[j]=true;//5i+4i*x一定不是H-素数 
    }
    for(int i=1;i<=cnt;i++)
     for(int j=1;j<=i&&hp[i]*hp[j]<=M;j++)
      ish_m[hp[i]*hp[j]]=true;
    for(int i=1;i<=M;i++)
     ans[i]=ans[i-1]+ish_m[i];
    int h;
    while(1)
    {
        scanf("%d",&h);
        if(h==0) break;
        printf("%d %d\n",h,ans[h]);
    }
    return 0;
}