poj3292 Semi-prime H-numbers

题意：

H-number 是一种4*n+1的数。

H-primes 是只可以被1和本身整除的H-numbers(在4n+1的数域内)。

H-semi-primes 是有且仅有两个H-numbers因子(除了1和本身)。

问给你一个n，求从1到n一共有多少个H-semi-primes数。

思路：

H[i*j]=H[i]+H[j]+1，H[i]=0表示i是H-primes，H[i]=1表示i是H-semi-primes。H[i]>1表示H[i]有多个只因子

#include<iostream>
using namespace std;
#define MAX 1000003
int  H[MAX];

int main(){
	int i,j;
	for(i=5;i<1000;i+=4){
		for(j=i;i*j<MAX;j+=4){
			H[i*j]=H[i]+H[j]+1;}
	}
	for(i=25;i<MAX;i++){
		H[i]=(H[i]==1)+H[i-1];
	}
	while(cin>>i&&i)
		cout<<i<<' '<<H[i]<<endl;
}