POJ-3292 Semi-prime H-numbers

虽然这题被分到了同余模运算，但是貌似用不到，打个表就行。

开始理解错了题意，以为Hprime是素数，打了素数表，后来发现 只要不被Hnumber整除就好。

而且先是开了Prime[maxn]，RE了orz，改成了Prime[maxn+10]过了。

代码如下：

#include <algorithm>
#include <iostream>
#include <cstring>
#include <string.h>
#include <string>
#include <cstdio>
#include <vector>
#include <cmath>
#include <stack>
#include <map>
#include <set>
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1e6+500;
int Prime[maxn+10];
int ans[maxn];
void Init()
{
    memset(Prime,0,sizeof(Prime));
    for(int i=5;i<=maxn;i+=4)
    {
        for(int j=5;j<=maxn;j+=4)
        {
            int cul=i*j;
            if(cul>maxn)
                break;
            if(Prime[i]==0&&Prime[j]==0)
                Prime[cul]=1;
            else Prime[cul]=-1;
        }
    }
    int cnt=0;
    for(int i=1;i<maxn;i++)
    {
        if(Prime[i]==1)
            cnt++;
        ans[i]=cnt;
    }
    return ;
}

int main()
{
    ios::sync_with_stdio(false);
    Init();
    int h;
    while(scanf("%d",&h)!=EOF&&h)
    {
        printf("%d %d\n",h,ans[h]);
    }
    return 0;
}