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
);
}
return 0;
}