牛半仙的妹子gcd【数学】

>Link

牛客1022普及T2

---

>Description
给出一个整数$n$，求出${\sum }_{i=1}^n{\sum }_{j=1}^n{\sum }_{k=1}^{n}gcd(i,j,k)$

---

>解题思路
数据有点水比赛的时候打了一个$n^3$的优化模拟都过了（我是菜鸡==）

正解如下。
将三个∑简化成两个，已知$gcd(i,j,k)=gcd(gcd(i,j),k)$
第一个$n^2$循环：枚举$i$和$j$，用一个桶累计相同的gcd出现的次数
第二个$n^2$循环：枚举不同的$gcd(i,j)$和$k$，累加上预处理的出现个数$\ast$三个数的$gcd$

---

>代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 1005
using namespace std;

int n, a[N], ans;

int gcd (int a, int b)
{
	if (!b) return a;
	return gcd (b, a % b);
}

int main()
{
	scanf ("%d", &n);
	for (int i = 1; i <= n; i++)
	  for (int j = 1; j <= n; j++)
	    a[gcd(i, j)]++;
	for (int G = 1; G <= n; G++)
	  for (int k = 1; k <= n; k++)
	    ans += a[G] * gcd (G, k);
	printf ("%d", ans);
	return 0;
}