【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-compositeint 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]);    }}