快速莫比乌斯变换(FMT)与莫比乌斯反演 例题:树上lcm

快速莫比乌斯变换

数学公式

  • 记\(S\)为全集,\(T\)为其子集

\[\begin{align} &zeta变换:F(S)=\sum_{T\subseteq S}f(T)\\ \\ &莫比乌斯反演:f(S)=\sum_{T\subseteq S}(-1)^{|S|-|T|}F(T) \end{align} \]

  • \(zeta\)变换与\(sosdp\)中的子集求和、超集求和一致,只不过这里的集合范畴不仅限于二进制集合
  • 莫比乌斯反演的作用就在于将超集和或者子集和反演为原来的函数

代码实现(二进制集合)

\(zeta\)变换(子集和)

伪代码:

r 复制代码
for 每一位 bit i:
    for 每个 mask:
        如果 mask 有第 i 位:
            F[mask] += F[mask 去掉第 i 位]

代码:

cpp 复制代码
void zeta(vector<ll> &f, int n) {
    for (int i = 0; i < n; i++) { // 枚举每一位
        for (int mask = 0; mask < (1 << n); mask++) {
            if (mask & (1 << i)) { // 如果第 i 位是 1
                f[mask] += f[mask ^ (1 << i)];
            }
        }
    }
}

\(Mobius\)反演

伪代码:

r 复制代码
for 每一位 bit i:
    for 每个 mask:
        如果 mask 有第 i 位:
            F[mask] -= F[mask 去掉第 i 位]

代码:

cpp 复制代码
void mobius(vector<ll> &f, int n) {
    for (int i = 0; i < n; i++) {
        for (int mask = 0; mask < (1 << n); mask++) {
            if (mask & (1 << i)) {
                f[mask] -= f[mask ^ (1 << i)];
            }
        }
    }
}

数论中的莫比乌斯反演

若有:

\[F(n)=\sum_{d|n}f(d) \]

其中\(d|n\)表示\(d\)为\(n\)的因数

则有:

\[f(n)=\sum_{d|n}\mu(d)·F\left( \frac{n}{d} \right) \]

其中\(\mu(d)\)为数论莫比乌斯函数,对应着莫比乌斯反演中的\((-1)^{|S|-|T|}\);\(n\)对应着全集\(S\);\(\frac{n}{d}\)对应着子集\(T\)

莫比乌斯函数的定义:

\[\mu(n)= \begin{cases} \begin{align} &1,&n=1\\ \\ &(-1)^k &n是k个不同质数的乘积\\ \\ &0 &n有平方因子 \end{align} \end{cases} \]

基于欧拉筛的\(\mu\)函数代码:

cpp 复制代码
int mu[MX], vis[MX], p[MX], cnt;
void init() {
    mu[1] = 1;
    rep(i, 2, MX - 1) {
        if (!vis[i])p[++cnt] = i, mu[i] = -1;
        for (int j = 1; i * p[j] < MX; j++) {
            vis[i * p[j]] = 1;
            if (i % p[j] == 0)break;
            mu[i * p[j]] = -mu[i];
        }
    }
}

例题:树上lcm

思路

设\(f(x)\)为路径\(lcm\)为\(x\)的简单路径数

由于涉及因数与求和,联想到莫比乌斯反演

\[设f(x)=\sum_{d|x}\mu(d)g\left( \frac{x}{d} \right) \]

由莫比乌斯反演:

\[g(n)=\sum_{d|n}f(d) \]

因此得到\(g(n)\)的定义:路径\(lcm\)为\(n\)的因数的简单路径数

因此可以先将整棵树上不是\(x\)的因数的节点删去,此时每个连通块内的任意两点构成的简单路径的\(lcm\)都是\(x\)的因数!

dfs染色统计每个连通块的大小,计算有多少条简单路径即可

代码实现

cpp 复制代码
#include<iostream>
#include<cstdio>
#include<vector>
#include<algorithm>
#include<queue>
#include<cmath>
#include<unordered_map>
#include<set>
using namespace std;
using ll = long long;
#define rep(i, a, b) for(ll i = (a); i <= (b); i ++)
#define per(i, a, b) for(ll i = (a); i >= (b); i --)
//#define see(stl) for(auto&ele:stl)cout<<ele<<" "; cout<<'\n';
#define int ll
const int N = 1e5 + 5;
const int MX = 1e7 + 5;
int n, x;
struct node {
    vector<int>e;
    int val, tag, siz;
}a[N];

ll gcd(ll a, ll b) {
    if (b == 0)return a;
    return gcd(b, a % b);
}
ll lcm(ll a, ll b) {
    return a * b / gcd(a, b);
}

void dfs(int u, int fa, int fac) {
    if (lcm(a[u].val, fac) > fac)a[u].tag = -1;
    else a[u].tag = 0;
    for (auto& son : a[u].e) {
        if (son == fa)continue;
        dfs(son, u, fac);
    }
}

void dfs2(int u, int fa) {
    a[u].tag = 1; a[u].siz = 1;
    for (auto& son : a[u].e) {
        if (son == fa || a[son].tag != 0)continue;
        dfs2(son, u);
        a[u].siz += a[son].siz;
    }
}

int mu[MX], vis[MX], p[MX], cnt;
void init() {
    mu[1] = 1;
    rep(i, 2, MX - 1) {
        if (!vis[i])p[++cnt] = i, mu[i] = -1;
        for (int j = 1; i * p[j] < MX; j++) {
            vis[i * p[j]] = 1;
            if (i % p[j] == 0)break;
            mu[i * p[j]] = -mu[i];
        }
    }
}

void eachT() {
    set<int>fac;
    rep(i, 1, n)a[i].e.clear(), a[i].tag = 0,a[i].siz=1;

    cin >> n >> x;
    rep(i, 1, n - 1) {
        int u, v; cin >> u >> v;
        a[u].e.push_back(v), a[v].e.push_back(u);
    }
    rep(i, 1, n)cin >> a[i].val;

    for (int i = 1; i * i <= x; i++) {
        if (x % i == 0)fac.insert(i), fac.insert(x / i);
    }

    ll ans = 0;
    for (auto& f : fac) {
        ll now = 0;
        dfs(1, 0, f);
        rep(i, 1, n) {
            if (a[i].tag == 0) {
                dfs2(i, 0);
                int s = a[i].siz;
                now += s * (s - 1) / 2 + s;
            }
        }
        ans += now * mu[x / f];
    }
    cout << ans << '\n';
}
signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(0), cout.tie(0);
    init();
    ll t = 1;
    cin >> t;
    while (t--) { eachT(); }
}
相关推荐
CUC-MenG1 天前
2025牛客多校第五场 K.完美旅程 J.最快覆盖问题 E.神秘异或操作 个人题解
数学·dfs·bfs·优先队列·二分·位运算·fmt·曼哈顿距离·fwt
Always_away4 天前
26考研|数学分析:重积分
笔记·学习·考研·数学
屈臣4 天前
AtCoder Beginner Contest 417 (A-E题解)
动态规划·dfs·二分
一只小蒟蒻16 天前
DFS 迷宫问题 难度:★★★★☆
算法·深度优先·dfs·最短路·迷宫问题·找过程
w909519 天前
【9】斯特林数学习笔记
数学·学习笔记
墨风如雪20 天前
8B 模型吊打 671B?数学证明界“卷王”Goedel-Prover-V2 来了!
数学·aigc
minglie123 天前
高斯代数基本定理的一种证明
数学
minglie123 天前
代数基本定理最简短的证明
数学
极大理想23 天前
可数集与不可数集
数学·几何学·拓扑学