经典的数论啊,说明一下题意:H-number是4*n+1这样的数,如1,5,9,13...。H-primes是这样一个H-number:它只能唯一分解成1*它本身,而不能表示为其他两个H-number的乘积。一个H-semi-prime是一个这样的H-number:它正好能表示成两个H-primes的乘积(除了1*它本身),这种表示法可以不唯一,但它不能表示为3个或者以上H-primes的乘积。现在给出一个数n,要求区间[1,n]内有多少个H-semi-prime。
所以,直接打表会快不少,表在init()函数里。
代码如下
#include<stdio.h>
#include<stdlib.h>#define size 1000001
int H_Number[size+1];
void init()
{
int i,j,k,m,p;
memset(H_Number,0,sizeof(H_Number));
for(i=5;i<=size;i+=4)
{
for(j=5;j<=size;j+=4)
{
m=i*j;
if(m>size)
break;
if(H_Number[i]==0&&H_Number[j]==0)
H_Number[m]=1;
else
H_Number[m]=-1;
}
}
p=0;
for(k=1;k<=size;k++)
{
if(H_Number[k]==1)
p++;
H_Number[k]=p;
}
return;
}
void main()
{
int h;
init();
while(scanf("%d",&h)!=EOF&&h!=0)
{
printf("%d %d\n",h,H_Number[h]);
}
}