poj3292

链接：点击打开链接

题意：可以写成4i+1的数为H-number，只能分解成1*本身，不能分解成其他的H-number的数为H-primes，能够恰好分解成两个H-primes的乘积，且只能是两个数的乘积的数为H-semi-primes，输出1-N之间有多少个H-semi-primes

代码：

#include <set>
#include <cstdio>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
bool sign[1000005];
int num[1000005];
set<int> s;
set<int>::iterator ite;
int main(){
    int i,j,k,n,ans;
    k=0;
    memset(sign,0,sizeof(sign));
    for(i=5;i<=200005;i+=4)
    for(j=5;j<=200005;j+=4){
    if(i*j>1000001)
    break;
    sign[i*j]=1;
    }                           //先把H-number的乘积都标记上
    for(i=5;i<=200005;i+=4)
    for(j=5;j<=200005;j+=4){
        if(i*j>1000001)
        break;
        if(!sign[i]&&!sign[j]){ //记录两个H-primes的乘积
        num[k++]=i*j;
        }
    }
    sort(num,num+k);
    k=unique(num,num+k)-num;
//    cout<<k<<endl;
//    for(i=0;i<100;i++)
//    cout<<num[i]<<" ";
//    cout<<endl;
    while(~scanf("%d",&n)!=EOF&&n){
        ans=lower_bound(num,num+k,n)-num;
        if(num[ans]==n)         //二分出最后答案
        printf("%d %d\n",n,ans+1);
        else
        printf("%d %d\n",n,ans);
    }
    return 0;
}