HAOI2011 Problem B

题目链接：https://www.luogu.com.cn/problem/P2522
这是一道莫比乌斯反演模板题

我们要求的是:$$\sum _{i=a}^b\sum _{i=c}^d[gcd(i,j)=k]$$
利用二维前缀和，我们发现只要求出$$\sum _{i=1}^{n/k}\sum _{j=1}^{m/k}[gcd(i,j)=1]$$

运用莫比乌斯反演的知识，原式等价于：$$\sum _{i=1}^{n/k}\sum _{j=1}^{m/k}\sum _{d|gcd(i,j)}\mu (d)$$

将这个和式变换位置：$$\sum _{d=1}\mu (d)\sum _{i=1}^{n/k}\sum _{j=1}^{m/k}1$$
$$=\sum _{d=1}\mu (d)[\frac{n}{kd}][\frac{m}{kd}]$$
数论分块即可.

$Code$

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#define int long long
const int inf = 1e9 + 7;
using namespace std;

signed main(){
	freopen ("std.in", "r", stdin);
	freopen ("std.out", "w", stdout);
	int n, k;
	scanf("%lld%lld", &n, &k);
	long long ans = n * k;
	for (register int l = 1, r; l <= n; l = r + 1){
		if (k / l != 0)	r = k / (k / l);
		else    r = inf;
		r = min(r, n);
		ans -= (l + r) * (r - l + 1) / 2 * (k / l);
	}
	cout << ans << endl;
	return 0;
}