正在实现,不知道对不对,但是先放这,哪个大佬发现问题了和我说下

\[f(l)=\sum\cdots\sum[\gcd=1,\text{lcm}=l]
\]
\[g(l)=\sum\cdots\sum[\gcd=1,\text{lcm}\mid l]
\]
\[h(l)=\sum\cdots\sum[\text{lcm}\mid l]
\]

\[g(l)=\sum_{l\mid d}f(d)
\]
\[f(l)=\sum_{l\mid d}\mu(\frac{d}{l})g(d)
\]
\[h(l)=\sum_{d}g(d)=\sum_{1\mid d}g(d)
\]
\[g(l)=\sum_{1\mid d}\mu(d)h(d)=\sum_{d}\mu(d)h(d)
\]
\[f(l)=\sum_{l\mid d}\mu(\frac{d}{l})\sum_{e}\mu(e)h(e)
\]

UPD:确实不对,下面这版是对的

\[f(g,l)=\sum\cdots\sum[\gcd=g,\text{lcm}=l]
\]
\[g(g,l)=\sum\cdots\sum[\gcd=g,\text{lcm}\mid l]
\]
\[h(l)=\sum\cdots\sum[\text{lcm}\mid l]
\]

\[g(1,l)=\sum_{d\mid l}f(1,d)
\]
\[f(1,l)=\sum_{d\mid l}\mu(d)g(1,\lfloor{\frac{l}{d}}\rfloor)
\]
\[\color{white}h(l)=\sum_{i}g(i,l)\rightarrow \sum\cdots\sum[\text{lcm}\mid l]\rightarrow l=k\text{lcm},\forall\gcd=i\rightarrow l=\frac{1}{i}k\text{lcm}\rightarrow il=k\text{lcm}\rightarrow h(l)=\sum_{i}g(1,il)
\]

内层求法:设

\[l=\sum_{i}p_{i}^{a_{i}}
\]

\[h(l)=C^{n}_{\sum_{i}a_{i}}
\]

外层求法:

\[\gcd\mid m
\]
\[\text{lcm}=m-\gcd
\]
\[\begin{cases}\gcd=g\\\text{lcm}=m-g\end{cases}\Longleftrightarrow\begin{cases}\gcd=1\\\text{lcm}=\frac{m-g}{g}\end{cases}
\]
\[\text{ans}=\sum_{k\mid m}f(\frac{m-k}{k})
\]

总之就是 先处理外层,再分解因子,提前处理好 \(\sum_{e}\mu(e)h(e)\),\(h(e)\) 内部分解质因数求和,再统计答案即可. 复杂度有点吓人

T2 的莫反式子的更多相关文章

  1. Codeforces 548E(莫反、容斥)

    转化为质数域上的操作,如果用莫反的话,记录因数的cnt. 其实莫反的推式子最后和容斥做法殊途同归了,容斥的系数就是莫比乌斯函数. const int maxn = 2e5 + 5, maxa = 5e ...

  2. Codeforces 585E - Present for Vitalik the Philatelist(简单莫反+狄利克雷前缀和)

    Codeforces 题目传送门 & 洛谷题目传送门 一道不算太难的 D1E 罢--虽然我不会做/kk u1s1 似乎这场 Div1 挺水的?F 就是个 AC 自动机板子还被评到了 3k2-- ...

  3. NFLSOJ #917 -「lych_cys模拟题2018」橘子树(树剖+ODT+莫反统计贡献的思想+动态开点线段树)

    题面传送门 sb 出题人不在题面里写 \(b_i=0\) 导致我挂成零蛋/fn/fn 首先考虑树链剖分将路径问题转化为序列上的问题,因此下文中简称"位置 \(i\)"表示 DFS ...

  4. Codeforces 1603D - Artistic Partition(莫反+线段树优化 dp)

    Codeforces 题面传送门 & 洛谷题面传送门 学 whk 时比较无聊开了道题做做发现是道神题( 介绍一种不太一样的做法,不观察出决策单调性也可以做. 首先一个很 trivial 的 o ...

  5. 难度并不NOIP的NOIP模拟赛

    今天老师请了前几届的学长来讲课,可是讲课为什么要考试呢... 学长说难度是NOIP,于是我就naive的跟着参加了,然而T3难度并不友好,感觉确实不是很适合我们现在做......不过课本来也不是给我们 ...

  6. bjwc Day1 暴力大战

    今天终于有题了... 题目是COCI2016/2017 Round #4 T1一看就是NP问题,k<=50,开始想暴力,想了个n^4的,大概能过,就没去管它 T2想得太naive,丢了100分给 ...

  7. NOI 2016 Day1 题解

    今天写了NOI2016Day1的题,来写一发题解. T2 网格 题目传送门 Description \(T\) 次询问,每次给出一个 \(n\times m\) 的传送门,上面有 \(c\) 个位置是 ...

  8. ZJOI 2019 游记

    一轮 现在才知道大家都写了2333,现在补上不晚吧. Day 0 跟Crossing打了一路的王者,丝毫没有困意.颁奖仪式看到rank 1 又是xj的zyz,QWQ.被冯缘的热情四射的演讲给吓到了.然 ...

  9. HZNU ACM一日游 2019.3.17 【2,4,6-三硝基甲苯(TNT)】

    Travel Diary 早上8:00到HG,听说hjc20032003在等我. 然后他竟然鸽我...最后还是勉强在8:30坐上去偏僻的HZNU的地铁. 到文新,然后带上fjl,打滴滴,一行人来到了H ...

  10. 题解-bzoj2154Crash的数字表格 & bzoj2693 jzptab

    Problem bzoj2818-单组询问-无权限 bzoj2693-多组询问-需权限 洛谷1829-单组询问-无权限 \(T\)组询问(如果有),给定 \(n,m\),求 \[\sum_{i=1}^ ...

随机推荐

  1. Known框架实战演练——进销存框架搭建

    本文介绍如何使用Known开发框架搭建进销存管理系统的项目结构,以及开发前的一些配置和基础代码. 项目代码:JxcLite 开源地址: https://gitee.com/known/JxcLite ...

  2. 一文全解:LVM(逻辑卷管理器)

    前两篇文章已经讲了关于磁盘分区和磁盘阵列的相关内容: 一文全懂:Linux磁盘分区 一文全懂:独立冗余磁盘阵列(RAID) 但是磁盘分区完后再想扩容或者缩容就比较麻烦了,甚至很多时候不能扩容或者缩容, ...

  3. java开发,json转list集合,原生实现

    java 是一门面象对象的语言,对象需要先定义,但是在外理网络请求时候会用到json 转成java 对象,虽然现代开发框架中也提供了很多工具和方法直接转换, 但是作为学习者了解 一下底层实现,更能灵活 ...

  4. Python 实现行为驱动开发 (BDD) 自动化测试详解

    ​ 在当今的软件开发领域,行为驱动开发(Behavior Driven Development,BDD)作为一种新兴的测试方法,逐渐受到越来越多开发者的关注和青睐.Python作为一门功能强大且易于使 ...

  5. Scratch作品-巴黎2024奥运会

    ​ <Scratch作品-巴黎2024奥运会>是一款以巴黎2024年奥运会为主题的互动作品,专为儿童和青少年设计.通过Scratch编程语言,这个作品生动地再现了奥运会的精彩瞬间,结合了动 ...

  6. DNS在架构中的使用

    1 介绍 DNS(Domain Name System,域名系统)是一种服务,它是域名和IP地址相互映射的一个分布式数据库,能够使人更方便的访问互联网,而不用去记住能够被机器直接读取的IP地址数串. ...

  7. vue3源码学习api-vue-sfc文件编译

    vue 最有代表性质的就是.VUE 的文件,每一个vue文件都是一个组件,那么vue 组件的编译过程是什么样的呢 Vue 单文件组件 (SFC)和指令 ast 语法树 一个 Vue 单文件组件 (SF ...

  8. 12、SpringMVC之拦截器

    12.1.环境搭建 创建名为spring_mvc_interceptor的新module,过程参考9.1节和9.5节 12.1.1.页面请求示例 <a th:href="@{/test ...

  9. 【SpringMVC】IDEA 不识别webapp的解决办法

    参考地址: https://blog.csdn.net/omrleft123/article/details/70237205

  10. 公共的网络云盘的存储真的安全吗?—— 百度云盘上的PDF文件无故被改名

    在百度云盘上上传了一个PDF文件,内容: 本来是没有啥问题的,但是今天使用百度云盘发现这个PDF文件居然被改名,被取消掉了扩展名: 简直是离谱离了一个大谱,太可怕了,看来这种公共云盘真的不太靠谱,虽然 ...