POJ3292 Semi-prime H-numbers

定义H数为4n+1,n>=0的数，H素数为H数中只能拆成1*x，x为H数的数，求1-n<=1000000中H素数的个数。

筛法。

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<algorithm>
#include<iostream>
using namespace std;

int n;
#define maxn 1000011
int prime[maxn],lp=0;bool notp[maxn];
void makeprime()
{
    memset(notp,0,sizeof(notp));
    notp[1]=0;
    for (int i=5;i<=n;i+=4)
        if (!notp[i])
        {
            prime[++lp]=i;
            for (int j=i*5;j<=n;j+=(i<<2)) notp[j]=1;
        }
}
bool isans[maxn];int ans[maxn];
int main()
{
    n=1000001;
    makeprime();
    memset(isans,0,sizeof(isans));
    for (int i=1;i<=lp;i++)
        for (int j=1;prime[j]*prime[i]<=n;j++)
            isans[prime[i]*prime[j]]=1;
    ans[0]=0;
    for (int i=1;i<=n;i++) ans[i]=ans[i-1]+isans[i];
    while (scanf("%d",&n) && n)
        printf("%d %d
",n,ans[n]);
    return 0;
}