11105 - Semi-prime H-numbers （生成特定的H素数表)

记MAXN 为最大的H范围 最大的H数n为(MAXN-1)/4;

则对所有的H数枚举其H数倍削掉其倍数可得到H素数表，然后枚举生成H半素数即可

#include <cstdio>
#include <vector>
#include <cmath>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
const int maxn = 1000001;
const int lim  = (maxn-1)/4;
int a[maxn+10],vis[lim+10];
vector<int> primes;
int get_prime(){
    memset(vis,0,sizeof(vis));
    int m = sqrt(lim+1);
    for(int i=1;i<=m;i++)if(!vis[i])
       for(int j=i;(4*i+1)*(4*j+1)<=maxn;j++){
                vis[4*i*j+i+j]=1;
    }
    for(int i=1;i<=lim;i++){
          if(!vis[i]) primes.push_back(4*i+1);
    }
    memset(a,0,sizeof(a));
    for(int i=0;i<primes.size();i++)
       for(int j=0;j<primes.size();j++){
            if(primes[i]*primes[j] > maxn) break;
            a[primes[i]*primes[j]]=1;
    }
    for(int i=1;i<=maxn;i++) a[i]+=a[i-1];
}
int main()
{
    get_prime();
    int n;
    while(scanf("%d",&n)==1&&n){
           printf("%d %d\n",n,a[n]);
    }
    return 0;
}