莫比乌斯反演入门(三)--洛谷P3327

题目大意：

给出$n,m$,求${\sum }_{i=1}^n{\sum }_{j=1}^{m}d(i\ast j)$,$d(i\ast j)代表i\ast j的约数个数.$

$T,n,m\le 5e4$

题目思路:

首先得知道一个关于约数个数和的特殊性质:
$d(i,j)={\sum }_{x|i}^{}{\sum }_{y|j}^{}[gcd(x,y)==1]$

$i$的约数A * $j$的约数B $\in$ $i\ast j$的约数C。
但是直接这样枚举问题是会出现重复的情况.

规定一下:当$C$的素因子$p$的指数 小于等于 $i$的素因子$p$的指数时，从$i$中拿，否则就取反(全拿走，然后在$j$中去掉多余的$p$).这样就保证了一个方案中用C拿出来的$i,j$的因子A，B满足$gcd(A,B)=1$.同时，任意一对合法的$A\ast B$得到的结果$C$是唯一的。即一一映射。

---

那么开始推公式了~

$原式={\sum }_{i=1}^n{\sum }_{j=1}^m{\sum }_{x|i}^n{\sum }_{y|j}^m[gcd(x,y)==1]$

经典套路，看不得下面带整除符号，直接换枚举顺序，先枚举约数x,y(由于我们枚举了$x,y,i,j$,所以上式子的第3/4个求和式就变成了判断式(如下所示))

$={\sum }_{x=1}^n{\sum }_{y=1}^m[gcd(x,y)==1]{\sum }_{i=1}^n{\sum }_{j=1}^m[x|i\wedge y|j]$

$={\sum }_{x=1}^n{\sum }_{y=1}^m[gcd(x,y)==1]\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor$

$={\sum }_{x=1}^n{\sum }_{y=1}^m{\sum }_{d|gcd(i,j)}^{}\mu (d)\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor$

$={\sum }_{d=1}^n\mu (d){\sum }_{d|x}^n{\sum }_{d|y}^m\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor$

$={\sum }_{d=1}^n\mu (d){\sum }_i^{\lfloor \frac{n}{d}\rfloor }{\sum }_j^{\lfloor \frac{m}{d}\rfloor }\lfloor \frac{n}{id}\rfloor \lfloor \frac{m}{jd}\rfloor$

$={\sum }_{d=1}^n\mu (d){\sum }_i^{\lfloor \frac{n}{d}\rfloor }\lfloor \frac{n}{id}\rfloor {\sum }_j^{\lfloor \frac{m}{d}\rfloor }\lfloor \frac{m}{jd}\rfloor$

预处理莫比乌斯函数和$[1,n]$的整除分块结果。然后对于每个询问跑一边整除分块即可。

时间复杂度:$O(T\sqrt{n})$

AC代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<int,int>
#define pb push_back
#define mp make_pair
#define vi vector<int>
#define vl vector<ll>
#define fi first
#define se second
const int maxn = 5e4 + 5;
const int mod = 1e9 + 7;
int a[maxn];
int bk[maxn] , p[maxn] , cnt , u[maxn];
void init (){
    u[1] = 1;
    for (int i = 2 ; i < maxn ; i++){
        if (!bk[i]){
            p[++cnt] = i;
            u[i] = -1;
        }
        for (int j = 1 ; j <= cnt && i * p[j] < maxn ; j++){
            bk[i * p[j]] = 1;
            if (i % p[j] == 0){
                u[i * p[j]] = 0;
                break;
            }
            u[i * p[j]] = -u[i];
        }
    }
    for (int i = 1 ; i < maxn ; i++){
        u[i] += u[i - 1];
    }
}
ll ff[maxn];
ll f (ll n){
    ll ans = 0;
    for (ll l = 1 , r; l <= n ; l = r + 1){
        r = n / (n / l);
        ans += (r - l + 1) * (n / l);
    }
    return ans;
}
ll solve (ll n , ll m){
    ll ans = 0;
    for (ll l = 1 , r; l <= n ; l = r + 1){
        r = min(n / (n / l) , m / (m / l));
        ll a = u[r] - u[l - 1];
        ans += a * ff[n / l] * ff[m / l];
    }
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);
    init();
    for (int i = 1 ; i < maxn ; i++){
        ff[i] = f(i);
    }
    int t; cin >> t;
    while (t--){
        int n , m;
        cin >> n >> m;
        if (n > m) swap(n , m);
        cout << solve(n , m) << endl;
    }
    return 0;
}