[NOI2010]能量采集[容斥,莫比乌斯反演]

最新推荐文章于 2019-11-27 21:38:08 发布

Hardict

最新推荐文章于 2019-11-27 21:38:08 发布

阅读量194

点赞数

分类专栏：数论-各种数论函数文章标签： gcd 容斥莫比乌斯反演

本文链接：https://blog.csdn.net/Hardict/article/details/82831486

版权

数论-各种数论函数专栏收录该内容

9 篇文章 0 订阅

订阅专栏

题意

计算 $\sum (2*i-1)f[i]$ , f[i] 表示在 $1 \leq x \leq N, 1\leq y \leq M$ 范围内 gcd(x,y)=i 的个数

数据范围: $N,M\leq 10^5$

题解

这里介绍容斥和莫比乌斯反演两种方法(如果想系统了解这种题目的各种大体思路可以看两两gcd求和的4种方法)
下面先设 F[x],f[x] 分别为 gcd(i,j) $(1 \leq i \leq N,1 \leq j \leq M)$ 为倍数和等于的个数,有 $F[n]=\sum_{n|d}f[d]=\lfloor \frac{N}{n} \rfloor\lfloor \frac{M}{n} \rfloor$

容斥(筛法)

容斥是利用 $F[n]=\sum_{n|d}f[d]=\lfloor \frac{N}{n} \rfloor\lfloor \frac{M}{n} \rfloor$ 先计算到每个 F[i] ,再通过从后向前处理得到每个 f[i] 最后针对问题选取需要的求解

  for(i = 1; i <= top; i++) F[i] = (LL)(N / i) * (M / i);
  for(i = top; i >= 1; i--){
    f[i] = F[i];
    for(j = i * 2; j <= top; j += i)
      f[i] -= f[j];
  }

莫比乌斯反演

由莫比乌斯反演知: $f[n]=\sum_{n|d}\mu (\frac{d}{n})F[d]$

问题的 $ans=\sum (2*i-1) \sum_{i|d}\mu (\frac{d}{i})[\frac{N}{d}][\frac{M}{d}] =\sum (2*i-1) \sum_{j=1}^{min([\frac{N}{i}],[\frac{M}{i}])}\mu (j)[\frac{N}{ij}][\frac{M}{ij}]$

这里如果是单次询问可以不进一步化解直接求和,时间复杂度为 O(N)

  for(i = 1; i <= top; i++) pre[i] = pre[i - 1] + mu[i];
  for(i = 1; i <= top; i++) sum[i] = sum[i - 1] + (2 * i - 1);
  int l, r; LL ans = 0;
  for(l = 1; l <= top; l = r + 1){
    r = min(N / (N / l), M / (M / l));
    ans += (LL)solve(N / l, M / l) * (sum[r] - sum[l - 1]);
  }

若进一步化解

$\\ ans=\sum_{i=1}^{min(N,M)}\sum_{j=1}^{min([\frac{N}{i}],[\frac{M}{i}])} (2i-1)\mu (j) [\frac{N}{ij}][\frac{M}{ij}]\\ =\sum_{k=1}^{min(N,M)} [\frac{N}{k}][\frac{M}{k}] \sum_{i|k} (2i-1)\mu (\frac{k}{i})$

这种化解一般都额外需要一个 O(NlogN) 的前缀和(不化解为 O(N) ),而每次计算只需要 $O(\sqrt N)$ ,一般针对多次询问

  for(i = 1; i <= top; i++)
    for(j = i; j <= top; j += i)
      sum[j] += mu[j / i] * (2 * i - 1);
  for(i = 1; i <= top; i++) sum[i] += sum[i - 1];
  int l, r; LL ans = 0;
  for(l = 1; l <= top; l = r + 1){
    r = min(N / (N / l), M / (M / l));
    ans += (LL)(N / l) * (M / l) * (sum[r] - sum[l - 1]);
  }

Hardict

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
[NOI2010]能量采集[容斥,莫比乌斯反演]

题意计算,表示在范围内的个数数据范围:题解这里介绍容斥和莫比乌斯反演两种方法(如果想系统了解这种题目的各种大体思路可以看两两gcd求和的4种方法)下面先设分别为为倍数和等于的个数,有容斥(筛法)容斥是利用先计算到每个,再通过从后向前处理得到每个最后针对问题选取需要的求解 for(i = 1; i &lt;= top; i++) F[i] = (L...
复制链接

扫一扫