AcWing 201. 可见的点(线性筛求欧拉函数)

该题是一个很经典的欧拉函数的题，中间用到了线性筛求欧拉函数和一定的数学知识。
题目链接：
AcWing 201. 可见的点

对于该题，先理解为什么能用欧拉函数进行求解：

• 先看横纵坐标的关系，例如：(6,4) 这个坐标中 6 和 4 有 2 这个公约数，那么我们就可以知道它俩前面有 (3,2) 这个坐标，所以如果一个点 $(a,b)$ ,如果 $a,b$ 互质的话，那么 $(a,b)$ 前面就一定没有与这个点在一个线上的点。

对于该题的思路：

1. 对于该题的观察可以看出，左上角和右下角是相对称的，那么我们只考虑进行计算右下角的，那么结果乘 2 就好了。
2. 我们从左往右扫描，看该横坐标 $x$ 在 $1 \sim x-1 中有多个数字与 $x$ 互质，那么求得的结果就是横坐标对应每一列的点的个数。
3. 该题的关键就是线性筛求欧拉函数。如果用试除法直接求每个 $1\sim N$ 之间的值，数据一定会超时的。

关于欧拉筛求欧拉值我也写在另一个博客中，不懂的前去观看：
线性筛求欧拉函数

代码实现：

#include<iostream>
using namespace std;
const int N = 1005;
int primes[N],used[N],cnt;
int phis[N];

void get_phi(int x)
{
    cnt =0;
    phis[1] = 1;
    for(int i = 2; i <= x; i++)
    {
        if(!used[i])
        {
            primes[++cnt] = i;
            phis[i] = i-1;
        }
        for(int j = 1; i * primes[j] <= x; j++)
        {
            used[i * primes[j]] = 1;
            if(i % primes[j] == 0)
            {
                phis[i * primes[j]] = primes[j] * phis[i];
                break;
            }
            else phis[i * primes[j]] = (primes[j]-1) * phis[i];
        }
    }
}

int main (void)
{
    get_phi(1000);
    int T;
    cin >> T;
    for(int cases = 1; cases <= T; cases++)
    {
        int n;
        cin >> n;
        int res = 3;
        for(int i = 2; i <= n; i++)
        res += phis[i] * 2;
        cout << cases << " " << n << " " << res << endl;
    }
    return 0;
}