Visible Lattice Points[POJ3090]

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

传送门

思路

详细参考

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;
const int N=1010;
const int inf=1000;
int prime[N],m,phi[N];bool v[N];
int f[N];
void g_p()
{
	memset(v,false,sizeof(v));m=0;
	for(int i=2;i<=inf;i++)
	{
		if(!v[i])prime[++m]=i,phi[i]=i-1;
		for(int j=1;j<=m&&prime[j]*i<=inf;j++)
		{
			v[i*prime[j]]=1;
			if(i%prime[j]==0){phi[i*prime[j]]=phi[i]*prime[j];break;}
			else phi[i*prime[j]]=phi[i]*(prime[j]-1);
		}
	}
}
int main()
{
	g_p();f[1]=0;
	for(int i=2;i<=inf;i++)f[i]=f[i-1]+phi[i]*2;
	int t;scanf("%d",&t);
	for(int i=1;i<=t;i++)
	{
		int x;scanf("%d",&x);
		printf("%d %d %d\n",i,x,f[x]+3);
	}
	return 0;
}