poj 3292 Semi-prime H-numbers

#include<stdio.h>char n[250010]={0};int num[250010];void f(){int i,j;for(i=1;i<=251;++i){if(!n[i]){/*a,b是素數，a*b=c，則c是和數，故(4m+1)，(4n+1)是素數，(4q+1)是和數(4m+1)*(4n+1)=4q+116mn+4m+n4+1=4q+14mn+m+n=q為其下標關係*/for(int j=i;((((i*j)<<2)+i+j))<=250000;++j)n[((i*j)<<2)+i+j]=1;}}for(i=1;i<=251;i++){if(n[i]==0){for(j=i;((((i*j)<<2)+i+j))<=250000;j++){if(n[j]==0)n[((i*j)<<2)+i+j]=2;}}}num[1]=0;for(i=2;i<=250001;i++)if(n[i]==2)num[i]=num[i-1]+1;else num[i]=num[i-1];}int main(){f();int a,asd;while(scanf("%d",&a),a){asd=(a-1)/4;printf("%d %d\n",a,num[asd]);}return 0;}

該題和找素素數、和數類似，用篩選法

1.篩選出H-primes；

2.篩選出H-semi-prime

3.遍曆，算出各下標對應的H-semi-prime個數。