牛半仙的妹子gcd

$$\mathrm{牛半仙的妹子gcd}$$

题目链接：$\mathrm{n}\mathrm{o}\mathrm{w}\mathrm{c}\mathrm{o}\mathrm{d}\mathrm{e}\mathrm{r}212919$

关于这场比赛

——>点我可以查看其它题目（目录）<——

到牛客看：

——>点我跳转<——

题目

牛半仙有 $n$ 个妹子，每个妹子有一个属性值，第 $i$ 个妹子的属性值为 $n-i+1$。
牛半仙认为 $k$ 个妹子 { b 1 , b 2 … b k } \{b_1,b_2…b_k\} {b1​,b2​…bk​} 的相同度为 $gcd(a_{b_1},a_{b_2}\dots a_{b_k})$。
牛半仙想知道 ${\sum }_{i=1}^n{\sum }_{j=1}^n{\sum }_{k=1}^{n}gcd(i,j,k)$ 的值。

输入

一行一个整数 $n$。

输出

一行一个整数，表示任意三个妹子的相同度之和。

样例输入

2

样例输出

9

数据范围

对于 $40\%$ 的数据，$n\le 200$
对于 $100\%$ 的数据，$n\le 1000$
$gcd(a_1,a_2\dots a_k)$ 表示 { a 1 , a 2 … a k } \{a_1,a_2…a_k\} {a1​,a2​…ak​} 的最大公约数。

思路

这道题有很多种做法，我选择了最简单的预处理（因为只会这个
还有莫比乌斯反演和欧拉两种做法。

预处理的主要思想就是先预处理出 gcd 中一边是 $i$，另一边从 $1$ 枚举到 $n$ 的 $j$，我们就是预处理数所有 $gcd(i,j)$ 的和。

那与处理完之后，我们就可以枚举题目的 $i$ 跟 $j$，因为前面预处理了，就可以省去枚举 $k$ 的部分。

时间复杂度是 $O(n^{2}logn)$，勉强能过。

比赛时

想了一会就想到了预处理的方法。

在看了一下复杂度应该可行就打完跑了。

代码

#include<cstdio>

using namespace std;

int n, a[1001], ans;

int gcd(int x, int y) {
	if (!y) return x;
	return gcd(y, x % y);
}

int main() {
	scanf("%d", &n);
	
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= n; j++)
			a[i] += gcd(i, j);//预处理出一对
	
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= n; j++)
			ans += a[gcd(i, j)];//统计
	
	printf("%d", ans);
	
	return 0;
}