POJ 3292 数论 类似素数筛

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
using namespace std;
const int Max=1000010;
int vis[1000010];
int sum[1000010];
int main()
{
int n,t;
memset(vis,0,sizeof(vis));
for(int i=5;i<=1001;i+=4)
{
if(!vis[i])
for(int j=5;(t=j*i)<=Max;j+=4)
{
if(!vis[j])                    <span id="transmark"></span>
vis[t]=1;
else                       //存在多个A-prime两两相乘的情况
vis[t]=-1;
}
}
sum[0]=0;
for(int i=1;i<1000010;i++)
sum[i]=sum[i-1]+(vis[i]==1);
while(~scanf("%d",&n),n)
printf("%d %d\n",n,sum
);
return 0;
}