题目大意
有如下几个定义:
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[n] << endl;
}
}