POJ 3090 欧拉函数

#include <cstdio>
#include <cstring>
const int maxn = 1010;
int phi[maxn], T, n, kase;
void init()
{
	memset(phi, 0, sizeof(phi));
	for (int i = 2; i < maxn; i++)
		if (!phi[i])
			for (int j = i; j < maxn; j += i)
			{
				if (!phi[j]) phi[j] = j;
				phi[j] = phi[j] / i * (i - 1);
			}
}
int main(int argc, char const *argv[])
{
	init();
	scanf("%d", &T);
	while (T--)
	{
		scanf("%d", &n);
		int ans = 3;
		for (int i = 2; i <= n; i++)
			ans += 2 * (phi[i]);
		printf("%d %d %d\n", ++kase, n, ans);
	}
	return 0;
}

大小为n的图（有（n+1）^2个格点），求出图中可视点的个数，欧拉函数不考虑a=b和为0的情况，所以特殊考虑。