
莫比乌斯反演
three trees
三木无林,却以成森
展开
-
洛谷P3312 [SDOI2014]数表 积性函数+反演+树状数组
分析我们设sd(i)sd(i)sd(i)为iii的所有约数之和如果不考虑a对答案的影响:ans=∑i=1n∑j=1msd(gcd(i,j))ans=\sum_{i=1}^{n}\sum_{j=1}^{m}sd(gcd(i,j))ans=i=1∑nj=1∑msd(gcd(i,j))正常的思路,我们枚举gcd(i,j)=dgcd(i,j)=dgcd(i,j)=dans=∑d=1min(n...原创 2019-10-31 18:56:35 · 123 阅读 · 0 评论 -
hdu4675 GCD of Sequence 组合数学+反演
分析我们设f(d)f(d)f(d)为gcd=d\gcd=dgcd=d的数量,F(d)F(d)F(d)为gcd\gcdgcd为ddd或者ddd的倍数的数量我们有反演(假设n=3):f(d)=∑d∣tμ(td)F(t)f(d)=\sum_{d \mid t}\mu(\frac{t}{d})F(t)f(d)=∑d∣tμ(dt)F(t)所以这道题难在F(t)F(t)F(t)怎么求:我们设n...原创 2019-10-17 13:01:02 · 161 阅读 · 0 评论 -
hdu 4746 Mophues 质因子数量小于p
推式子设:num[x]num[x]num[x]表示xxx唯一分解后所有质因子指数之和ans=∑i=1n∑j=1m[num[gcd(i,j)]<=p]\begin{aligned}ans=\sum_{i=1}^{n}\sum_{j=1}^{m}[num[\gcd(i,j)]<=p]\end{aligned}ans=i=1∑nj=1∑m[num[gcd(i,j)]<...原创 2019-10-15 21:58:37 · 110 阅读 · 0 评论 -
洛谷P3327-约数个数和-莫比乌斯反演+约数个数函数前缀和
推式子ans=∑i=1N∑j=1Md(ij)\begin{aligned}ans=\sum_{i=1}^{N}\sum_{j=1}^{M}d(ij)\end{aligned}ans=i=1∑Nj=1∑Md(ij)现在有个式子:d(ij)=∑x∣i∑y∣j[gcd(x,y)==1]\begin{aligned}d(ij)=\sum_{x \mid i}\sum_{y \...原创 2019-10-09 21:15:09 · 305 阅读 · 0 评论 -
洛谷-P3172 选数 杜教筛+莫比乌斯反演
题意从[L,H][L,H][L,H]选n个数,且这个n个数的gcd\gcdgcd为k的方案数推式子假设n=3ans=∑x=LH∑y=LH∑z=LH[gcd(x,y,z)==k]\begin{aligned}ans=\sum_{x=L}^{H}\sum_{y=L}^{H}\sum_{z=L}^{H}[\gcd(x,y,z)==k]\end{aligned}ans=x=L∑Hy...原创 2019-10-08 20:57:31 · 115 阅读 · 0 评论 -
gym-101982B 2018-2019 ACM-ICPC Pacific Northwest 反演
推式子ans=∑i=ab∑j=cd[gcd(i,j)==1]\begin{aligned}ans=\sum_{i=a}^{b}\sum_{j=c}^{d}[\gcd(i,j)==1]\end{aligned}ans=i=a∑bj=c∑d[gcd(i,j)==1]我们给(a,b)(a,b)(a,b)拆成(1,b)−(1,a−1)(1,b)-(1,a-1)(1,b)−(1,a−1...原创 2019-10-06 22:03:11 · 145 阅读 · 0 评论 -
洛谷--P3768 简单的数学题--杜教筛
先贴个代码#include <bits/stdc++.h>using namespace std;typedef long long ll;typedef unsigned long long ull;template <typename T>void out(T x) { cout << x << endl; }const int N...原创 2019-09-28 21:55:49 · 114 阅读 · 0 评论 -
莫比乌斯学习小结
前置知识数论分块,莫比乌斯函数,和式的相关推导(积累)数论分块这个我写过一篇博客,下面是链接https://blog.youkuaiyun.com/qq_43101466/article/details/100999784莫比乌斯函数μ(n)={1n为1(−1)kk为n的质因子个数0n存在非平方因子{\displaystyle \mu (n)={\begin{cases}1 \quad n为...原创 2019-09-27 19:42:41 · 165 阅读 · 0 评论 -
洛谷-P1829 Crash的数字表格--积性函数
推式子ans=∑i=1N∑j=1Mlcm(i,j)=∑i=1N∑j=1Mi∗jgcd(i,j)\begin{aligned}ans = \sum_{i=1}^{N}\sum_{j=1}^{M}lcm(i,j) =\sum_{i=1}^{N}\sum_{j=1}^{M}\frac{i*j}{\gcd(i,j)}\end{aligned}ans=i=1∑Nj=1∑Mlcm(i,j)...原创 2019-09-25 21:00:44 · 227 阅读 · 1 评论 -
hdu 1695 GCD
开始推式子ans=∑x=ab∑y=cd[gcd(x,y)=k]\begin{aligned}ans = \sum_{x=a}^{b}\sum_{y=c}^{d}[\gcd(x,y)=k]\end{aligned}ans=x=a∑by=c∑d[gcd(x,y)=k]还是老套路,我们令f(k)=∑x=ab∑y=cd[gcd(x,y)=k]f(k)=\sum_{x=a}^{b}...原创 2019-09-24 21:58:27 · 116 阅读 · 0 评论 -
hdu 6428 Calculate(狄利克雷卷积 + 莫比乌斯反演)
开始推式子ans=∑i=1A∑j=1B∑k=1Cϕ(gcd(i,j2,k3))\begin{aligned}ans = \sum_{i = 1}^{A}\sum_{j = 1}^{B}\sum_{k = 1}^{C}\phi(\gcd(i, j^{2}, k^{3}))\end{aligned} ans=i=1∑Aj=1∑Bk=1∑Cϕ(gcd(i,j2,k3))我们令gc...原创 2019-09-24 18:22:51 · 165 阅读 · 0 评论 -
洛谷-P4450双亲数
题意求 ans=∑i=1A∑j=1B[gcd(i,j)==d]ans = \sum_{i = 1}^{A}\sum_{j = 1}^{B}[\gcd(i, j) == d]ans=∑i=1A∑j=1B[gcd(i,j)==d]分析我们发现酷似洛谷P2257,题题解链接:https://blog.youkuaiyun.com/qq_43101466/article/details/101...原创 2019-09-22 21:36:20 · 180 阅读 · 0 评论 -
洛谷-P2257 YY的GCD
洛谷-P2257分析本题需要前置知识:数论分块,莫比乌斯反演数论分块,可看https://blog.youkuaiyun.com/qq_43101466/article/details/100999784先介绍一下莫比乌斯函数μ(n)={1n为1(−1)kk为n的质因子个数0n存在非平方因子{\displaystyle \mu (n)={\begin{cases}1 \quad n为1\\(...原创 2019-09-20 18:45:24 · 169 阅读 · 0 评论