201 可见的点（线性筛求欧拉函数）

1. 问题描述：

在一个平面直角坐标系的第一象限内，如果一个点 (x，y) 与原点 (0，0) 的连线中没有通过其他任何点，则称该点在原点处是可见的。例如，点 (4，2) 就是不可见的，因为它与原点的连线会通过点 (2，1)。部分可见点与原点的连线如下图所示：

编写一个程序，计算给定整数 N 的情况下，满足 0 ≤ x，y ≤ N 的可见点 (x，y) 的数量（可见点不包括原点）。

输入格式

第一行包含整数 C，表示共有 C 组测试数据。每组测试数据占一行，包含一个整数 N。

输出格式

每组测试数据的输出占据一行。应包括：测试数据的编号（从 1 开始），该组测试数据对应的 N 以及可见点的数量。同行数据之间用空格隔开。

数据范围

1 ≤ N，C ≤ 1000

输入样例：

4
2
4
5
231

输出样例：

1 2 5
2 4 13
3 5 21
4 231 32549
来源：https://www.acwing.com/problem/content/203/

2. 思路分析：

首先我们需要看一下什么样的点可以被照到，可以发现照到的点都是从原点出发在直线y = kx（1 <= x，y <= N）上的第一个点，若y = kx直线上的第一个点是(x0，y0)，可以发现(x0，y0)互为质数，所以(x0，y0)是y = kx直线上的第一个点等价于(x0，y0)互为质数，如果不互为质数那么（x0，y0）> 1，则（x0 / d，y0 / d）在直线上并且会将(x0，y0)挡住就矛盾了，所以(x0，y0)一定是互质的，并且由直线y = kx的上边与下边是对称的，所以我们只需要求一边最终将答案乘以2即可，以下边为例，可以发现求解的其实是欧拉函数φ(x)的值（欧拉函数是小于等于n的正整数中与n互质的数的数目），所以我们求解出1~n的欧拉函数的值即可，可以使用线性筛法求解欧拉函数，主要分为以下几种情况求解：

• 当i = 1的时候结合题目可以发现φ(1) = 1
• 当i为质数的时候那么φ(i) = i - 1
• 当pj是i的最小质因子的时候，φ(i * pj) = φ(i) * primes[j]
• 当pj不是i的最小质因子的时候，φ(i * pj) = φ(i) * (primes[j] - 1)

3. 代码如下：

from typing import List


class Solution:
    count = 0
    
    # 线性筛法求解欧拉函数
    def init(self, n: int, primes: List[int], st: List[int], phi: List[int]):
        phi[1] = 1
        for i in range(2, n):
            if st[i] == 0:
                primes[self.count] = i
                self.count += 1
                phi[i] = i - 1
            j = 0
            while i * primes[j] < n:
                st[i * primes[j]] = 1
                # pj是i的最小质因子
                if i % primes[j] == 0:
                    phi[i * primes[j]] = phi[i] * primes[j]
                    break
                # pj不是i的最小质因子
                phi[i * primes[j]] = phi[i] * (primes[j] - 1)
                j += 1

    def process(self):
        N = 1010
        primes, st, phi = [0] * N, [0] * N, [0] * N
        self.init(N, primes, st, phi)
        self.count = 0
        t = int(input())
        for i in range(1, t + 1):
            n = int(input())
            # 中间那个点
            res = 1
            for j in range(1, n + 1):
                res += phi[j] * 2
            print(i, n, res)


if __name__ == '__main__':
    Solution().process()