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

\[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. Oracle 序列学习与使用总结

    Oracle序列学习与使用总结 by:授客 QQ:1033553122 简述 序列是oracle提供的用于生成一系列数字的数据库对象,序列会自动生成顺序递增的序列号,可用于提供唯一的自动递增主键.序列 ...

  2. Docker安装及操作

    目录 docker 安装: 官方文档方法 CentOS Ubuntu docker-compose 单独安装 centos7 ubuntu22.04 docker 容器操作: docker启动与停止 ...

  3. 【MySQL】二级MySQL考试 救场帮助表

    周六去考二级,应用第一题就是添加外键约束 草,写了半天老说语法不对,然后急中生智,觉得默认的库里应该有文档说明表 以下是SQL查询过程: -- 猜测是在mysql库里面 mysql> USE m ...

  4. 【Java】自制查找工具

    需求:查找后台代码中写的SQL是否包含拆分表,如果存在,则返回那些表名 Context.txt 粘贴我们找到的DAO层代码,因为所有方法封装的SQL都在DAO层里[就理解为筛查的字符串] Dictio ...

  5. 【Vue】Re21 VueX 第二部分(Mutations)

    一.Mutations携带参数处理 Store状态的更新唯一方式:提交Mutations Mutations包含两部分: 1.字符串的事件类型[TYPE] 2.一个回调函数[HANDLER] 该函数的 ...

  6. Google的TPU的逐元素操作特性 —— TPU数据操作时最基本的操作原语(primitive)

    地址: https://jax.readthedocs.io/en/latest/pallas/tpu.html TPU是原来进行矩阵计算加速的一个科学计算的设备,功能大致等同于NVIDIA的GPU的 ...

  7. 对于强化学习算法中的AC算法(Actor-Critic算法) 的一些理解

    AC算法(Actor-Critic算法)最早是由<Neuronlike Adaptive Elements That Can Solve Difficult Learning Control P ...

  8. baselines算法库common/tile_images.py模块分析

    该模块只有一个函数,全部内容: import numpy as np def tile_images(img_nhwc): """ Tile N images into ...

  9. git clone 如何通过proxy进行远程代码仓库拷贝下载

    git使用proxy的方式和ssh的情况是差不多的,给出借鉴: SSH如何通过proxy进行服务器连接 ------------------------------------------------ ...

  10. SMU Spring 2023 Contest Round 4(第 21 届上海大学程序设计联赛 春季赛)

    A. Antiamuny wants to learn binary search 签到题. #include <map> #include <set> #include &l ...