Crash的数字表格 [莫比乌斯反演][神题]

传送门

发现有点棘手, 不妨把k放出来

是时候反演了

令

有

则

可以整除分块做, 只不过要预处理mu[k]*k^2

这里又可以整除分块...

---

#include<bits/stdc++.h>
#define N 10000050
#define LL long long
#define Mod 20101009
using namespace std;
int prim[N],isp[N],mu[N],tot;
LL val[N],ans; // mu[i] * i ^ 2
int n,m;
void prework(){
	mu[1] = val[1] = 1;
	for(int i=2;i<=N-50;i++){
		if(!isp[i]) prim[++tot] = i, mu[i] = -1;
		for(int j=1;j<=tot;j++){
			if(prim[j] * i > N - 50) break;
			isp[prim[j] * i] = 1;
			if(i % prim[j] == 0) break;
			mu[i * prim[j]] = -mu[i];
		}
	}
	for(int i=1;i<=N-50;i++){
		val[i] = (LL)mu[i] * (LL)i * (LL)i;
		val[i] %= Mod;
		val[i] += val[i-1]; 
		val[i] %= Mod;
	}
}
LL calc(int x,int y){
	LL ans = 0;
	if(x>y) swap(x,y);
	for(int l=1,r; l<=x; l=r+1){
		int v1 = x/l, v2 = y/l;
		r = min(x/v1, y/v2);
		LL tmp = ((((LL)v1 * (LL)(v1+1) / 2) % Mod) * (((LL)v2 * (LL)(v2+1) / 2) % Mod)) % Mod;
		ans += ((val[r] - val[l-1] + Mod) % Mod) * tmp; ans %= Mod;
	} return ans % Mod;
}
int main(){
	prework();
	scanf("%d%d",&n,&m);
	if(n>m) swap(n,m);
	for(int l=1,r; l<=n; l=r+1){
		int v1 = n/l, v2 = m/l;
		r = min(n/v1, m/v2);
		ans += (((LL)(l+r) * (LL)(r-l+1) / 2) % Mod) * calc(v1, v2);
		ans %= Mod;
	} printf("%lld",ans); return 0;
}