落谷 P3327 [SDOI2015]约数个数和 (莫比乌斯反演+分块)*

题目链接:https://www.luogu.org/problem/P3327题目描述

在这里插入图片描述
设d(x)为x的约数个数,给定N、M,求 ∑i=1N∑j=1Md(ij)\sumN_{i=1}\sumM_{j=1}d(ij)∑i=1N​∑j=1M​d(ij)
输入格式

输入文件包含多组测试数据。第一行,一个整数T,表示测试数据的组数。接下来的T行,每行两个整数N、M。
输出格式

T行,每行一个整数,表示你所求的答案。
输入输出样例
输入 #1

2
7 4
5 6

输出 #1

110
121

说明/提示

1<=N, M<=50000

1<=T<=50000


分析:

求 ∑ i = 1 N ∑ j = 1 M d ( i j ) 求\sum_{i=1}^N\sum_{j=1}^M{d(ij)} i=1Nj=1Md(ij)
对于d函数有个性质,在最上面题目链接页面中的解题报告中有它的证明,不过能力不够,不太理解,先记着吧
d ( i j ) = ∑ x ∣ i ∑ y ∣ j [ g c d ( x , y ) = 1 ] d(ij)=\sum_{x|i}\sum_{y|j}{[gcd(x,y)=1]} d(ij)=xiyj[gcd(x,y)=1]
注:[gcd(x,y)=1]表示gcd(x,y)等于1时值为1,其余情况为0
那么所求为
∑ i = 1 N ∑ j = 1 M ∑ x ∣ i ∑ y ∣ j [ g c d ( x , y ) = 1 ] \sum_{i=1}^N\sum_{j=1}^M{\sum_{x|i}\sum_{y|j}{[gcd(x,y)=1]}} i=1Nj=1Mxiyj[gcd(x,y)=1]
上式中,假如以x,y为自变量,i和j满足分别是x,y的倍数,对应的个数分别为 N x , M y \frac{N}{x},\frac{M}{y} xN,yM(注:像这种除法的,若没有特别说明,都表示向下取整),由此改变枚举顺序得
∑ x = 1 N ∑ y = 1 M [ g c d ( x , y ) = 1 ] ⋅ N x ⋅ M y \sum_{x=1}^N\sum_{y=1}^M{[gcd(x,y)=1]\cdot \frac{N}{x}\cdot \frac{M}{y} } x=1Ny=1M[gcd(x,y)=1]xNyM
定义f(n)函数,和g(n)函数:
f ( n ) : = ∑ x = 1 N ∑ y = 1 M [ g c d ( x , y ) = n ] ⋅ N x ⋅ M y g

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值