莫比乌斯反演&整数分块

前言

参考资料

1. peng-ym-莫比乌斯反演

莫比乌斯反演

预备知识：整数分块（感觉就是一个思想）

1. 参考博客：peng-ym-整除分块
   

总结

1. 一般形式：${\sum }_{i=1}^n[n/i]$
2. 过程：r=n/(n/l)，然后不断枚举l，求得r，[l,r]部分的值是一样的，不必遍历，直接得到长度*数。
3. 复杂度：由$O(n)$变为$O(\sqrt{n})$

莫比乌斯函数$\mu$（数论函数也可以视作一个数列）

总结

1. 定义：很重要
2. 性质1：很重要
3. 性质2还需要回顾以下欧拉函数才行。

莫比乌斯反演

题目训练

题目1（oi-wiki推荐）：求满足gcd(x,y)=k且x,y分别在区间[a,b],[c,d]内的个数（莫比乌斯反演&容斥定理&整数分块）

1. 题目：P2522 [HAOI2011]Problem b
   
   

2. 题解：

   1. 首先由容斥定理把题目简化为求四个${\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)=1]$的形式。
   2. 其中$[gcd(i,j)=1]$表示满足内部条件时返回true，否则返回false。
   3. 由莫比乌斯函数性质1：${\sum }_{d|n}\mu (d)=[n=1]$
   4. 带入得到：${\sum }_{i=1}^n{\sum }_{j=1}^m{\sum }_{d|gcd(i,j)}\mu (d)$（即条件转化为具体求和式子）
   5. 变形得到：${\sum }_{d=1}^{min(n,m)}{\sum }_{i=1}^n[d|i]{\sum }_{j=1}^m[d|j]$。
   6. 最后用整数分块解决以下就ok了（不过这里怎么就敢说复杂度为$O(\sqrt{max(n,m})$呢？还是雀食不是根号，但是也大不了多少）。

3. 代码：

/*
4. 莫比乌斯反演
*/
#include <bits/stdc++.h>
using namespace std;
const int N = 5e4 + 5;
int a, b, c, d, k;
int mu[N], p[N], vis[N];

void Moblus(int n) {
    memset(vis, 0, sizeof(vis));
    mu[1] = 1;
    int tot = 0;
    for (int i = 2; i <= n; i++) {
        if (!vis[i]) p[tot++] = i, mu[i] = -1;
        for (int j = 0; j < tot; j++) {
            if (i * p[j] > n) break;
            vis[i * p[j]] = 1;
            if (i % p[j] == 0) {
                mu[i * p[j]] = 0;
                break;
            } else
                mu[i * p[j]] = -mu[i];
        }
    }
    for (int i = 1; i <= n; i++) mu[i] += mu[i - 1];
}
int solve(int n, int m) {
    int res = 0;
    for (int i = 1, j; i <= min(n, m); i = j + 1) {
        j = min(
            n / (n / i),
            m / (m / i));  // emmm，不敢多问，这你敢说是$O(\sqrt{n})$复杂度？
            // 得到一个正确的限制条件，然后模式化（无脑）枚举就好emm
        res += (mu[j] - mu[i - 1]) * (n / i) * (m / i);
    }
    return res;
}
//容斥定理
int f(int x, int k) {
    if (x % k == 0) return x / k - 1;
    return x / k;
}
signed main() {
    Moblus(N - 1);
    int T;
    cin >> T;
    while (T--) {
        scanf("%d%d%d%d%d", &a, &b, &c, &d, &k);
        int ans = solve(b / k, d / k) - solve(b / k, f(c, k)) -
                  solve(f(a, k), d / k) + solve(f(a, k), f(c, k));
        printf("%d\n", ans);
    }
    return 0;
}

题目2（cf2200）：