该题是一个很经典的欧拉函数的题,中间用到了线性筛求欧拉函数和一定的数学知识。
题目链接:
AcWing 201. 可见的点
对于该题,先理解为什么能用欧拉函数进行求解:
- 先看横纵坐标的关系,例如:(6,4) 这个坐标中 6 和 4 有 2 这个公约数,那么我们就可以知道它俩前面有 (3,2) 这个坐标,所以如果一个点 ( a , b ) (a,b) (a,b) ,如果 a , b a,b a,b 互质的话,那么 ( a , b ) (a,b) (a,b) 前面就一定没有与这个点在一个线上的点。
对于该题的思路:
- 对于该题的观察可以看出,左上角和右下角是相对称的,那么我们只考虑进行计算右下角的,那么结果乘 2 就好了。
- 我们从左往右扫描,看该横坐标 x x x 在 $1 \sim x-1 中有多个数字与 x x x 互质,那么求得的结果就是横坐标对应每一列的点的个数。
- 该题的关键就是线性筛求欧拉函数。如果用试除法直接求每个 1 ∼ N 1 \sim N 1∼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;
}