莫比乌斯反演
ixRic
ixRic
展开
-
[51Nod 1584] 加权约数和(约数和函数性质 + 莫比乌斯反演) | 错题本
文章目录题目分析代码题目[51Nod 1584] 加权约数和分析先把 max{i,j}\max\{i, j\}max{i,j} 去掉:原式=2∑i=1N∑j=1iiσ(ij)−∑i=1Niσ(i2)原式=2\sum_{i = 1}^{N}\sum_{j = 1}^i i\sigma(ij)-\sum_{i = 1}^N i\sigma(i^2)原式=2i=1∑Nj=1∑iiσ(ij)−i=1∑Niσ(i2)推式子前先介绍约数和函数有一个重要性质(和约数个数一样可以转化为 gcd\gc原创 2020-09-17 10:50:44 · 225 阅读 · 0 评论 -
[SDOI 2014] 数表(莫比乌斯反演 + 树状数组) | 错题本
文章目录题目分析题目[SDOI 2014] 数表分析∑i=1n∑j=1m⌊alcm(i,j)⌋=∑i=1n∑j=1m⌊alcm(i,j)ij⌋\begin{aligned}&\sum_{i = 1}^n\sum_{j = 1}^m \left\lfloor\frac{a}{\text {lcm}(i, j)}\right\rfloor \\ =&\sum_{i = 1}^n\sum_{j = 1}^m \left\lfloor\frac{a\text{lcm}(i, j)}{ij}\原创 2020-09-11 10:58:10 · 143 阅读 · 0 评论 -
[国家集训队] Crash的数字表格 / JZPTAB(莫比乌斯反演) | 错题本
文章目录题目分析代码题目[国家集训队] Crash的数字表格 / JZPTAB分析不妨设 n≥mn \geq mn≥m:∑i=1n∑j=1mlcm(i,j)=∑i=1n∑j=1mijgcd(i,j)=∑d=1m∑i=1⌊nd⌋∑j=1⌊md⌋ijd[gcd(i,j)=1]=∑d=1m1d∑k=1⌊md⌋μ(k)∑i=1⌊ndk⌋∑j=1⌊mdk⌋(idk)(jdk)=∑d=1md∑k=1⌊md⌋μ(k)k2∑i=1⌊ndk⌋∑j=1⌊mdk⌋ij\begin{aligned} &\su原创 2020-09-09 11:21:46 · 182 阅读 · 0 评论 -
[RC - 02] GCD(莫比乌斯反演 + 杜教筛) | 错题本
文章目录题目分析代码题目[RC - 02] GCD分析枚举 ppp 到 ⌊Nj⌋\left\lfloor\frac{N}{j}\right\rfloor⌊jN⌋ 实际上就是在枚举 jpjpjp,于是原式变成∑i=1N∑j=1N∑p=1N∑q=1N[gcd(i,j)=1][gcd(p,q)=j]=∑i=1N∑j=1N∑p=1N∑q=1N[gcd(i,p,q)=1]=∑d=1N⌊Nd⌋3μ(d)\begin{aligned}&\sum_{i = 1}^N\sum_{j = 1}^N\s原创 2020-09-09 10:19:01 · 199 阅读 · 0 评论 -
[洛谷 P4450] 双亲数(莫比乌斯反演) | 错题本
文章目录题目分析代码题目[洛谷 P4450] 双亲数分析同 [POI2007]ZAP-Queries。代码#include <bits/stdc++.h>typedef long long LL;const int MAXN = 1000000;int Mu[MAXN + 5];bool Vis[MAXN + 5];std::vector<int> Primes;void Init(int n) { Mu[1] = 1; for (int i =原创 2020-09-08 11:35:22 · 145 阅读 · 0 评论 -
[POI 2007] ZAP-Queries(莫比乌斯反演) | 错题本
文章目录题目分析代码题目[POI2007]ZAP-Queries分析同 [HAOI 2011] Problem b。代码很迷惑的是筛莫比乌斯函数的时候要把 MAXN 开大一点点不然只有 80 pts。#include <bits/stdc++.h>typedef long long LL;const int MAXN = 60000;int Mu[MAXN + 5];bool Vis[MAXN + 5];std::vector<int> Primes;原创 2020-09-08 11:28:49 · 156 阅读 · 0 评论 -
[SDOI 2015] 约数个数和(约数个数函数性质 + 莫比乌斯反演) | 错题本
文章目录题目分析代码题目[SDOI 2015] 约数个数和分析d(xy)=∑d1∣x∑d2∣y[gcd(d1,d2)=1]d(xy)=\sum_{d_1|x}\sum_{d_2|y}[\gcd(d_1,d_2) = 1]d(xy)=d1∣x∑d2∣y∑[gcd(d1,d2)=1] 我们证明每个因子与 gcd(d1,d2)=1\gcd(d_1,d_2) = 1gcd(d1,d2)=1 的双射关系:设 x=∏i=1npiai,y=∏i=1npibi (ai+bi≠0)x原创 2020-09-08 11:08:46 · 239 阅读 · 0 评论 -
[CQOI 2015] 选数(莫比乌斯反演 + 杜教筛) | 错题本
文章目录题目分析题目[CQOI 2015] 选数分析∑a1=LH∑a2=LH⋯∑aN=LH[gcdi=1Nai=K]=∑a1=⌈LK⌉⌊HK⌋∑a2=⌈LK⌉⌊HK⌋⋯∑aN=⌈LK⌉⌊HK⌋[gcdi=1Nai=1]=∑d=1⌊HK⌋μ(d)f(d)\begin{aligned}&\sum_{a_1= L}^{H}\sum_{a_2 = L}^{H} \cdots \sum_{a_N = L}^{H}[\gcd_{i = 1}^N a_i = K] \\ =&\sum_{a_原创 2020-09-07 12:02:39 · 174 阅读 · 0 评论 -
[HAOI 2011] Problem b(莫比乌斯反演) | 错题本
文章目录题目分析代码题目[HAOI 2011] Problem b分析转换成求四块的答案,然后和 YY的GCD 一样。代码#include <bits/stdc++.h>typedef long long LL;const int MAXN = 50000;int Mu[MAXN + 5];bool Vis[MAXN + 5];std::vector<int> Primes;void Init(int n) { Mu[1] = 1; for (in原创 2020-09-05 11:34:33 · 216 阅读 · 0 评论 -
[洛谷 P2257] YY的GCD(莫比乌斯反演) | 错题本
文章目录题目分析代码题目[洛谷 P2257] YY的GCD分析不妨设 N≥MN \geq MN≥M,PPP 是素数集,枚举 d=gcd(x,y),i=xd,j=ydd = \gcd(x, y), i = \frac{x}{d}, j = \frac{y}{d}d=gcd(x,y),i=dx,j=dy,要求的即 ∑d=1M[d∈P]∑i=1⌊Nd⌋∑j=1⌊Md⌋[gcd(i,j)=1]=∑d=1M[d∈P]∑i=1⌊Nd⌋∑j=1⌊Md⌋∑k∣gcd(i,j)μ(k)=∑d=1M[d∈P]原创 2020-09-05 11:17:00 · 199 阅读 · 0 评论 -
[CodeForces 839D] Winter is here(莫比乌斯反演) | 错题本
文章目录题目分析代码题目[CodeForces 839D] Winter is here分析设 mmm 是最大的数,枚举 gcd=d\gcd = dgcd=d,子序列长度为 lll,要求的就是 ∑d=2m∑l=1nd⋅l∑1≤k1<k2<⋯<kl≤n[gcdi=1laki=d]=∑d=2m∑l=1nd⋅l∑1≤k1<k2<⋯<kl≤n(∏i=1l[d∣aki])[gcdi=1lakid=1]\begin{aligned} &\sum_{d = 2}^原创 2020-09-05 09:47:33 · 206 阅读 · 0 评论 -
[CodeForces 547C] Mike and Foam(莫比乌斯反演) | 错题本
文章目录题目分析代码题目[CodeForces 547C] Mike and Foam分析设 cic_ici 表示 iii 的个数,要求 i<ji < ji<j 比较麻烦,考虑求 t=∑i=1n∑j=1n[gcd(ai,aj)=1]t = \sum_{i = 1}^n\sum_{j = 1}^{n}[\gcd(a_i, a_j) = 1]t=i=1∑nj=1∑n[gcd(ai,aj)=1] 则答案为 t−c12\frac{t - c_1}{2}2t−c1。设 mmm原创 2020-09-04 22:14:21 · 236 阅读 · 0 评论 -
[CodeForces 439E] Devu and Birthday Celebration(莫比乌斯反演) | 错题本
文章目录题目分析代码题目[CodeForces 439E] Devu and Birthday Celebration分析莫比乌斯函数比较重要的性质:μ∗1=ε\mu * 1 = \varepsilonμ∗1=ε 即 ∑d∣nμ(d)=[n=1]\sum_{d | n}\mu(d) = [n = 1]d∣n∑μ(d)=[n=1]证明:设 n=∏i=1mpiαin = \prod_{i = 1}^{m} p_i^{\alpha_i}n=∏i=1mpiαi,n′=∏i=1mpin' = \原创 2020-09-04 12:33:29 · 164 阅读 · 0 评论