POJ - 3292 Semi-prime H-numbers（素数筛法）

http://poj.org/problem?id=3292

先求所有的H-pirme，所有的H-prime两两相乘打表。

注：

1、H-numbers是所有除以4余1的数，而H-prime则是只能在这些H-numbers中分解因式只得到1*本身的数（1除外）。所以9也是H-prime（原来理解错了）。

2、 H-semi-primes可能有两种不同的分解方法，如441=21*21=9*49 。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int MAX_N = 1000002, INF = 0x3f3f3f3f;
int N;
int prime[MAX_N];
bool isprime[MAX_N];
bool hp[MAX_N];
int a[MAX_N];
int ans;
void sieve(int n)
{
	int p = 0;
	ans = 0;
	for(int i = 0; i <= n; i++) isprime[i] = true;
	isprime[0] = isprime[1] = false;
	for(int i = 2; i <= n; i++)
		if(isprime[i] && i%4 == 1)
		{
			prime[p++] = i;
			for(int j = 2*i; j <= n; j+=i)
				isprime[j] = false;
		}
		for(int i = 0; i < p; i++)
			for(int j = i; j < p && (long long)prime[i]*prime[j]<=n; j++)
				if(prime[i]*prime[j] % 4 == 1)
					ans++, hp[prime[i]*prime[j]] = true;
		for(int i = 1; i <= n; i++)
			a[i] = a[i-1] + hp[i];
}
int main()
{
	sieve(MAX_N-1);
	while(scanf("%d", &N), N)
	{
		printf("%d %d\n", N, a[N]);
	}
	return 0;
}