bzoj4659 Lcm （莫比乌斯反演）

题意

求 $\sum\limits_{i=1}^n\sum\limits_{j=1}^m lcm(i,j) (gcd(i,j) 无平方质因子)$

题解

无平方质因子，即 $\mu^2=1$
因此式子可以转化为

$$\begin{aligned} &\sum\limits_{i=1}^n\sum\limits_{j=1}^m lcm(i,j) (gcd(i,j) 无平方质因子)\\ &=\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)\mu^2(gcd(i,j))\\ &=\sum_{i=1}^n\sum_{j=1}^m\frac{ij}{gcd(i,j)}\mu^2(gcd(i,j))\\ &=\sum_{d=1}^{min(n,m)}\frac{1}{d}\mu^2(d)\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{m}{d}}ijd^2[gcd(i,j)==1] \qquad 枚举gcd\\ &=\sum_{d=1}^{min(n,m)}d\mu^2(d)\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{m}{d}}ij[gcd(i,j)==1]\\ &=\sum_{d=1}^{min(n,m)}d\mu^2(d)\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{m}{d}}\sum_{k|gcd(i,j)}ij\mu(k)\\ &=\sum_{d=1}^{min(n,m)}d\mu^2(d)\sum_{k=1}^{min(\frac{n}{d},\frac{m}{d})}\sum_{i=1}^{\frac{n}{kd}}\sum_{j=1}^{\frac{m}{kd}}ijk^2\mu(k)\qquad 先枚举gcd\\ &=\sum_{d=1}^{min(n,m)}d\mu^2(d)\sum_{k=1}^{min(\frac{n}{d},\frac{m}{d})}k^2\mu(k)\sum_{i=1}^{\frac{n}{kd}}i\sum_{j=1}^{\frac{m}{kd}}j\\ &=\sum_{D=1}^{min(n,m)}\sum_{d|D}d\mu^2(d)\mu(\frac{D}{d})(\frac{D}{d})^2\sum_{i=1}^{\frac{n}{D}}i\sum_{j=1}^{\frac{m}{D}}j\qquad 令D=kd\\ &=\sum_{D=1}^{min(n,m)}\sum_{d|D}\mu^2(d)\mu(\frac{D}{d})\frac{D}{d}D\sum_{i=1}^{\frac{n}{D}}i\sum_{j=1}^{\frac{m}{D}}j\\ &=\sum_{D=1}^{min(n,m)}D\sum_{d|D}\mu^2(d)\mu(\frac{D}{d})\frac{D}{d}\sum_{i=1}^{\frac{n}{D}}i\sum_{j=1}^{\frac{m}{D}}j\\ &=\sum_{D=1}^{min(n,m)}D\sum_{d|D}\mu^2(d)\mu(\frac{D}{d})\frac{D}{d}sum(1,\frac{n}{D})sum(1,\frac{m}{D})\\ \end{aligned}$$

令 $f(x)=\sum\limits_{d|x}\mu^2(d)\mu(\frac{x}{d})\frac{x}{d}$

$g(x)=\mu^2(x)$

$h(x)=\mu(x)\cdot x$

因此 $f(x)=g(x)*h(x)$ (卷积)

因为 $g(x)$ 与 $h(x)$ 为积性函数，因此 $f(x)$ 也为积性函数，可以线筛（虽然…我是看不出来的 再见）

在写程序的时候呢，其实可以线筛 $D\cdot f(x)$ ,会好做一些

代码

#include <cstdio>
#include <algorithm>
using namespace std;
#define ll long long
#define N 4000010
#define U unsigned int
int tot=0,T,notprime[N],prime[N];
U f[N];
inline void sieve(){
    f[1]=1;
    for(int i=2;i<=4000000;i++){
        if(!notprime[i]) prime[++tot]=i,f[i]=1-i;
        for(int j=1;j<=tot && prime[j]*i<=4000000;j++){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0){
                if(i%(prime[j]*prime[j])==0) f[i*prime[j]]=0;
                else f[i*prime[j]]=-f[i/prime[j]]*prime[j];
                break;
            }f[prime[j]*i]=f[i]*(1-prime[j]);
        }
    }
    for(int i=2;i<=4000000;i++) f[i]=f[i]*i+f[i-1];
}
inline U sum(int x){
    return (U)x*(x+1)/2;
}
int main(){
    scanf("%d",&T);sieve();
    while(T--){
        int n,m;scanf("%d%d",&n,&m);
        if(n>m) swap(n,m);
        U ans=0;
        for(int i=1,last;i<=n;i=last+1){
            last=min(n/(n/i),m/(m/i));
            ans+=(f[last]-f[i-1])*sum(n/i)*sum(m/i);
        }
        printf("%u\n",ans&0x3fffffff);
    }
    return 0;
}