H数

如果i是H素数，那么5i+4i*x一定是H数且不是H素数
预处理处所有的H素数就可以啦
题目

#include <cstdio>
#include <iostream>
using namespace std;
const int H_MAX=1000001+16;
int f_prime[H_MAX],mix_prime[H_MAX];
int H_prime[H_MAX];
int cnt[H_MAX];
int tot;
void cl()
{
    int j,i;
    for(i=5;i<H_MAX;i+=4)
    {
        if(f_prime[i]) continue;
        H_prime[++tot]=i;
        for(j=i*5;j<H_MAX;j+=i*4)
         f_prime[j]=1;
    }

    for(i=1;i<=tot;i++)
     for(j=1;j<=i&&H_prime[i]*H_prime[j]<H_MAX;j++)
      mix_prime[H_prime[i]*H_prime[j]]=1;

    for(i=1;i<H_MAX;i++)
     cnt[i]=cnt[i-1]+mix_prime[i];
}
int main()
{
    cl();

    while(1)
    {
        int x;

        cin>>x;

        if(x) printf("%d %d\n",x,cnt[x]);
        else return 0;
    }
}