BZOJ DZY Loves Math系列
⑤(BZOJ 3560)
$\Sigma_{i_1|a_1}\Sigma_{i_2|a_2}\Sigma_{i_3|a_3}\Sigma_{i_4|a_4}...\Sigma_{i_n|a_n}\phi(i_1i_2i_3i_4...i_n)$
$\phi()$是积性函数
$\phi(p^k)=p^{k-1}*(p-1)$
设当前质数为p,对于第i个数,假设它分解质因数后p的次数为ai,那么p的答案就是
$[(1+p^1+...+p^{a1})(1+p^1+...+p^{a2})...(1+p^1+...+p^{an})-1]\frac{p-1}{p}+1$
乘起来就好了.....
$\Sigma_{i=1}^n\Sigma_{j=1}^mlcm(i,j)^{gcd(i,j)}$
$=\Sigma_{i=1}^n\Sigma_{j=1}^m (\frac{i*j}{gcd(i,j)})^{gcd(i,j)}$
枚举gcd(i,j)=d
$=\Sigma_{d=1}^n\Sigma_{i=1}^{\lfloor \frac{n}{d}\rfloor}\Sigma_{j=1}^{\lfloor \frac{m}{d}\rfloor}(d*i*j)^d*(gcd(i,j)==1)$
$=\Sigma_{d=1}^n\Sigma_{i=1}^{\lfloor \frac{n}{d}\rfloor}\Sigma_{j=1}^{\lfloor \frac{m}{d}\rfloor}\Sigma_{k|i且k|j}(d*i*j)^d$
$=\Sigma_{d=1}^nd^d\Sigma_{t=1}^{\lfloor\frac{n}{d}\rfloor}\mu(t)[\Sigma_{i=1}^{\lfloor\frac{n}{dt}\rfloor}(it)^d\Sigma_{j=1}^{\lfloor\frac{m}{dt}\rfloor}(jt)^d]$
$=\Sigma_{d=1}^nd^d\Sigma_{t=1}^{\lfloor\frac{n}{d}\rfloor}\mu(t)*t^{2d}[\Sigma_{i=1}^{\lfloor\frac{n}{dt}\rfloor}i^d\Sigma_{j=1}^{\lfloor\frac{m}{dt}\rfloor}j^d]$
$\Sigma _{i=1}^n\Sigma _{j=1}^i\mu(lcm(i,j)^{gcd(i,j)})$
$=\Sigma_{k=1}^n\Sigma_{i=1}^{\lfloor\frac{n}{k}\rfloor}\Sigma_{j=1}^i\mu((ijk)^{k}*gcd(i,j)==1)$
$∵$k>1时 $\mu(x^k)=0$
$∴ =\Sigma_{i=1}^n\Sigma_{j=1}^i\mu(ij)*e(gcd(i,j))$
$∵gcd(i,j)==1$
$∴\mu(ij)=\mu(i)*\mu(j)$
$=\Sigma_{i=1}^n\mu(i)*\Sigma_{j=1}^i\mu(j)*\Sigma_{k|i且k|j}\mu(k)$
$=\Sigma_{i=1}^n\mu(i)*\Sigma_{k|i}\mu(k)\Sigma_{j=1}^{\lfloor\frac{i}{k}\rfloor}\mu(jk)$
$\mu(i)≠0$时 再枚举k是i的约数 发现数量只有$5*10^7$
复杂度变成了
什么复杂度
O(能过)就好了...
BZOJ DZY Loves Math系列的更多相关文章
- [BZOJ] DZY Loves Math 系列 I && II
为了让自己看起来有点事干 ,做个套题吧..不然老是东翻翻西翻翻也不知道在干嘛... \(\bf 3309: DZY \ Loves \ Math\) 令 \(h=f*\mu\) 很明显题目要求的就是\ ...
- DZY Loves Math 系列详细题解
BZOJ 3309: DZY Loves Math I 题意 \(f(n)\) 为 \(n\) 幂指数的最大值. \[ \sum_{i = 1}^{a} \sum_{j = 1}^{b} f(\gcd ...
- DZY Loves Math系列
link 好久没写数学题了,再这样下去吃枣药丸啊. 找一套应该还比较有意思的数学题来做. [bzoj3309]DZY Loves Math 简单推一下. \[\sum_{i=1}^n\sum_{j=1 ...
- bzoj DZY Loves Math V
Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 509 Solved: 284[Submit][Status][Discuss] Descriptio ...
- ●BZOJ 3309 DZY Loves Math
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3309 题解: 莫比乌斯反演,线筛 化一化式子: f(x)表示x的质因子分解中的最大幂指数 $ ...
- BZOJ 3561 DZY Loves Math VI
BZOJ 3561 DZY Loves Math VI 求\(\sum_{i=1}^{n}\sum_{j=1}^{m}\text{lcm}(i,j)^{\gcd(i,j)}\),钦定\(n\leq m ...
- BZOJ 3309: DZY Loves Math
3309: DZY Loves Math Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 761 Solved: 401[Submit][Status ...
- 【BZOJ】3309: DZY Loves Math 莫比乌斯反演优化
3309: DZY Loves Math Description 对于正整数n,定义f(n)为n所含质因子的最大幂指数.例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007) ...
- BZOJ 3512: DZY Loves Math IV [杜教筛]
3512: DZY Loves Math IV 题意:求\(\sum_{i=1}^n \sum_{j=1}^m \varphi(ij)\),\(n \le 10^5, m \le 10^9\) n较小 ...
随机推荐
- elasticsearch数据库使用
elasticsearch的一个最为显著的优点:快速全文检索.关于elasticsearch 全文检索的原理,请看:https://blog.csdn.net/wolfcode_cn/article/ ...
- MVC和MVVM的区别
现在是市面上使用MVVM模型越来越多,相关的前端框架的选择也相应增多. 那么什么是MVVM模型,它和传统MVC模型有什么区别,接下来我们来总结探讨一下. MVVM即Model-View-ViewMod ...
- MySQL基于域名做高可用切换(Consul初试)
一,Consul功能介绍 服务发现 - Consul的客户端可用提供一个服务,比如 api 或者mysql ,另外一些客户端可用使用Consul去发现一个指定服务的提供者.通过DNS或者HTTP应用程 ...
- MySQL Workbench查看和修改表字段的Comment值
查看: 选择单个表->[右键]->[Table Inspector] 再选择Columns选项卡即可,把表格拉倒最后一列. 编辑: 选择单个表->[右键]->[Alter Ta ...
- SpringMVC断言--Assert
Web 应用在接受表单提交的数据后都需要对其进行合法性检查,如果表单数据不合法,请求将被驳回.类似的,当我们在编写类的方法时,也常常需要对方法入参进行合 法性检查,如果入参不符合要求,方法将通过抛出异 ...
- 1.4-动态路由协议OSPF②
LAB2.通过反掩码控制有哪些接口,在运行OSPF ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ step1:启动OSPF,并宣告网络: R1(config ...
- 重啓ubuntu后 VNC 自動運行
Vino-Server是Ubuntu自带的有個缺点:重启后不能自動運行(可能是基於安全吧!) 親身測試对象:windows & ubuntu 10.04已安装图形桌面gnome ***wind ...
- wikioi 1029 中序遍历总数
题意:给出前序遍历和后序遍历,求总共同拥有多少种中序遍历的可能. 思路: 对于一个节点.当且仅当它仅有一棵子树时,在保证先序和后序同样的前提下,才可能有不同的中序(它的子树可在左或右,所以有2种): ...
- 如何删除 mac 系统信息下的安装记录
方法:在finder中:Command+Shift+G 打开 /Library/Receipts/InstallHistory.plist 所有的记录都在 InstallHistory.plist 文 ...
- Qt Quick综合实例之文件查看器
假设你基于Qt SDK 5.3.1来创建一个Qt Quick App项目,项目模板为你准备的main.qml文档的根元素是ApplicationWindow或Window.这次我们就以Applicat ...