[BZOJ2705] [SDOI2012] Longge的问题 [莫比乌斯反演&欧拉函数]

[$\mathfrak{L}\mathfrak{i}\mathfrak{n}\mathfrak{k}$]

---

求 $\sum (i,n)$ ， $n\le 2^{32}$ 。
所有小于 $n$ 的数都可以对答案产生贡献，考虑利用 $gcd$ 缩小数据范围
那么按照 $gcd$ 分类？

$(i,n)=1$ 的数有 $\varphi (n)$ 个。

$(i,n)=k$ 的数有几个？显然就是 $(i/k,n/k)=1$ 的数的个数 $\varphi (n/k)$
但是只是这样还不够，因为我们要算的是贡献（
所以要求的应该是 $k\cdot \varphi (n/k)$

于是只需要求 $\sum _{d|n}(n/d)\cdot \varphi (d)$
把 $\varphi (d)$ 拆开继续化式子还可以进一步降低复杂度。

---

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cmath>
using namespace std;
unsigned long long n, sq, ans, xsq;
unsigned long long phi(unsigned long long x)
{
	unsigned long long ret = x;
	xsq = 1ull * sqrt(1.0 * x);
	for (unsigned long long i = 2ull; i <= xsq; ++i)
	{
		if (x % i == 0) ret = ret / i * (i - 1);
		while (x % i == 0) x /= i;
	}
	return (x > 1ull) ? (ret / x * (x - 1ull)) : ret;
}
int main()
{
	scanf("%lld", &n);
	sq = 1ull * sqrt(1.0 * n);
	for (unsigned long long i = 1ull; i <= sq; ++i)
	{
		if (n % i == 0) 
		{
			ans += (n / i) * phi(i);
			if (i * i != n) ans += i * phi(n / i);
		}
	}
	printf("%lld", ans);
	return 0;
}