[SDOI2015]约数个数和

一、题目

点此看题

二、解法

首先证明一个结论：
$$d(i,j)=\sum _{x|i}\sum _{y|j}[\gcd (x,y)=1]$$证明考虑单个质因数$p$，设$i$中有$p^a$，$j$中有$p^b$，等式右边的贡献是$a+b+1$，等式左边考虑$x$中没有$p$，贡献为$b$；$y$中没有$p$，贡献为$a$；$x,y$都没有$p$，贡献为$1$，总贡献也是$a+b+1$，等式右边的每个质因数互相独立，可以理解为枚举了所有情况相当于所有贡献的乘积，证毕。

---

现在用上面的结论来推式子：
$$\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}[\gcd (x,y)=1]$$我们枚举$x,y$，这样$i,j$的取值我们就知道了：
$$\sum _{x=1}^n\sum _{y=1}^m\frac{n}{x}\times \frac{m}{y}[\gcd (x,y)=1]$$$\gcd$是很容易反演的，搞一波莫比乌斯反演：
$$f(i)=\sum _{x=1}^n\sum _{y=1}^m\frac{n}{x}\times \frac{m}{y}[\gcd (x,y)=i]$$$$F(i)=\sum _{x=1}^n\sum _{y=1}^m\frac{n}{x}\times \frac{m}{y}[i|\gcd (x,y)]=\sum _{i|d}f(d)$$$$f(i)=\sum _{i|d}\mu (\frac{d}{i})F(d)$$我们要求的是$f(1)$：
$$f(1)=\sum _{i=1}^n\mu (i)F(i)$$问题变成了如何求$F(i)$，继续推式子，考虑消去$\gcd$的影响：
$$F(i)=\sum _{x=1}^{n/i}\sum _{y=1}^{m/i}\frac{n}{xi}\frac{m}{yi}$$只要我们知道了$s(x)={\sum }_{i=1}^x\frac{x}{i}$，那么$F(i)$也就迎刃而解了，$s(i)$是可以通过分块$O(n\sqrt{n})$预处理的，至于询问也可以用分块解决，时间复杂度$O(T\sqrt{n})$，贴个代码$qwq$

#include <cstdio>
#include <iostream>
using namespace std;
#define int long long
const int M = 50005;
int read()
{
    int x=0,flag=1;
    char c;
    while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
    while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
    return x*flag;
}
int T,n,m,cnt,p[M],vis[M],mu[M],s[M];
void init(int n)
{
    mu[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!vis[i])
        {
            p[++cnt]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=cnt && i*p[j]<=n;j++)
        {
            vis[i*p[j]]=1;
            if(i%p[j]==0) break;
            mu[i*p[j]]=-mu[i];
        }
    }
    for(int i=1;i<=n;i++)
    {
        mu[i]+=mu[i-1];
        for(int l=1,r;l<=i;l=r+1)
        {
            r=i/(i/l);
            s[i]+=(r-l+1)*(i/l);
        }
    }
}
signed main()
{
    T=read();
    init(5e4);
    while(T--)
    {
        n=read();m=read();
        if(n>m) swap(n,m);
        int ans=0;
        for(int l=1,r;l<=n;l=r+1)
        {
            r=min(n/(n/l),m/(m/l));
            ans+=(mu[r]-mu[l-1])*s[n/l]*s[m/l];
        }
        printf("%lld\n",ans);
    }
}