题意:
题解:
我们反过来,先求出 gcd(i,j) 是平方数的组合个数,因为感觉这个更好算,然后用 n*m 减去这些个数就是答案了。
设 f(x) 为 gcd(i,j)=x 的个数
那么,因为我们要算出全部的在 n 内的 gcd 等于平方数的个数。
就是说,我们要求出,这样子还是比较难算,我们交换一下求和的次序。
关于这一块的打表跟之前写的 BZOJ2818 有点像,同样用线性筛可以筛出结果 (当然也可以不用线性筛,两遍循环也同样可以不超时的筛出结果)。
① 当 d 是素数的时候,只有 x=1 符合,就是说 T(d)=μ(d)
② d%prim[j] != 0 意味着 d 中不包含 prim[j] 但是 d*prim[j] 中包含了,而莫比乌斯函数是个积性函数。而且 prim[j] 仅有一个。
③ d%prim[j] == 0 意味着 d 中包含了 pirm[j] 且大于 1 个
而莫比乌斯的值,只要 μ(x) 出现其中一个素因子的指数大于 1 ,那么 μ(x)=0 ,就是说这一项是等于 0 的。
能打表出这个结果,然后T(d)的前缀和,我们就可以用分块加速来算出最终结果了。
,
#pragma comment(linker, "/STACK:102400000,102400000")
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include