
Mobius 反演与杜教筛
YangHao5
这个作者很懒,什么都没留下…
展开
-
bzoj 3930: [CQOI2015]选数
题目描述 今有 N,K,L,H∈N∗N,K,L,H\in\N^*N,K,L,H∈N∗,在 [L,H][L,H][L,H] 中选出 NNN 个数,使它们的最大公约数为 KKK。求方案数。 Solution 不妨把问题转化成:在 [⌊LK⌋,⌊HK⌋][\lfloor\frac LK\rfloor,\lfloor\frac HK\rfloor][⌊KL⌋,⌊KH⌋] 中取 NNN 个数使它们的最...原创 2019-07-07 16:42:50 · 182 阅读 · 0 评论 -
bzoj 4916: 神犇和蒟蒻
题目描述 给你正整数 nnn,求∑i=1nμ(i2)f(n)=∑i=1nφ(i2)\sum_{i=1}^{n}\mu(i^2)\\ f(n)=\sum_{i=1}^{n}\varphi(i^2)i=1∑nμ(i2)f(n)=i=1∑nφ(i2)的值。 Solution 首先,第一行是搞笑的,因为 ∀i∈[2,∞)\forall i\in[2,\infty)∀i∈[2,∞) 都有 μ(i2)=...原创 2019-07-08 10:28:27 · 217 阅读 · 0 评论 -
bzoj 3994: [SDOI2015]约数个数和
题目描述 给你 N,MN,MN,M,求∑i=1N∑j=1Md(ij)\sum_{i=1}^{N}\sum_{j=1}^{M}d(ij)i=1∑Nj=1∑Md(ij)你需要回答 TTT 个询问。1≤N,M,T≤50 0001\leq N,M,T\leq50\ 0001≤N,M,T≤50 000。 Solution 今有结论(*)d(ij)=∑x∣i∑y∣j[gcd(x,y...原创 2019-07-09 08:38:41 · 175 阅读 · 0 评论 -
bzoj 2005: [Noi2010]能量采集
题目描述 给你 n,m∈N∗n,m\in\N^*n,m∈N∗,求∑i=1n∑j=1m(gcd(i,j)−1)\sum_{i=1}^{n}\sum_{j=1}^{m}(\gcd(i,j)-1)i=1∑nj=1∑m(gcd(i,j)−1) Solution 设f=∑i=1n∑j=1mgcd(i,j)=∑i=1n∑j=1m∑d∣gcd(i,j)φ(d)=∑d=1min(n,m)φ(d)(∑...原创 2019-07-09 10:45:17 · 198 阅读 · 0 评论 -
bzoj 3512: DZY Loves Math IV
题目描述 今有 n,m∈N∗n,m\in\N^*n,m∈N∗,求∑i=1n∑j=1mφ(ij)mod  (109+7)\sum_{i=1}^{n}\sum_{j=1}^{m}\varphi(ij)\mod(10^9+7)i=1∑nj=1∑mφ(ij)mod(109+7) 的值,n≤105,m≤109n\leq10^5,m\leq10^9n...原创 2019-07-05 10:42:31 · 622 阅读 · 0 评论 -
luogu P3768 简单的数学题
题目描述 今有 n,p∈N∗n,p\in\N^*n,p∈N∗,求∑i=1n∑j=1nij⋅gcd(i,j)\sum_{i=1}^{n}{\sum_{j=1}^{n}{ij·\gcd(i,j)}}i=1∑nj=1∑nij⋅gcd(i,j) n≤1010n\leq10^{10}n≤1010。 Solution 设S(n)=∑i=1n∑j=1nij⋅gcd(i,j)=∑k=1n∑i=1n∑j=...原创 2019-07-04 15:09:17 · 208 阅读 · 0 评论 -
Mobius 反演与杜教筛
杜教筛 杜教筛 是在低于线性的时间复杂度,求一个积性函数前缀和的算法。 (找不到定义,只好自己编一个) 积性函数 积性函数 指对于所有互质的整数 aaa 和 bbb 有性质 f(ab)=f(a)f(b)f(ab)=f(a)f(b)f(ab)=f(a)f(b) 的数论函数。 特别地,若所有的整数 aaa 和 bbb 有性质 f(ab)=f(a)f(b)f(ab)=f(a)f(b)f(ab)=f(a...原创 2019-07-03 11:15:16 · 267 阅读 · 0 评论 -
bzoj 4805: 欧拉函数求和
背景 删代码就能AC另一题的感觉真好。 题目描述 今有 n∈[1,2×109]n\in[1,2\times10^9]n∈[1,2×109],求∑i=1nφ(i)\sum_{i=1}^{n}{\varphi(i)}i=1∑nφ(i) Solution原创 2019-07-06 16:46:28 · 272 阅读 · 0 评论