BZOJ_2693_jzptab_莫比乌斯反演

Description

Input

一个正整数T表示数据组数

接下来T行 每行两个正整数 表示N、M

Output

T行 每行一个整数 表示第i组数据的结果

Sample Input

1

4 5

Sample Output

122

HINT
T <= 10000

N, M<=10000000


$\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}lcm(i,j)$

$=\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}\frac{i*j}{gcd(i,j)}$

$=\sum\limits_{p=1}^{n}\sum\limits_{i=1}^{\lfloor\frac{n}{p} \rfloor}
\sum\limits_{j=1}^{\lfloor\frac{m}{p} \rfloor} i*j*p*[gcd(i,j)=1]$

$=\sum\limits_{p=1}^{n}p\sum\limits_{i=1}^{\lfloor\frac{n}{p} \rfloor}
\sum\limits_{j=1}^{\lfloor\frac{m}{p} \rfloor} i*j
\sum\limits_{d|gcd(i,j)}\mu(d)$

$=\sum\limits_{p=1}^{n}p
\sum\limits_{d=1}^{n/p}\mu(d)*d^{2}
\sum\limits_{i=1}^{\lfloor\frac{n/p}{d} \rfloor}
\sum\limits_{j=1}^{\lfloor\frac{m/p}{d} \rfloor} i*j
$

设$s[n]=\sum\limits_{i=1}^{n}i$

$=\sum\limits_{p=1}^{n}p
\sum\limits_{d=1}^{n/p}\mu(d)*d^{2}*
s[\lfloor\frac{n/p}{d} \rfloor]*
s[\lfloor\frac{m/p}{d} \rfloor]
$

设$Q=d*p,先枚举Q$

$=\sum\limits_{Q=1}^{n}
s[\lfloor\frac{n}{Q} \rfloor]*
s[\lfloor\frac{m}{Q} \rfloor]
\sum\limits_{d|Q}\mu(d)*d^{2}*\lfloor\frac{Q}{d} \rfloor
$

设$f[n]=\sum\limits_{d|n}\mu(d)*d^{2}*\lfloor\frac{n}{d} \rfloor
=n\sum\limits_{d|n}\mu(d)*d$

$=\sum\limits_{Q=1}^{n}
s[\lfloor\frac{n}{Q} \rfloor]*
s[\lfloor\frac{m}{Q} \rfloor]*f[Q]
$

$然后发现f[n]=n*g[n],g[n]为 id卷\mu 的积性函数$

$我们可以处理出f[n]的前缀和,然后O(\sqrt{n})处理即可$

$mdlswl$

BZOJ_2693_jzptab_莫比乌斯反演的更多相关文章

  1. hdu1695 GCD(莫比乌斯反演)

    题意:求(1,b)区间和(1,d)区间里面gcd(x, y) = k的数的对数(1<=x<=b , 1<= y <= d). 知识点: 莫比乌斯反演/*12*/ 线性筛求莫比乌 ...

  2. BZOJ 2154: Crash的数字表格 [莫比乌斯反演]

    2154: Crash的数字表格 Time Limit: 20 Sec  Memory Limit: 259 MBSubmit: 2924  Solved: 1091[Submit][Status][ ...

  3. BZOJ2301: [HAOI2011]Problem b[莫比乌斯反演 容斥原理]【学习笔记】

    2301: [HAOI2011]Problem b Time Limit: 50 Sec  Memory Limit: 256 MBSubmit: 4032  Solved: 1817[Submit] ...

  4. Bzoj2154 Crash的数字表格 乘法逆元+莫比乌斯反演(TLE)

    题意:求sigma{lcm(i,j)},1<=i<=n,1<=j<=m 不妨令n<=m 首先把lcm(i,j)转成i*j/gcd(i,j) 正解不会...总之最后化出来的 ...

  5. 莫比乌斯函数筛法 & 莫比乌斯反演

    模板: int p[MAXN],pcnt=0,mu[MAXN]; bool notp[MAXN]; void shai(int n){ mu[1]=1; for(int i=2;i<=n;++i ...

  6. 【BZOJ-2440】完全平方数 容斥原理 + 线性筛莫比乌斯反演函数 + 二分判定

    2440: [中山市选2011]完全平方数 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2371  Solved: 1143[Submit][Sta ...

  7. POI2007_zap 莫比乌斯反演

    题意:http://hzwer.com/4205.html 同hdu1695 #include <iostream> #include <cstring> #include & ...

  8. hdu.5212.Code(莫比乌斯反演 && 埃氏筛)

    Code Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submi ...

  9. CSU 1325 莫比乌斯反演

    题目大意: 一.有多少个有序数对(x,y)满足1<=x<=A,1<=y<=B,并且gcd(x,y)为p的一个约数: 二.有多少个有序数对(x,y)满足1<=x<=A ...

随机推荐

  1. 小议 HashMap

    大家都知道,在Java里对对象的操作是基于引用的.而当我们需要对一组对象操作的时候,就需要有接收这一组引用的容器.平时我们最常用的就是数组.在Java里可以定义一个对象数组来完成许多操作.可是,数组长 ...

  2. Centos7 时区的设置

    Linux 系统(我特指发行版, 没说内核) 下大部分软件的风格就是不会仔细去考虑向后 的兼容性, 比如你上个版本能用这种程序配置, 没准到了下一个版本, 该程序已经不见了. 比如 sysvinit ...

  3. Day19 Django

    老师代码博客: http://www.cnblogs.com/yuanchenqi/articles/7552333.html 上节内容回顾: class Book(models.Model): ti ...

  4. jsoup 使用总结1--添加header

    jsoup 使用总结1--添加header 大部分时候,我们使用jsoup解析网页的是否都是直接Jsoup.connect(url); 孰不知,我们可以添加header,例如 Host, Refere ...

  5. (转)java之Spring(IOC)注解装配Bean详解

    java之Spring(IOC)注解装配Bean详解   在这里我们要详细说明一下利用Annotation-注解来装配Bean. 因为如果你学会了注解,你就再也不愿意去手动配置xml文件了,下面就看看 ...

  6. jsonp跨域获取数据小解

    jsonp跨域获取数据小解 由于浏览器有同源策略,所以要想获取非同源(协议,域名,端口三者有一不同都算非同源)的页面的数据,就得进行跨域 (1) jsonp原理 由于script标签的src属性可以访 ...

  7. 部署网站: 配置项目到iis上运行报目录错误

    配置项目到iis上运行报目录错误 以下三种方法可使用: 1.添加文件访问权限 everyone  (线上环境慎用) 2.重新注册iis 3.web.config 加一段话 在<system.we ...

  8. ArcEngine小问题解决

    最近开始使用VS2012,在引用COM组件的时候,出现了无法嵌入互操作类型“……”,请改用适用的接口的错误提示. 找到解决方案:选中项目中引入的dll,鼠标右键,选择属性,把“嵌入互操作类型”设置为F ...

  9. 【转】asp.net获取当前页面的url地址

    设当前页完整地址是:http://www.jb51.net/aaa/bbb.aspx?id=5&name=kelli "http://"是协议名 "www.jb5 ...

  10. cug-day02

    昨天就主要就熟悉了项目,今天就简单完成了 后台代码,但是没有进行测试,前端的vue也不是很懂,所以还有很多需要学习的地方. 目前主要就是需要理解项目运行流程, 项目采用的是ssh ,用了vo和po , ...