poj3090--Visible Lattice Points（欧拉函数）

题目地址

✌，it‘s me

题目大意

给出一个n, 找到从坐标系原点（0， 0）出发并且和整数坐标连线，要求除了两个端点外线段上不能有整数坐标的点。问这样的线段有多少条。

解题思路

首先我们知道，gcd(x₂-x₁, y₂-y₁)就是(x₁, y₁)到(x₂, y₂)连线上有多少个整数点。
本题是源点，所以直接就是找到所有gcd(x, y) == 1的点有多少个。容易知道gcd(x, y) = gcd(y, x)，但是对于坐标来说，(x, y)在x!=y的情况下是并不相等的，正好本题的x, y一定是不相等的。所以求得的结果sum还要乘上2.
知道了大致思路后就要开始考虑边界问题了，对于（1， 0）， （0， 1），（1，1），这三个点与原点的连线是符合条件的，但是他们不符合gcd(x, y) == 1，所以要单独计算。所以结果为sum*2+3。
那么sum我们要如何去计算呢？
暴力去找到所有gcd(x, y) = 1的点肯定是太费时间了。我们可以想到使用欧拉函数的线性筛来求解本问题，本题数据最大为1000，可以先把1000以内的欧拉函数值求出来保存下来，每次直接使用。如果欧拉函数还没学请先看一下我的这篇博客，可以学习一下。

AC代码

#include <iostream>

using namespace std;

int prime[2000];
int phi[2000];
int vis[2000];
int ans = 0;

void Euler(int n)
{	
	phi[1] = 1;
	
	for (int i=2; i<=n; i++)
	{
		if (!vis[i])
		{
			prime[ans++] = i;
			phi[i] = i-1;
		}
		
		for (int j=0; j<ans && i*prime[j] <= n; j++)
		{
			vis[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);
			}
		}
	}
}

long long Sum(int n)
{
	long long sum = 0;
	
	for (int i=2; i<=n; i++)
	{
		sum+=phi[i];
	}
	
	return sum * 2 + 3;
}

int main()
{
	ios::sync_with_stdio(false);
	
	int t, num = 1;
	
	Euler(1000);
	
	cin >> t;
	
	while (t--)
	{
		int n;
		
		cin >> n;
		
		cout << num++ << " " << n << " " << Sum(n) << endl;
	}
	
	return 0;
}