POJ 3292 Semi-prime H-numbers

解题思路:

1 首先判断prime：H-numbers, 下标i表示4*i+1,j表示4*j+1,则(4*i+1)*(4*j+1)=4(4*i*j+i+j)+1,即下表i,j的乘积为下标4*i*j+i+j，判断prime方法类似判断普通prime方法，Line 10~12

2 判断Semi-prime，prime以及prime与no-prime的乘积肯定为非Semi-prime，Line 15~22

3 统计

#include <iostream>
 using namespace std;
#define MAXN 250001
bool IsPri[MAXN],IsSPri[MAXN];
int col[MAXN],pri[200];
int main()
{
    int n,i,j,t,p,ans;
    memset(IsPri,0x01,sizeof(IsPri));
    for(IsPri[0]=p=0,i=1;i<251;i++)
        if(IsPri[i])
            for(pri[p++]=i,j=i;(4*i*j+i+j)<MAXN;j++)
                IsPri[4*i*j+i+j]=false;
    memset(col, 0, sizeof(col));
    for(i=0;i<MAXN;i++)IsSPri[i]=IsPri[i]?0:1;
    for(IsSPri[0]=i=0;i<p;i++)
        for(t=0,j=pri[i]+1;t<MAXN;j++)
        {
            if(IsPri[j])continue;
            t=j*pri[i]*4+j+pri[i];
            if(t<MAXN)IsSPri[t]=0;
        }
    for(i=1;i<MAXN;i++)
        if(IsSPri[i])col[i]=col[i-1]+1;
        else col[i]=col[i-1];
    while(scanf("%d", &n)&&n)
        printf("%d %d\n",n, col[(n-1)/4]);
    return 0;
}