poj解题报告——3292

        经典的数论啊，说明一下题意：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]);
    }
}