[BZOJ2693]jzptab
Description
求\(\sum\limits_{i=1}^n\sum\limits_{j=1}^m lcm(i,j)\),答案模1e9+9输出,多组询问
Input
一个正整数T表示数据组数
接下来T行 每行两个正整数 表示N、M
Output
T行 每行一个整数 表示第i组数据的结果
Sample Input
1
4 5
Sample Output
122
HINT
T <= 10000
N, M<=10000000
我们令n<m,然后将柿子化简
\[\sum\limits_{i=1}^n\sum\limits_{j=1}^m lcm(i,j)\]
\[\sum\limits_{i=1}^n\sum\limits_{j=1}^m \dfrac{i\times j}{\gcd(i,j)}\]
\[\sum\limits_{d=1}^n\sum\limits_{i=1}^n\sum\limits_{j=1}^m\dfrac{i\times j}{d}[\gcd(i,j)=d]\]
\[\sum\limits_{d=1}^n\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum\limits_{j=1}^{\lfloor\frac{m}{d}\rfloor}\dfrac{d^2\times i\times j}{d}[\gcd(i,j)=1]\]
\[\sum\limits_{d=1}^n d\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum\limits_{j=1}^{\lfloor\frac{m}{d}\rfloor}i\times j\sum\limits_{x|i,x|j}\mu(x)\]
\[\sum\limits_{d=1}^n d\sum\limits_{x=1}^{\lfloor\frac{n}{d}\rfloor}\mu(x)\times x^2\sum\limits_{i=1}^{\lfloor\frac{n}{dx}\rfloor}\sum\limits_{j=1}^{\lfloor\frac{m}{dx}\rfloor}i\times j\]
我们发现最后那个是等差数列,继续化简
\[\sum\limits_{d=1}^n d\sum\limits_{x=1}^{\lfloor\frac{n}{d}\rfloor}\mu(x)\times x^2\dfrac{\lfloor\frac{n}{dx}\rfloor(\lfloor\frac{n}{dx}\rfloor+1)\lfloor\frac{m}{dx}\rfloor(\lfloor\frac{m}{dx}\rfloor+1)}{4}\]
然后我们令\(T=dx\),那么得到
\[\sum\limits_{T=1}^n\dfrac{\lfloor\frac{n}{T}\rfloor(\lfloor\frac{n}{T}\rfloor+1)\lfloor\frac{m}{T}\rfloor(\lfloor\frac{m}{T}\rfloor+1)}{4}T\sum\limits_{x|T}\mu(x)x\]
我们设\(f(T)=\sum\limits_{x|T}\mu(x)x\),预处理出f,就可以分块了
设\(g(x)=\mu(x)x\),当a,b互质,\(g(a)\times g(b)=ab\mu(a)\mu(b)=ab\mu(ab)=g(ab)\),所以g是积性函数,根据莫比乌斯反演的性质,f也是积性函数
令\(T=\prod\limits_{i=1}^k P_i^{x_i}\),\(f(P_i^{x_i})=(1-P_i)\),那么
\[f(T)=\prod\limits_{i=1}^k(1-P_i)\]
这样子我们可以在\(O(n)\)时间内线筛出来,然后维护一下\(T\times f(T)\)的前缀和,然后就可以在\(O(\sqrt N)\)的时间内完成每次询问
/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline int read(){int x=0,f=1;char ch=getchar();for (;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') f=-1;for (;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+ch-'0';return x*f;
}
inline void print(int x){if (x>=10) print(x/10);putchar(x%10+'0');
}
const int N=1e7,p=1e8+9;
int prime[N+10],f[N+10],tot;
bool inprime[N+10];
void prepare(){f[1]=1;for (int i=2;i<=N;i++){if (!inprime[i]) prime[++tot]=i,f[i]=1-i+p;for (int j=1;j<=tot&&i*prime[j]<=N;j++){inprime[i*prime[j]]=1;if (i%prime[j]==0){f[i*prime[j]]=f[i];break;}f[i*prime[j]]=1ll*f[i]*f[prime[j]]%p;}}for (int i=1;i<=N;i++) f[i]=(f[i-1]+1ll*i*f[i]%p)%p;
}
int get(int x){return (1ll*x*(x+1)>>1)%p;}
int main(){prepare();for (int Data=read();Data;Data--){int n=read(),m=read(),pos,Ans=0;if (n>m) swap(n,m);for (int T=1;T<=n;T=pos+1){pos=min(n/(n/T),m/(m/T));Ans=(Ans+1ll*get(n/T)*get(m/T)%p*(f[pos]-f[T-1]+p)%p)%p;}printf("%d\n",Ans);}return 0;
}
转载于:https://www.cnblogs.com/Wolfycz/p/9489894.html
[BZOJ2693]jzptab相关推荐
- BZOJ2154: Crash的数字表格 BZOJ2693: jzptab
[传送门:BZOJ2154&BZOJ2693] 简要题意: 给出n,m,求$\sum_{i=1}^{n}\sum_{j=1}^{m}LCM(i,j)$ 题解: 莫比乌斯反演(因为BZOJ269 ...
- 【莫比乌斯反演】[HYSBZ/BZOJ2693]jzptab
题目大意就是求∑ni=1∑mj=1lcm(i,j),但是有多组输入数据. 如果之前做过[莫比乌斯反演][HYSBZ\BZOJ2154]Crash的数字表格,就会发现,对于每一个询问,有O(n)的做法, ...
- 【bzoj2693】jzptab 莫比乌斯反演+线性筛
题目描述 输入 一个正整数T表示数据组数 接下来T行 每行两个正整数 表示N.M 输出 T行 每行一个整数 表示第i组数据的结果 样例输入 1 4 5 样例输出 122 题解 莫比乌斯反演+线性筛 由 ...
- BZOJ 2154 [国家集训队]Crash的数字表格 / JZPTAB(莫比乌斯反演,经典好题)(Luogu P1829)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 P1829 [国家集训队]Crash的数字表格 / JZPTAB(反演,经典好题) Problem S ...
- #2693. jzptab
jzptab ∑i=1n∑j=1mlcm(i,j)∑i=1n∑j=1mijgcd(i,j)∑d=1nd∑i=1nd∑j=1mdij[gcd(i,j)=1]∑d=1nd∑k=1ndk2μ(k)∑i= ...
- P1829 [国家集训队]Crash的数字表格 / JZPTAB
P1829 [国家集训队]Crash的数字表格 / JZPTAB 题意: 求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)∑i=1n∑j ...
- P1829 [国家集训队]Crash的数字表格 / JZPTAB(莫比乌斯反演)
[国家集训队]Crash的数字表格 / JZPTAB 题目描述 今天的数学课上,Crash 小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数 a a a 和 b b ...
- P1829-[国家集训队]Crash的数字表格/JZPTAB【莫比乌斯反演】
正题 题目链接:https://www.luogu.com.cn/problem/P1829 题目大意 给出n,mn,mn,m求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^n\sum_{ ...
- 【数论】Crash的数字表格 / JZPTAB(P1829)
正题 P1829 题目大意 给出n,m,求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)i=1∑nj=1∑mlcm(i,j) 解题思路 ∑i= ...
最新文章
- 在Eclipse中导入Java程序
- nginx 访问日志分析工具 goacess
- python是什么怎么用-什么是python items函数?怎么使用它?
- 洛谷3933 Chtholly Nota Seniorious 二分答案+贪心
- 初探linux子系统集之led子系统(一)【转】
- Shell的sort、uniq、tr、cut、命令和 正则表达式
- Java垃圾回收(4)
- mybatis generator生成example_[Springboot系列] SpringBoot与Mybatis结合
- vb.net中递归退到最外层_数组中的逆序对
- 《如何搭建小微企业风控模型》第六节节选 特征工程(下)
- 阿里云蒋江伟:什么是真正的云原生?
- 微信小程序转盘实现(真心话大冒险)
- mysql 拼音模糊查询_mysql中文字段转拼音首字母,以及中文拼音模糊查询
- 小丸工具箱\FFmpeg踩坑记录
- 这一年,我所阐释的Linux哲学
- 邮政挂号信和包裹查询网站(非EMS)
- c语言编程一张纸对折,【动手试试】把一张纸对折
- 计算机视觉注意力机制-Attention
- 集成微控制器使太阳能微型逆变器设计成本有效
- Unity Editor Window Zooming