bzoj4659LCM

总感觉我的做法是假的莫比乌斯反演。。

只用到了莫比乌斯函数的性质

复杂度是O(nlogn)的，这题正好p是2^30可以用int自然溢出，不用取模，但是最后的答案还是要+p在%p，防止是负数

#include <bits/stdc++.h>
using namespace std;
const int P = 1<<30;
#define N 4000005
int T,m,n,cnt;
bool vis[N];
int sum[N],miu[N];
int prime[N];
int ans;
void Miu()
{
	miu[1] = 1;
	cnt = 0;
	for (int i=2;i<=N-5;i++)
	{
		if (!vis[i]) {prime[++cnt] = i;miu[i] = -1;}
		for (int j=1;j<=cnt && prime[j] * i <= N-5;j++)
		{
			vis[i * prime[j]] = 1;
			if (i % prime[j] == 0) break;
			else miu[i * prime[j]] = -miu[i];
		}
	}
}
void Sq()
{
	for (int i=1;i<=N-5;i++)
		if (miu[i] != 0)
			for (int j=i;j<=N-5;j+=i)
				sum[j] = sum[j] + j *(j/i)*miu[j/i];
	for (int i=1;i<=N-5;i++)
		sum[i] = sum[i-1] + sum[i];
}
int Sum(int x)
{
	return (x+1) * x / 2;
}
int main()
{
	Miu();
	Sq();
	scanf("%d", &T);
	while (T--)
	{
		scanf("%d%d", &n, &m);
		if (n > m) swap(n,m);
		ans = 0;
		for (int l=1,r;l<=n;l = r+1)
		{
			r = min(n / (n / l), m / (m / l));
			ans = ans + Sum(n/l) * Sum(m/l) * (sum[r] - sum[l-1]);
		}
		printf("%d\n", (ans % P + P) % P);
	}
}