【题解】Crash的数字表格 BZOJ 2154 莫比乌斯反演
题目传送门 http://www.lydsy.com/JudgeOnline/problem.php?id=2154
人生中第一道自己做出来的莫比乌斯反演
人生中第一篇用LaTeX写数学公式的博客
大家别看公式多就害怕了啊,这里面的公式大多是很显然的
首先,题目要我们求
$\Large\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}lcm(i,j)$
这个$lcm$很不好办,我们想办法转化成$gcd$,然后尝试搞莫比乌斯反演的套路
那么因为
$\Large lcm(i,j)=\frac{i \cdot j}{gcd(i,j)}$
所以我们要求的就是
$\Large\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}\frac{i \cdot j}{gcd(i,j)}$
于是我们可以枚举每个最大公约数$g$,对于每个$g$,我们求
$\Large\sum\limits_{g=gcd(i,j)}\frac{i \cdot j}{g}$
最后把这些都加起来就好了
继续变形,得到
$\Large\sum\limits_{g=gcd(i,j)}\frac{i}{g} \cdot \frac{j}{g} \cdot {g}$
然后我们设函数
$\Large f(a,b,x)=\sum\limits_{x=gcd(i,j),i \le a,j \le b}\frac{i}{x} \cdot \frac{j}{x}$
于是我们要最终的答案就是
$\Large \sum\limits_{g=1}^{min(n,m)}f(n,m,g) \cdot g$
注意到
$\Large f(a,b,x)=f( \lfloor \frac{a}{x} \rfloor, \lfloor \frac{b}{x} \rfloor,1)$
这一步的原因请读者自行思考
所以最终的答案是
$\Large\sum\limits_{g=1}^{min(n,m)}f( \lfloor \frac{n}{g} \rfloor, \lfloor \frac{m}{g} \rfloor,1) \cdot g$
因为$\lfloor \frac{n}{g} \rfloor$和$\lfloor \frac{m}{g} \rfloor$都是整除,于是这一步我们就可以分块优化,不会分块优化的请先学习 http://sxysxy.org/blogs/77
既然有了分块优化,这一步就可以在$O(\sqrt{n})$的时间内解决,我们只要在$O(\sqrt{n})$的时间内求出来$f(a,b,1)$就好了
我们发现$f(a,b,1)$是和$gcd$相关的函数,于是可以很自然的想到莫比乌斯反演的套路,不会套路的请先学习 http://sxysxy.org/blogs/77
按照套路,我们设
$\Large g(a,b,x)=\sum\limits_{x|gcd(i,j), i \le a, j \le b} \frac{i}{x} \cdot \frac{j}{x}$
于是就有
$\Large g(a,b,x)=\sum\limits_{i=1}^{ \lfloor \frac{min(a,b)}{x} \rfloor }f(a,b,i \cdot x)$
于是按照套路,我们有莫比乌斯反演
$\Large f(a,b,x)=\sum\limits_{i=1}^{ \lfloor \frac{min(a,b)}{x} \rfloor }\mu(i)g(a,b,i \cdot x)$
不会这个套路的请学习 http://sxysxy.org/blogs/77,或者学习WC2016课件,这个也可以从容斥的角度理解
同样对于$g(a,b,x)$,我们有
$\Large g(a,b,x)=g( \lfloor \frac{a}{x} \rfloor, \lfloor \frac{b}{x} \rfloor,1) \cdot x^2$
请读者思考这一步
那么就有
$\Large f(a,b,x)=\sum\limits_{i=1}^{ \lfloor \frac{min(a,b)}{x} \rfloor }\mu(i)g( \lfloor \frac{a}{i} \rfloor, \lfloor \frac{b}{i} \rfloor,x) \cdot i^2$
于是$\lfloor\frac{a}{i}\rfloor$和$\lfloor\frac{b}{i}\rfloor$也都是整除,又可以 前缀和+分块 优化,我们就能在$O(\sqrt{n})$的时间内求出每个$f(a,b,x)$啦
于是我们只要在$O(1)$的时间内求出$g(a,b,x)$,总时间复杂度就是$O(n)$啦!
显然
$\Large g(a,b,x)=g( \lfloor \frac{a}{x} \rfloor, \lfloor \frac{b}{x} \rfloor,1) \cdot x^2=\frac{ \lfloor \frac{a}{x} \rfloor \cdot ( \lfloor \frac{a}{x} \rfloor +1)}{2} \cdot \frac{ \lfloor \frac{b}{x} \rfloor \cdot ( \lfloor \frac{b}{x} \rfloor +1)}{2} \cdot x^2$
为什么显然请读者思考
于是任务大功告成!总时间复杂度为$O(n)$
实际上对于函数$g$和$f$,最后的那个参数$x$是没必要传的,因为在我们的分析中,$x$都被化作$1$了
上代码:
#include <cstring>
#include <cstdio>
#include <algorithm>
#define register
#define inline using namespace std;
typedef long long ll;
const int MOD = ;
const int MAXN = ; int n, m; bool vis[MAXN] = {};
int prime[MAXN], pidx = , miu[MAXN] = {}, pf[MAXN] = {};
void prelude_miu() {
miu[] = ;
for( register int i = ; i <= n; ++i ) {
if( !vis[i] ) {
prime[pidx++] = i;
miu[i] = -;
}
for( register int j = ; j < pidx; ++j ) {
register int k = i * prime[j];
if( k > n ) break;
vis[k] = true;
if( i % prime[j] ) miu[k] = -miu[i];
else break;
}
}
for( register int i = ; i <= n; ++i ) // 计算miu[i] * i^2的前缀和,注意取模
pf[i] = (pf[i-] + ((ll)miu[i] * i * i % MOD + MOD) % MOD) % MOD;
} inline int g( register int n, register int m ) { // 函数g
return (ll(n)*(n+)/ % MOD) * (ll(m)*(m+)/ % MOD) % MOD;
}
int f( register int n, register int m ) { // 函数f
register int ans = ;
for( register int i = , j; i <= n; i = j+ ) {
j = min( n/(n/i), m/(m/i) ); // 分块优化
// printf( "g(%d,%d) = %lld\n", n/j, m/j, g(n/j,m/j) );
ans = (ans + g(n/j,m/j) * ll((pf[j] - pf[i-] + MOD) % MOD) % MOD) % MOD;
}
return ans;
}
void solve() {
register int ans = ;
for( register int i = , j; i <= n; i = j+ ) { // 枚举最大公因数i
j = min( n/(n/i), m/(m/i) ); // 分块优化
// printf( "f(%d,%d) = %lld\n", n/j, m/j, f(n/j,m/j) );
ans = (ans + f(n/j,m/j) * (ll(i+j)*(j-i+)/ % MOD) % MOD) % MOD;
}
printf( "%d\n", ans );
} int main() {
scanf( "%d%d", &n, &m );
prelude_miu();
if( n > m ) swap(n,m); // 保证n < m,不用每次取min
solve();
return ;
}
【题解】Crash的数字表格 BZOJ 2154 莫比乌斯反演的更多相关文章
- Bzoj2154 Crash的数字表格 乘法逆元+莫比乌斯反演(TLE)
题意:求sigma{lcm(i,j)},1<=i<=n,1<=j<=m 不妨令n<=m 首先把lcm(i,j)转成i*j/gcd(i,j) 正解不会...总之最后化出来的 ...
- Crash的数字表格 BZOJ 2154 / jzptab BZOJ 2693
jzptab [问题描述] 求: 多组询问 [输入格式] 一个正整数T表示数据组数 接下来T行 每行两个正整数 表示N.M [输出格式] T行 每行一个整数 表示第i组数据的结果 [样例输入] 1 4 ...
- [luogu1829][bzoj2154][国家集训队]Crash的数字表格 / JZPTAB【莫比乌斯反演】
传送门:洛谷,bzoj 题目描述 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b)表示能同时整除a和b的最小正整 ...
- 洛谷P1829 [国家集训队]Crash的数字表格 / JZPTAB(莫比乌斯反演)
传送门 式子好麻烦orz……大佬好腻害orz->这里 //minamoto #include<iostream> #include<cstdio> #define ll ...
- 洛谷 P1829 [国家集训队]Crash的数字表格 / JZPTAB(莫比乌斯反演)
题意:求$\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)$. 开始开心(自闭)化简: $\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)$ =$\su ...
- 【莫比乌斯反演】关于Mobius反演与lcm的一些关系与问题简化(BZOJ 2154 crash的数字表格&&BZOJ 2693 jzptab)
BZOJ 2154 crash的数字表格 Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b ...
- bzoj 2154 莫比乌斯反演求lcm的和
题目大意: 表格中每一个位置(i,j)填的值是lcm(i,j) , 求n*m的表格值有多大 论文贾志鹏线性筛中过程讲的很好 最后的逆元我利用的是欧拉定理求解的 我这个最后线性扫了一遍,勉强过了,效率不 ...
- [BZOJ 2154]Crash的数字表格(莫比乌斯反演+数论分块)
[BZOJ 2154]Crash的数字表格(莫比乌斯反演+数论分块) 题面 求 \[\sum_{i=1}^{n} \sum_{j=1}^{m} \mathrm{lcm}(i,j)\] 分析 \[\su ...
- BZOJ 2154: Crash的数字表格 [莫比乌斯反演]
2154: Crash的数字表格 Time Limit: 20 Sec Memory Limit: 259 MBSubmit: 2924 Solved: 1091[Submit][Status][ ...
随机推荐
- Java学习 · 初识 面向对象深入一
面向对象深入 1.面向对象三大特征 a) 继承 inheritance 子类可以从父类继承属性和方法 子类可以提供自己的属性方法 b) 封装 encapsulation 对外隐藏某些属性和方法 对外公 ...
- Wordcount -- MapReduce example -- Mapper
Mapper maps input key/value pairs into intermediate key/value pairs. E.g. Input: (docID, doc) Output ...
- zabbix 2.2.2 安装部署
zabbix 2.2.2版本与1.8.3版本安装过程略有不同,下面为实施步骤: 服务端:172.16.1.61 客户端:172.16.1.8 搭建zbbix软件 安装LAMP环境及依赖包 [root@ ...
- 十 Writing YARN Applications
本节介绍: 使用yarn 高级提交写yarn应用程序.其实已经yarn底层API.MR计算框架对底层的API实现了封装. 高级提交指直接使用yarn的三种接口来提交应用程序: 1)YarnCl ...
- Rightmost Digit(最后一位数字)
Description Given a positive integer N, you should output the most right digit of N^N. Input The ...
- 王者荣耀交流协会beta冲刺贡献分分配结果
根据原来约定的(贡献分分配规则){http://www.cnblogs.com/gaoyb348/p/7828523.html} 总共20个任务,一共35分,那么一个任务1.75分. 冉华完成了9个任 ...
- 软件工程 part4 评价3作品 修改
作品1 抢答器 地址: https://modao.cc/app/ylGTXobcMU7ePNi6tY53gG4iraLl0md评价: 老师有意见,我重新评价,这个作品是我测试最方便的,地址点进去直接 ...
- Calculation PartⅡ
GitHub/object-oriented 误删内容--周末修复
- IOS开发NSBundle使用
bundle是一个目录,其中包含了程序会使用到的资源. 这些资源包含了如图像,声音,编译好的代码,nib文件(用户也会把bundle称为plug-in). 对应bundle,cocoa提供了类NSBu ...
- <Android>资源的访问,颜色、字符串、尺寸、XML、DRAWABLES资源分使用
1.资源的访问 代码中使用Context的getResources()方法得到Resources对象,访问自己定义的资源R.资源文件类型.资源文件名称,访问系统定义的资源android.R. 资源文件 ...