poj 3292 Semi-prime H-numbers （筛选法思想）

链接：poj 3292

题意：一个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;
}