本题输入范围为[1,1000001],要求计算出小于任意H number的HSemiPrime的个数,易知可以对输入范围进行打表。因为题目只在Hnumber的范围内考虑输入输出,所以首先就可以考虑优化:h = 4*k+1,0<=k<=250000。
不管是HPrime还是HSemiPrime都要考虑两个Hnumber相乘的情况,假设(4*n+1)(4*m+1) = 4*q+1,易推出4mn+m+n = q (1)
,也即是说两个Hnumber相乘必定是另一个Hnumber,且它们对应的k有上述关系。
因为HSemiNumber是两个HPrime的乘积,可以首先利用筛选法筛选出其中的Hprime对应的k:类似于素数筛选的方法,初始化k = [1,250000]所对应的Hnumber全为Hprime。再根据等式一所确定关系,逐步减少范围内的Hprime的个数。其中可类似于素数中只选择sqrt(x)作为最大因子的做法,因为若大于此数的数作为因子,则另一个因子也即是小于sqrt(x),则在前面已经枚举过,所以筛选法中m的枚举范围为[1,sqrt((q-2)/4)] = [1,250]。
当HPrime的打表以后,然后就可以用同样的方法求出其中的HSemiPrime。
#include<cstdio>
#include<cstring>
using namespace std;
#define maxK 250001
#define maxH 89075
char IsHPrime[maxK];
int HSemiPrimeNum[maxK];
int HPrimeNum;
int HPrimeList[maxH];
int main()
{
freopen("C:\\Users\\ifuding\\Desktop\\in.txt","r",stdin);
int i,j,k,n,num;
memset(IsHPrime,1,sizeof(IsHPrime));
memset(HSemiPrimeNum,0,sizeof(HSemiPrimeNum));
/***筛选出[1,1000001]之间的Hprime***/
for(i = 1;i <= 250;i++)
{
if(IsHPrime[i])
{
k = ((i*i)<<2)+(i<<1);
while(k < maxK)
{
IsHPrime[k] = 0;
k += (i<<2)+1;
}
}
}
HPrimeNum = 0;
for(i = 1;i < maxK;i++)
{
if(IsHPrime[i])
HPrimeList[HPrimeNum++] = i;
}
memset(IsHPrime,0,sizeof(IsHPrime));
/***筛选出[1,1000001]之间的HSemiPrime***/
for(i = 0;i < HPrimeNum&&HPrimeList[i] <= 250;i++)
{
j = HPrimeList[i];
n = i;
k = 4*j*j+j+j;
while(k < maxK)
{
IsHPrime[k] = 1;
n++;
k = 4*j*HPrimeList[n]+j+HPrimeList[n];
}
}
for(i = 1;i < maxK;i++)
{
HSemiPrimeNum[i] = HSemiPrimeNum[i-1];
if(IsHPrime[i])
HSemiPrimeNum[i]++;
}
while(~scanf("%d",&num)&&num)
printf("%d %d\n",num,HSemiPrimeNum[num/4]);
return 0;
}