POJ 3292 删数

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define M 1000000
typedef long long ll;
ll n;
ll p[M];

void init()
{
	memset(p,0,sizeof(p));
	ll i,j;
	ll tot=0;
	for(i=5;i<M;i+=4)
	{
		for(j=5;i*j<M;j+=4)
		{
			if(p[i]==0 && p[j]==0)
			{
				p[i*j]=1;
			}
			else
			{
				p[i*j]=-1;
			}
		}
	}
	ll sum=0;
	for(i=1;i<M;i++)
	{
		if(p[i]==1)
		{
			sum++;
		}
		p[i]=sum;
	}
}
int main()
{
	init();
	while(scanf("%lld",&n),n)
	{
		printf("%lld %lld\n",n,p[n]);
	}
	return 0;
}