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

\[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. mysql 忘记root密码怎么办?

    忘记root可以跳过grant table来登录 1.打开命令行输入以下命令 mysqld -nt --grant-skip-tables 2.在打开一个新命令行,输入以下命令可以登录, mysql ...

  2. 一键语法错误增强工具 ChineseErrorCorrector

    一键语法错误增强工具 欢迎使用我最近开源的使用一键语法错误增强工具,该工具可以进行14种语法错误的增强,不同行业可以根据自己的数据进行错误替换,来训练自己的语法和拼写模型,希望推动行业文本纠错的发展, ...

  3. mysql执行步骤口诀:发连缓分析,优化执行器

    mysql执行步骤口诀:发连缓分析,优化执行器 mysql执行步骤: 1.我们在客户端发起一个SQL的查询: 2.连接器判断用户登录以及用户权限: 3.缓存命中,走缓存(mysql缓存不是redis缓 ...

  4. Python列表、元组、字典和集合的用法

    1.列表 标志符号是[],元素可以修改.删除和新增 1.1 提取元素(索引从0开始计算) testList=['A','B','C',1,'D'] print(testList[1]) #打印索引区间 ...

  5. (一)MOngoDB的安装与配置(Windows和Centos)

    其实Windows的安装没有什么难的,主要是后的Centos下的安装,我是在centos7上安装测试的,不管是手动安装还是脚本安装我都是成功才发出来的,有什么的问题可以留言,写博客只是为了记录和简单的 ...

  6. 2018年视频,路径规划:层次化路径规划系统——hierarchical pathfinding system —— Hierarchical Dynamic Pathfinding for Large Voxel Worlds (续)

    前文: 2018年视频,路径规划:层次化路径规划系统--hierarchical pathfinding system -- Hierarchical Dynamic Pathfinding for ...

  7. java零基础到架构师学习线路(附视频教程)

    1.背景 很多人都在问,如何学java,要学那些内容,感觉学起来很痛苦,没得方向,学到什么程度才可以去找工作等, 在这里我以自己的学习经验工作经验和辅导学生的经验给大家梳理了一个学习线路和准备了我自己 ...

  8. 海豚调度调优 | 如何解决任务被禁用出现的Bug

    本系列文章是 DolphinScheduler 由浅入深的教程,涵盖搭建.二开迭代.核心原理解读.运维和管理等一系列内容.适用于想对 DolphinScheduler了解或想要加深理解的读者. 祝开卷 ...

  9. 【故障排查】10分钟解决Quartz重复调度的疑难杂症

    我司使用Apache DolphinScheduler作为调度框架很久了,感兴趣的小伙伴可以看看这些干货文章: 因为之前监控到会出现重复的调度的问题,所以此文记录排查重复调度问题的全过程,希望对社区其 ...

  10. # Apache SeaTunnel 究竟是什么?

    作者 | Shawn Gordon 翻译 | Debra Chen 原文链接 | What the Heck is Apache SeaTunnel? 我在2023年初开始注意到Apache SeaT ...