【poj 3292】Semi-prime H-numbers

题目链接：
http://poj.org/problem?id=3292
题意：
定义一种数叫H-numbers，它是所有能除以四余一的数。

在H-numbers中分三种数：

1、H-primes，这种数只能被1和它本身整除，不能被其他的H-number整除，例如9是一个H-number，能被1，3，9整除，但3不是H-number，所以他是H-primes。

2、H-semi-primes是由两个H-primes相乘得出的。

3、剩下的是H-composite。

问给一个数，求1到这个数之间有多少个H-semi-primes。
题解：
注意他所定义的H-primes，并不是普通的素数，是由特殊的素数（23333），所以可以用特殊的素数筛法，直接一边将答案筛出来，简单说就是打表。
代码：

#include<iostream>
#include<algorithm>
#include<stdio.h>
#define maxn (1000500)
using namespace std;
int n,sum,vis[maxn],pri[maxn];
//vis[i] 为0表示是H-primes 为1表示为H-semi-primes 为-1表示是H-composite
int main()
{
    for (int i=5;i<=maxn;i+=4)
    for (int j=5;j<=maxn;j+=4)
    if (i*j>maxn) break;
    else if (vis[i]==0&&vis[j]==0)
    vis[i*j]=1;
    else vis[i*j]=-1;
    for (int i=1;i<=maxn;i++)
    {
        if (vis[i]==1) sum++;
        pri[i]=sum;
    }
    while(scanf("%d",&n))
    {
        if (n==0) return 0;
        printf("%d %d\n",n,pri[n]);
    }
}