[POJ3292]Semi-prime H-numbers

原题链接

又是一道英文题
总之就
上翻译器吧

大体解释一下原题
形如4n+1的数叫做H数
乘法在H数组成的集合内是封闭的
H数中只能被1和本身整除的数叫H素数
剩下的叫H合数
其中能且仅能分解成两个H素数的H合数叫做H合成数

给你若干个数h
求在0~h范围内的H合成数的个数

利用筛选法求出H素数
[若i是H素数，那么5i+4i*x一定是H数但不是H素数]

然后两两相乘
处理出H合成数

预处理出每个范围内的H合成数的数量
直接输出

#include<iostream>#include<cstring>#include<algorithm>#include<cstdio>#include<cmath>#include<queue>#include<vector>#include<climits>#include<string>#include<cstdlib>#include<ctime>#define MAX 1000001#define LL long longusing namespace std;int i,j,Hprime[MAX+10],Hprime_num[MAX+10],Hmix[MAX+10],Hnum[MAX+10],n,tot;int main(){    for(i=5;i<=MAX;i+=4)    {        if(Hprime[i]) continue;        tot++;        Hprime_num[tot]=i;        for(j=i*5;j<=MAX;j+=i*4)            Hprime[j]=1;    }    for(i=1;i<=tot;i++)        for(j=1;j<=i&&Hprime_num[i]*Hprime_num[j]<=MAX;j++)        {            int t=Hprime_num[i]*Hprime_num[j];            Hmix[t]=1;        }    for(i=1;i<=MAX;i++)        Hnum[i]=Hnum[i-1]+Hmix[i];    scanf("%d",&n);    while(n)    {        printf("%d %d\n",n,Hnum[n]);        scanf("%d",&n);    }       return 0;}