BZOJ 3994 Sdoi2015 约数个数和 莫比乌斯反演

最新推荐文章于 2019-11-06 20:09:08 发布
最新推荐文章于 2019-11-06 20:09:08 发布
阅读量4.8k 收藏 1
点赞数 1
分类专栏： BZOJ 莫比乌斯反演 文章标签： BZOJ BZOJ3994 莫比乌斯反演 数论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/popoqqq/article/details/45078079
版权
该博客探讨了BZOJ 3994题目的解决方案，涉及数论中的一个重要结论：求所有整数对(i, j)的约数个数和。作者通过莫比乌斯反演证明了特定等式的正确性，并详细解释了每个质数对答案的贡献，展示了如何利用数论技巧解决这个问题。" 119281173,7800548,使用Plotly进行双坐标轴图表可视化,"['数据可视化', 'Python编程', '数据科学', '机器学习', '深度学习']
摘要由CSDN通过智能技术生成

题目大意：求 ∑ni=1∑mj=1d(ij)
首先我们有一个很神的结论：
∑ni=1∑mj=1d(ij)=∑ni=1∑

最低0.47元/天 解锁文章
PoPoQQQ
关注
专栏目录
BZOJ 3930 CQOI2015 选数 莫比乌斯反演
世界
04-07 2万+
题目见 http://pan.baidu.com/s/1o6zajc2 此外不知道H-L #include #include #include #include #include #define M 10001000 #define INF 0x3f3f3f3f #define MOD 1000000007 using namespace std; int mu[M],p
BZOJ 3994 约束个数和（约数函数性质+莫比乌斯反演+筛积性函数）*
满月照亮的路
09-11 345
题目链接：https://cn.vjudge.net/problem/HYSBZ-3994 #include<bits/stdc++.h> #pragma comment(linker,"/STACK:1024000000,1024000000") using namespace std; #define debug puts("YES"); #define rep(x,y,z)...
bzoj3994【SDOI2005】约数个数和
AaronPolaris
04-05 1859
莫比乌斯反演(依旧没题解...hhh)
【SDOI2015】约数个数和
banchangyue3928的博客
11-22 141
题面 求\(\sum_{i=1}^n\sum_{j=1}^md(ij)\) \(\leq 50000\)组数据，\(1\leq n,m\leq 50000\)。 题目分析 首先，你需要知道一个结论： \[ d(ij)=\sum\limits_{x|i}\sum\limits_{y|j}[gcd(x,y)==1] \] 你可以认为\(x,y\)表示你选择的因数为\(\frac i x \c...
[BZOJ3994][SDOI2015]约数个数和（莫比乌斯反演）
zyf2000
05-02 1223
一个人走路，才是你和风景之间的单独私会。
Bzoj2694/Bzoj4659:莫比乌斯反演
weixin_30597269的博客
03-07 132
Bzoj2694/Bzoj4659:莫比乌斯反演先上题面:首先看到这数据范围显然是反演了，然而第三个限制条件十分不可做。于是我们暂且无视他，大不了补集转化算完再减是吧。于是我们有:这里我们定义:于是这个东西我们可以nlogn筛的说。也就是说，我们求出f的前缀和后，就可以O(sqrt(n)+sqrt(m))分块计算了。然而需要减去的东西怎么办呢？反演题最难的不是推公式，而是你推出了公式却不知道是否可...
莫比乌斯反演入门
mumei的博客
11-06 987
这个算是竞赛中用到的组合数学里最多的知识点，也是比较难理解的，在被将近一个星期的折磨后总算是理解了一些。 其实莫比乌斯反演就只有两个重要的公式，而且用到的最多的也只有一个（其实都差不多），但是需要一些预备知识。 莫比乌斯函数，质数线性筛，整除分块，积性函数以及一些其他的数论知识，下面会将其中一部分重要的算法， 目录 1.莫比乌斯函数 2.莫比乌斯函数的线性筛 3.整除分块 4....
BZOJ2693(BZOJ2154)——莫比乌斯反演经典例题
mumei的博客
11-03 455
题目链接：https://www.lydsy.com/JudgeOnline/problem.php?id=2154 题意理解： 给你n和m，求所有的lcm(i,j)之和，1<=i<=n,1<=j<=m; 很经典的莫比乌斯反演例题，从这里能真正发现它的牛皮之处。 学习这个算法需要有一点预备的知识，数论分块和唯一分解定理。 关于数论分块这里有一篇很不错的文章：htt...
莫比乌斯反演_王天懿.ppt
09-14
直接解决这个问题可能非常复杂，但是通过计算所有数对的约数和（即 \( F(n) \)），然后利用莫比乌斯反演公式，我们可以高效地求解出原问题的答案。 莫比乌斯反演的证明通常依赖于莫比乌斯函数的性质和积性函数的...
【SDOI 2015】约数个数和
forever_dreams的博客
02-17 306
算法标签：莫比乌斯反演+整除分块
「SDOI2015」约数个数和
tata_d2的博客
03-21 372
「SDOI2015」约数个数和 设d(x)d(x)d(x)为xxx的约数个数，求 ∑i=1n∑j=1md(ij)∑i=1n∑j=1md(ij)\sum_{i=1}^n \sum_{j=1}^m d(ij) n,m&lt;=50000n,m&lt;=50000n, m &lt;=50000&lt;=50000 ——————- 设ijijij唯一分解后为 ∏kqqik+qjkk∏kqkq...
bzoj 3994
clover_hxy的博客
02-16 961
设d(x)为x的约数个数，给定N,M，求 ∑i=1N∑j=1Md(ij) 。 Input 输入文件包含多组测试数据。 第一行，一个整数T，表示测试数据的组数。 接下来的T行，每行两个整数N,M。 Output T行，每行一个整数，表示你所求的答案。 Sample Input 2 7 4 5 6 Sample
BZOJ 3994
zzj的专栏
09-26 1298
BZOJ 3994Description：      ~~~~~~令 D(x)D(x) 表示 xx 的约数个数。求 ∑Ni=1∑Mj=1D(i∗j)\sum_{i=1}^{N}\sum_{j=1}^{M}D(i*j)。       ~~~~~~       ~~~~~~       ~~~~~~Solution：      ~~~~~~由对称性，不妨设 $N<=M$      ~~~~~~首先
bzoj3994/洛谷P3327 莫比乌斯反演
litble的成(tui)长(fei)史
06-24 1660
题目分析话说这是让我来学莫比乌斯反演的入坑题呢，然而学了莫比乌斯反演还是不会做=_=，连题解都看不懂QAQ 写的又臭又长是因为本蒟蒻每一步都理解了半天QAQ 数恐症患者慎入！第一步：预备的结论首先，我们要证明一条结论：d(nm)=∑i|n∑j|m[gcd(i,j)==1]d(nm)=\sum_{i|n}\sum_{j|m}[gcd(i,j)==1] 怎么证明呢？对于nmnm的每一个质因数pi
【SDOI 2015】【BZOJ 3994】约数个数和
Lcomyn的专栏
10-22 835
送个结论d(n∗m)=∑d|n∑k|ne(gcd(d,k))d(n*m)=\sum_{d|n}\sum_{k|n} e(gcd(d,k)) 记g(n)=∑ni=1⌊ni⌋g(n)=\sum_{i=1}^{n} \lfloor{n \over i} \rfloor 然后我们开始推 ∑ni=1∑mj=1d(i∗j)=∑ni=1∑mj=1⌊ni⌋⌊mj⌋e(gcd(i,j))=∑ni=1∑mj=1⌊
[SDOI2015] 约数个数和(莫比乌斯反演)
Tan_JX
08-28 188
[SDOI2015] 约数个数和 设d(x)为x的约数个数，给定N、M，求 ∑i=1N∑j=1Md(ij)\displaystyle\sum_{i=1}^N\displaystyle\sum_{j=1}^Md(ij)i=1∑N​j=1∑M​d(ij) Input 输入文件包含多组测试数据。 第一行，一个整数T，表示测试数据的组数。 接下来的T行，每行两个整数N、M。 Output T行，每行一个...
luogu P3327 [SDOI2015]约数个数和
zsyz_ZZY的博客
12-13 305
背景： 以下图片均来自我的PDFPDFPDF文件，谢绝转载。 题目传送门： https://www.luogu.org/problemnew/show/P3327 题意： 思路： 代码： #include&lt;cstdio&gt; #include&lt;cstring&gt; #include&lt;algorithm&gt; #define LL long long #...
【bzoj3994】[SDOI2015]约数个数和 莫比乌斯反演
DQS的树洞
02-23 1177
Description设d(x)为x的约数个数，给定N、M，求 Input输入文件包含多组测试数据。第一行，一个整数T，表示测试数据的组数。 接下来的T行，每行两个整数N、M。OutputT行，每行一个整数，表示你所求的答案。Sample Input27 45 6 Sample Output110121 HINT1<=N, M<=500001<=T<=50000SourceRound 1 感谢y
【bzoj3994】【SDOI2015】约数个数和【数论】【反演】
ZMOIYNLP的专栏
04-21 5756
虽然题目上写了反演但是我不知道什么是反演……如果你把Sigma调换位置叫做反演的话。 这道题题面非常简单： 设d(x)d(x)为xx的约数个数，给定N、MN、M， 求∑i=1N∑j=1Md(ij)\sum_{i=1}^N\sum_{j=1}^Md(ij) 当时我too naive，看到这玩意就默默地打50分暴力去了。。。 今天江苏神犇们做了这道题，我顺便听明白了~~ 首先它不知用什么精妙
BZOJ2627 伯努利数
最新发布
05-26
这是一道经典的组合数学题目，需要一些基本的组合数学知识。 首先，我们先来了解一下伯努利数。伯努利数是组合数学中的一类数列，它们的递推式如下： $$B_0 = 1, B_n = -\frac{1}{n+1}\sum_{i=0}^{n-1}{{n+1}\choose{i}}B_i$$ 其中 $n\geq 1$，${n\choose k}$ 表示组合数，$B_n$ 表示第 $n$ 个伯努利数。 接下来，我们考虑如何计算 $\sum_{i=0}^{n}{n+1\choose{i}}B_i$。根据二项式定理，我们有： $$(1+x)^{n+1} = \sum_{i=0}^{n+1}{n+1\choose{i}}x^i$$ 对上式两边求导可以得到： $$(n+1)(1+x)^n = \sum_{i=1}^{n+1}{n+1\choose{i}}ix^{i-1}$$ 将 $x=1$ 带入上式，得到： $$(n+1)\cdot 2^n = \sum_{i=1}^{n+1}{n+1\choose{i}}i$$ 注意到 $B_0=1$，我们可以对伯努利数的递推式进行变形： \begin{aligned} B_n &= -\frac{1}{n+1}\sum_{i=0}^{n-1}{{n+1}\choose{i}}B_i \\ &= -\frac{1}{n+1}\left({{n+1}\choose{0}}B_0 + \sum_{i=1}^{n}{{n+1}\choose{i}}B_i\right) \\ &= -\frac{1}{n+1}\left({{n+1}\choose{0}}B_0 + \sum_{i=1}^{n}{n+1\choose{i}}B_i\right) + \frac{1}{n+1}{{n+1}\choose{0}}B_0 \\ &= -\frac{1}{n+1}\sum_{i=0}^{n}{n+1\choose{i}}B_i + \frac{1}{n+1} \end{aligned} 因此，我们有： $$\sum_{i=0}^{n}{n+1\choose{i}}B_i = (n+1)(B_{n+1}-1)$$ 这个式子可以通过数学归纳法进行证明。 现在，我们已经得到了 $\sum_{i=0}^{n}{n+1\choose{i}}B_i$ 的计算公式，只需要预处理出前 $n$ 个伯努利数就可以在 $O(n)$ 的时间复杂度内解决这个问题了。
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值