牛客挑战赛47 F 简单题（莫比乌斯反演）

最新推荐文章于 2024-06-18 23:00:58 发布

Ilnay

最新推荐文章于 2024-06-18 23:00:58 发布

阅读量206

点赞数 1

分类专栏：题解

本文链接：https://blog.csdn.net/qq_51682942/article/details/119741584

版权

题解专栏收录该内容

8 篇文章 0 订阅

订阅专栏

这篇博客介绍了牛客挑战赛47F的一道题目，主要涉及数学概念莫比乌斯反演和狄利克雷卷积。通过解析题目，作者展示了如何利用积性函数的性质和狄利克雷后缀和优化算法，将时间复杂度从O(nlogn)降低到O(nloglogn)，以解决数论问题。代码部分演示了C++实现的详细过程。

摘要由CSDN通过智能技术生成

牛客挑战赛47 F 简单题（莫比乌斯反演）

在这里插入图片描述

题解

$当gcd(i,j)>1时：\mu(ij)中包含平方因子gcd(i,j)，此时\mu(ij)=0，因此只需枚举gcd(i,j)=1的情况：$
$\sum_{i=1}^{n}\sum_{j=1}^{m}\varphi(ij)\mu(ij)[gcd(i,j)==1]$ $由积性函数的性质可得：当 g c d (i, j) = 1 时：$ $\varphi(ij)=\varphi(i)\varphi(j)，\mu(ij)=\mu(i)\mu(j)$ $\sum_{i=1}^{n}\sum_{j=1}^{m}\varphi(i)\varphi(j)\mu(i)\mu(j)[gcd(i,j)==1]$ $=\sum_{i=1}^{n}\sum_{j=1}^{m}\varphi(i)\varphi(j)\mu(i)\mu(j)\sum_{d|gcd(i,j)}\mu(d)$ $=\sum_{i=1}^{n}\varphi(i)\mu(i)\sum_{j=1}^{m}\varphi(j)\mu(j)\sum_{d|i,d|j}\mu(d)$ $=\sum_{d=1}^{min(n,m)}\mu(d)(\sum_{i=1}^{n}\varphi(i)\mu(i)[d|i])(\sum_{j=1}^{m}\varphi(j)\mu(j)[d|j])$
$令f(d)=\sum_{i=1}^{n}\varphi(i)\mu(i)[d|i]，g(d)=\sum_{i=1}^{m}\varphi(i)\mu(i)[d|i]，最后就可以化简成：$ $\sum_{d=1}^{min(n,m)}\mu(d)f(d)g(d)$
$我们首先筛出函数\varphi和\mu，此时如果暴力枚举d，然后f(d)和g(d)，时间复杂度为O(nlogn)，会T，但是我们观察到函数f和g的形式就是狄利克雷后缀和的形式，所以可以用狄利克雷后缀和预处理出来函数f和g，再来枚举d，这样可以优化到O(nloglogn)，就可以过了$

代码

#include<iostream>
using namespace std;
using uint = unsigned int;
const int N = 2e7 + 5;
int T, n, m, p[N], phi[N], mu[N], cnt;
bool vis[N]; uint f[N], g[N];
void get(int n) {
    phi[1] = 1, mu[1] = 1;
    for(int i = 2; i <= n; i++) {
        if(!vis[i]) {
            p[cnt++] = i;
            phi[i] = i - 1;
            mu[i] = -1;
        }
        for(int j = 0; p[j] <= n / i; j++) {
            vis[i * p[j]] = 1;
            if(i % p[j] == 0) {
                phi[i * p[j]] = phi[i] * p[j];
                break;
            }
            phi[i * p[j]] = phi[i] * (p[j] - 1);
            mu[i * p[j]] = -mu[i];
        }
    }
}
int main() {
    cin >> T;
    get(N - 5);
    while(T--) {
        cin >> n >> m;
        for(int i = 1; i <= n; i++) f[i] = uint(phi[i] * mu[i]);
        for(int i = 1; i <= m; i++) g[i] = uint(phi[i] * mu[i]);
        for(int i = 0; i < cnt && p[i] <= n; i++)
            for(int j = n / p[i]; j; j--) f[j] += f[j * p[i]];
        for(int i = 0; i < cnt && p[i] <= m; i++)
            for(int j = m / p[i]; j; j--) g[j] += g[j * p[i]];
        uint res = 0;
        for(int i = 1; i <= min(n, m); i++)
            res += mu[i] * f[i] * g[i];
        cout << res << endl;
    }
    return 0;
}