POJ 3292 Semi-prime H-numbers

         题意：把形如 4n+1的数称作 H数 像 1 5 9 13····  类似素数的定义  定义了H素数 即 它的所有 H数约数 只有 1和这个数本身 例如 5 9 13 17 21 而 25 不是。又定义了semi-prime H数，如果一个数是semi-prime H数 那么 它最多只能拆成两个 H素数的积  例如 25=5*5  而 125=5*5*5 所以125不是。
        给一个数n，求出1~n之间有多少semi-prime H数

解法：  用筛法去找出H素数 然后去找semi-prime H数 注意 找semi-prime H数时 会有重复 如 441=9*49=21*21

代码如下：

/**

吉林大学
Jilin U

Author:     sinianluoye (JLU_LiChuang)
Date:       2014-07-26
Usage:      一水到底 C

**/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>

#define ll long long
#define eps 1e-8
#define ms(x,y) (memset(x,y,sizeof(x)))
#define fr(i,x,y) for(ll i=x;i<=y;i++)

using namespace std;
const ll maxn=1e6+10;

bool prime[maxn];
ll pri[maxn],cnt;
ll s[maxn];

void make()
{
cnt=0;
ms(prime,0);
for(ll i=1;(4*i+1)<maxn;i++)
prime[4*i+1]=1;
ll num=0;
for(ll i=0;(4*i+1)<maxn;i++)
s[num++]=i;
for(ll i=5;i<maxn;i+=4)
{
if(prime[i])
{
pri[cnt++]=i;
for(ll j=5*i,k=1;j<maxn;j=(s[++k]*i))
{
prime[j]=0;
}
}
}
}

ll ans[maxn];

void doit()
{
//cout<<"here"<<endl;

ll num=0;
for(ll i=0;i<cnt;i++)
for(ll j=0;j<=i&&((pri[i]*pri[j])<maxn);j++)
{
//cout<<i<<' '<<j<<endl;
ans[num++]=(pri[i]*pri[j]);
// if(num%10000==0)cout<<num<<endl;
}
sort(ans,ans+num);
int f=-1;
for(int i=0;i<num;i++)
{
if(ans[i]==f)ans[i]=maxn+1;
else f=ans[i];
}
sort(ans,ans+num);
}
int main()
{
make();
doit();
//cout<<cnt;
//for(ll i=0;i<100;i++)cout<<i<<' '<<ans[i]<<endl;
ll n;
while(cin>>n,n)
{
ll num;
for(num=0;ans[num]<=n;num++);
cout<<n<<' '<<num<<endl;
}
}
/*************copyright by sinianluoye (JLU_LiChuang)***********/