题意:一个H-number是所有的模四余一的数,如 1,5,9,13,17,21...
H-primes数是H-number数(1除外),且它的H-number因子除了1只有它本身,如5,9,13,17,21...
但65是H-number数,却不是H-primes数,因为 65=5*13.
H-semi-prime是H-number数,且等于2个H-primes的乘积.如65
给你一个数n,问1到n有多少个H-semi-prime数
分析:用筛选法的思想,将H-primes筛选出来,同时标记在范围内两个H-primes的
乘积,最后再统计满足条件的个数
#include<stdio.h>
#define M 1000000
int a[M+5]={0}; //0表示为 H-primes数
void semi_prime()
{
long long i,j,t=0; //防止i*j乘法溢出,定义为long long
for(i=5;i<=M+1;i+=4)
for(j=i;j<=M+1;j+=4){
if(i*j>M+1)
break;
if(a[i]==0&&a[j]==0)
a[i*j]=1; //1表示为 H-semi-prime数
else
a[i*j]=-1; //-1表示 不是 H-primes数
}
for(i=5;i<=M+1;i++){ //统计1-i满足条件的个数
if(a[i]==1)
t++;
a[i]=t;
}
}
int main()
{
int n;
semi_prime();
while(scanf("%d",&n)!=EOF){
if(n==0)
break;
printf("%d %d\n",n,a[n]);
}
return 0;
}