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

\[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. Day 11 - 模拟考

    WTP 的大洗牌 题目描述 \(\text{input1}\) 3 1 1 1 1 2 3 \(\text{output1}\) 10 0 数据范围 shuffle.zip 第一题 题目描述 输入格式 ...

  2. Java JVM——13. 垃圾回收相关算法

    1.生存还是死亡? 在堆里存放着几乎所有的 Java 对象实例,在 GC 执行垃圾回收之前,首先需要区分出内存中哪些是存活对象,哪些是已经死亡的对象.只有被标记为己经死亡的对象,GC 才会在执行垃圾回 ...

  3. ElementUI FORM结合Vue实现横向排列表单项

    结合Vue实现横向排列表单项 前言 默认的,ElementUI的Form表单项(控件)是垂直排列,即一行显示一个表单项.但是在实际应用中,有时候会需要一行显示多个表单项.针对这类需求,笔者提供以下解决 ...

  4. 题解:P7482 不条理狂诗曲

    题解:P7482 不条理狂诗曲 本题解借鉴 blossom_j 大佬思路,但这位大佬的题解似乎没放正确代码. 题意 对于每一个 \(a\) 的子区间 \(a_{l\dots r}\),求选择若干个不连 ...

  5. Spectre.Console.Cli注入服务的几种姿势

    Spectre.Console大家可能都不陌生,写控制台程序美化还是不错的,支持着色,表格,图标等相当nice,如果对这个库不熟悉我强烈推荐你了解一下,对于写一些CLI小工具还是相当方便的, 本文主要 ...

  6. LeetCode654. 最大二叉树

    题目链接:https://leetcode.cn/problems/maximum-binary-tree/description/ 题目叙述 给定一个不重复的整数数组 nums . 最大二叉树 可以 ...

  7. 分段树(segment tree)的实现 —— 强化学习中 "优先级回放机制" 的重要组成部分

    分段树(segment tree)是强化学习中 "优先级回放机制" 的重要组成部分.本文针对分段树(segment tree)的一个开源版本的实现来进行分析,代码地址: https ...

  8. 安卓如何设置开机启动某个程序?init.rc给你搞定

    一.如何设置开机启动某个程序? 1.需求描述 最近有个项目需要在Android开机启动之后,自动执行一个C语言编写的程序:pengd 该程序运行时需要修改网络ip地址及其他网络操作,所以需要root权 ...

  9. Figma 替代品 Penpot 安装和使用教程

    在设计领域,Figma 无疑是一个巨人.它彻底改变了设计流程,将协作带到了一个全新的高度.然而,随着 Adobe 收购 Figma 的消息传出,许多设计师和开发者开始担心:Figma 未来会如何演变? ...

  10. MarginNote 4 内存泄露?

    在床上用电脑的时候突然发现电脑风扇呼呼响,一摸很烫,以为是被子把出风口堵住了,于是调整角度继续用.结果一段时间之后风扇还是狂转不停,然后收到了这样的提示.不看不知道一看吓一跳,MarginNote 4 ...