P1447 [NOI2010] 能量采集（容斥原理）

题目链接：点击这里

题目大意：
给定 $n,m$ 求 ${\sum }_{i=1}^n{\sum }_{j=1}^{m}2\gcd (i,j)-1$

题目分析：
这个题可以通过反演求解，但我们换一种思路，从容斥原理的角度除法
设 $f_x$ 表示 $\gcd (i,j)=x$ 的数对个数，此时问题就转化成了求 $\sum (2x-1)f_x$ ，直接求这个函数不太好求我们考虑间接求该函数值
设 $g_x$ 表示 $x$ 是 $(i,j)$ 公因数的数对个数，那么显然 $g_x=(n/x)\times (m/x)$
有了 $g_x$ 后我们就可以通过容斥原理把 $x$ 的 $2,3,...$ 倍删掉，即 $f_x=g_x-{\sum }_{i=2}^{ix\le min(n,m)}{f}_i$
这样我们就可以倒着枚举求解 $f_x$ ，边求解边维护答案即可，时间复杂度为 $O(nlogn)$ （调和级数）

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<string>
#define int long long
using namespace std;
const int maxn = 1e5+5;
int n,m,ans,f[maxn];
signed main()
{
	cin>>n>>m;
	for(int i = n;i;i--)
	{
		f[i] = (n/i)*(m/i);
		for(int j = 2;i*j <= n;j++) f[i] -= f[i*j];
		ans += (2*i-1)*f[i];
	}
	cout<<ans<<endl;
	return 0;
}