[BZOJ 4659] Lcm

BZOJ传送门

题目描述

给出A，B，考虑所有满足$l\le a\le A$，$l\le b\le B$，且不存在$n>1$使得$n^2$同时整除$a$和$b$的有序数,对$(a,b)$,求其$lcm(a,b)$之和。答案模$2^{30}$。

输入输出格式

输入格式

第一行一个整数T表示数据组数。接下来T行每行两个整数A，B表示一组数据。

$T\le 2000$,$A,B\le 4\times 10^6$

输出格式

对每组数据输出一行一个整数表示答案模$2^{30}$的值

输入输出样例

输入样例#1

5
2 2
4 6
3 4
5 1
23333 33333

输出样例#1

7
148
48
15
451085813

解题分析

其实这道题不考虑那个因数没有平方数的条件的话和 **这道题**是一样的， 有这样一个式子：
$$\sum _{T=1}^n\frac{(1+\lfloor \frac{n}{T}\rfloor )\times \lfloor \frac{n}{T}\rfloor }{2}\times \frac{(1+\lfloor \frac{m}{T}\rfloor )\times \lfloor \frac{m}{T}\rfloor }{2}T\sum _{k|T}k\times \mu (k)$$
然后我们发现我们枚举的$T$似乎就是$gcd(i,j)$？？ 于是$O(Nln(N))$预处理，$O(T\sqrt{N})$回答即可。

是否因数中含有完全平方数用$\mu (i)$判定即可。

代码如下：

#include <cstdio>
#include <cctype>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>
#define R register
#define IN inline
#define gc getchar()
#define W while
#define MX 4000500
const int MOD = 1 << 30;
#define ll long long
int miu[MX], pri[MX], pcnt, n, m, T;
ll F[MX];
bool npr[MX];
void get()
{
	miu[1] = 1;
	R int i, j, tar;
	for (i = 2; i < MX; ++i)
	{
		if(!npr[i]) pri[++pcnt] = i, miu[i] = -1;
		for (j = 1; j <= pcnt; ++j)
		{
			tar = i * pri[j];
			if(tar > MX) break;
			npr[tar] = true;
			if(!(i % pri[j])) break;
			miu[tar] = -miu[i];
		}
	}
	for (i = 1; i < MX; ++i) if(miu[i])
	{
		for (j = i; j < MX; j += i)
		(F[j] += miu[j / i] * j % MOD * (j / i) % MOD) %= MOD;
	}
	for (i = 1; i < MX; ++i) (F[i] += F[i - 1]) %= MOD;
}
IN int getsum(R int val) {return (val + 1) * val >> 1;}
int main(void)
{
	scanf("%d", &T);
	get();
	ll ans;
	W (T--)
	{
		scanf("%d%d", &n, &m);
		if(n > m) std::swap(n, m); ans = 0;
		for (R int lef = 1, rig; lef <= n; lef = rig + 1)
		{
			rig = std::min(n / (n / lef), m / (m / lef));
			(ans += (F[rig] - F[lef - 1]) * getsum(n / lef) % MOD * getsum(m / lef) % MOD) %= MOD;
		}
		printf("%lld\n", (ans % MOD + MOD) % MOD);
	}
}