Luogu P3327 [SDOI2015]约数个数和___莫比乌斯函数

最新推荐文章于 2021-03-21 18:31:13 发布

disPlayLzy_

最新推荐文章于 2021-03-21 18:31:13 发布

阅读量167

点赞数

分类专栏： C++ 莫比乌斯反演

本文链接：https://blog.csdn.net/Gx_Man_VIP/article/details/98934078

版权

C++ 同时被 2 个专栏收录

368 篇文章 2 订阅

订阅专栏

莫比乌斯反演

5 篇文章 0 订阅

订阅专栏

题目大意：

在这里插入图片描述

分析：

$d(ij)=\sum_{x|i}\sum_{y|j}[gcd(x,j/y)=1]$
$ans=\sum_{i=1}^{n}\sum_{j=1}^{m}d(ij)$
$=\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{x|i}\sum_{y|j}[gcd(x,j/y)=1]$
$=\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{x|i}\sum_{y|j}\sum_{d|gcd(x,j/y)}μ(d)$
$=\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{d|gcd(i,j)}μ(d)\sum_{d|x,x|i}\sum_{d|(j/y),y|j}$
设 $x = x 1 * d$ , $i = s * d$ ,因为 $x ∣ i$ 所以 $x 1 ∣ s$ ,因此 $x 1 ∣ (i / d)$
$=\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{d|gcd(i,j)}μ(d)\sum_{x1|(i/d)}\sum_{y|(j/d)}$
设 $d (i)$ 表示 $i$ 的约数个数，
$=\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{d|gcd(i,j)}μ(d)d(i/d)d(j/d)$
$=\sum_{d|n}μ(d)\sum_{i=1}^{n/d}d(i)\sum_{j=1}^{m/d}d(j)$
$O (n l o g n)$ 处理出d()
对 $μ (d)$ 和 $d (i)$ 做前缀和，
然后每个询问整除分块即可
时间复杂度:O(nlogn+T $\sqrt n$ )

代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <queue>
#include <cstring> 
#include <algorithm>

#define rep(i, st, ed) for (int i = st; i <= ed; i++)
#define rwp(i, ed, st) for (int i = ed; i >= st; i--)

#define N 50005

using namespace std;

typedef long long ll;

int sum1[N], sum2[N], pri[N], mu[N], C[N], n, m, T, cnt;
bool vis[N];
ll ans;
 
void Pre_Work() {
    mu[1] = 1;
    rep(i, 2, N - 5) {
         if (!vis[i]) { pri[++cnt] = i; mu[i] = -1; }
         rep(j, 1, cnt) {
              if ((ll)pri[j] * i > N - 5) break;
              vis[pri[j] * i] = 1;
              if (i % pri[j] == 0) { mu[pri[j] * i] = 0; break; } else mu[pri[j] * i] = - mu[i];
         }
    }
    rep(i, 1, N - 5) sum1[i] = sum1[i - 1] + mu[i];
    rep(i, 1, N - 5)
       rep(j, 1, (N - 5) / i) C[i * j]++;
    rep(i, 1, N - 5) sum2[i] = sum2[i - 1] + C[i];
}

int main() { 
	Pre_Work();
	scanf("%d", &T);
	rep(i, 1, T)  {
		scanf("%d %d", &n, &m);
		if (n > m) swap(n, m);
		int L = 1, R = 0;
		ans = 0;
		while (L <= n) {
            R = min(n / (n / L), m / (m / L));
            ans = ans + 1ll * (sum1[R] - sum1[L - 1]) * sum2[n / L] * sum2[m / L];
            L = R + 1;
        }
        printf("%lld\n", ans);
	}
	return 0;  
}