POJ 3292 Semi-prime H-numbers

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

题意 ：表示这个题我真的没看懂。。。大意是讲：对于4*n+1(n为自然数)组成的集合H-numbers，乘法在该集合中是封闭的,在该集合内，只能被1和自身整除的为H-primes，然后由两个 H-primes的乘积为H-semi-prime，给你一个数k，让你求从1到k有多少个H-semi-prime。

思路 ：一开始没看懂题也就不知道意思，其实看来这道题就是枚举嘛。

#include <iostream>
#include <string.h>

using namespace std;

const int maxn = 1000201 ;
int h_prime[maxn] ;
int h_number[maxn] ;

void prime()
{
int cnt ;
for(int i = 5 ; i < maxn ; i += 4)
for(int j = 5 ; j < maxn ; j += 4)
{
cnt = i*j ;
if(cnt > maxn)
break ;
if(h_prime[i] == 0&&h_prime[j] == 0)
h_prime[cnt] = 1 ;
else h_prime[cnt] = 2 ;
}
}
int main()
{
memset(h_prime,0,sizeof(h_prime)) ;
prime() ;
int cnt = 0 ;
for(int i = 1 ; i < maxn ; i++)
{
if(h_prime[i] == 1)
cnt++ ;
h_prime[i] = cnt ;
}
int k ;
while(cin>>k && k)
{
cout<<k<<" "<<h_prime[k]<<endl ;
}
return 0;
}

View Code