POJ 3292 Semi-prime H-numbers

题目大意

有如下几个定义：
H-numbers为4n+1(n属于自然数)
H-primes必为H-numbers，且H-primes只能分解为以下两个H-numbers的乘积1 * H-primes == H-primes
H-semi-primes是刚好两个H-primes的乘积，并且为H-numbers
给出n<=1000001，要你求出1~n中有多少H-semi-primes

解题思路

纯英语的题目很容易漏过细节。一开始我以为H-primes为质数且为H-numbers，所以做这道题时花了不少功夫。
其次是，我很容易把题目写的比标解更复杂，不管标解是简单的还是复杂的，要自己想想怎样才能更接近标解的思想。
另外，标程中筛H-semi-primes的方法我认为不是很严谨。也许是我还没有证明出来，呵呵~

Code

#include <iostream>

using namespace std;

const long long MAXN = 1000002;

long long f[MAXN];

int main()
{
for (int i = 5; i < MAXN; i += 4)
{
for (int j = i; j < MAXN && (long long)i * (long long)j < MAXN; j += 4)
{
if (f[i] == 0 && f[j] == 0) f[i * j] = 1;
else f[i * j] = -1;

}
}
int cnt = 0;
for (int i = 1; i < MAXN; i++)
{
if (f[i] == 1) cnt++;
f[i] = cnt;
}
int n;
while (cin >> n)
{
if (!n) return 0;
cout << n << ' ' << f
<< endl;
}
}