洛谷P3327：[SDOI2015]约数个数和 （莫比乌斯反演）

题目传送门：https://www.luogu.org/problem/show?pid=3327

---

题目分析：这题我又没有自己想出来……
主要是本题要用到一个很神的结论：

$$d(ij)=\sum_{x|i}\sum_{y|j}[(x,y)=1]$$
这个是怎么推出来的呢？
我们考虑质数p对d(ij)的贡献，假设i的质因数分解中有k个p，j的质因数分解中q个p，那么d(ij)中就会有因数k+q+1，而：
$$k+q+1=\sum_{x=0}^{k}\sum_{y=0}^{q}[(p^x,p^y)=1]$$
这就相当于一个(k+1)*(q+1)的矩形，只取了左下角的”L”字型的一块。
现在我们假设 $i=p_1^{k_1}*p_2^{k_2}……p_n^{k_n},j=p_1^{q_1}*p_2^{q_2}……p_n^{q_n}$，那么我们逐个考虑每一个质数 $p_i$，得到一条式子：
$$d(ij)=\sum_{x_1=1}^{k_1}\sum_{y_1=1}^{q_1}[(p_1^{x_1},p_1^{y_1})=1]……\sum_{x_n=1}^{k_n}\sum_{y_n=1}^{q_n}[(p_n^{x_n},p_n^{y_n})=1]$$
现在我们记 $x=p_1^{x_1}p_2^{x_2}……p_n^{x_n},y=p_1^{y_1}p_2^{y_2}……p_n^{y_n}$，则变成：
$$d(ij)=\sum_{x|i}\sum_{y|j}[(x,y)=1]$$
是不是很神奇？

---

现在回到原题，经上面变形，我们要求：

$$\sum_{i=1}^n\sum_{j=1}^m\sum_{x|i}\sum_{y|j}[(x,y)=1]$$
我们将x,y的枚举拉到i,j的前面，得到：
$$\sum_{x=1}^n\sum_{y=1}^m[(x,y)=1]\lfloor\frac{n}{x}\rfloor\lfloor\frac{m}{y}\rfloor$$
我们用莫比乌斯函数将中括号内的东西替代掉，得到：
$$\sum_{x=1}^n\sum_{y=1}^m\sum_{k|x,k|y}\mu(k)\lfloor\frac{n}{x}\rfloor\lfloor\frac{m}{y}\rfloor$$
将k的枚举拉到前面：
$$\sum_{k=1}^{min(n,m)}\mu(k)\sum_{x=1}^{\lfloor\frac{n}{k}\rfloor}\lfloor\frac{n}{kx}\rfloor\sum_{y=1}^{\lfloor\frac{m}{k}\rfloor}\lfloor\frac{m}{ky}\rfloor$$
我们设 $g(n)=\sum_{i=1}^{n}\lfloor\frac{n}{i}\rfloor$，那么右边的部分变成 $g(\lfloor\frac{n}{k}\rfloor)g(\lfloor\frac{m}{k}\rfloor)$，于是我们左边用下底函数分块，右边预处理g函数即可。
我们可以用 $O(n\sqrt{n})$通过下底函数分块求出每一个g，但我们发现其实g函数就是约数个数的前缀和，直接线性筛即可。

---

CODE:

#include<iostream>
#include<string>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<stdio.h>
#include<algorithm>
using namespace std;

const int maxn=50005;
typedef long long LL;

int sum[maxn];
bool vis[maxn];
int prime[maxn];
int cur=0;

LL g[maxn];
int t,n,m;
LL ans;

void Preparation()
{
    for (int i=1; i<maxn; i++)
    {
        g[i]=0;
        int last;
        for (int j=1; j<=i; j=last+1)
        {
            last=i/(i/j);
            g[i]=g[i]+(long long)( i/j*(last-j+1) );
        }
    }
    sum[1]=1;
    for (int i=2; i<maxn; i++)
    {
        if (!vis[i]) sum[i]=-1,prime[++cur]=i;
        for (int j=1; j<=cur && i*prime[j]<maxn; j++)
        {
            int k=i*prime[j];
            vis[k]=true;
            if (i%prime[j]) sum[k]=-sum[i];
            else
            {
                sum[k]=0;
                break;
            }
        }
    }
    for (int i=2; i<maxn; i++) sum[i]+=sum[i-1];
}

int main()
{
    freopen("c.in","r",stdin);
    freopen("c.out","w",stdout);

    Preparation();
    scanf("%d",&t);
    while (t--)
    {
        scanf("%d%d",&n,&m);
        if (n>m) swap(n,m);
        ans=0;
        int last;
        for (int i=1; i<=n; i=last+1)
        {
            last=min( n/(n/i) , m/(m/i) );
            ans=ans+( (long long)(sum[last]-sum[i-1])*g[n/i]*g[m/i] );
        }
        cout<<ans<<endl;
    }

    return 0;
}