当前位置: 首页 > news >正文

P6210 「SWTR-4」Easy Math Problems 莫比乌斯反演(不完整.没推完)

没推完,以后再补充……

题目链接:https://www.luogu.com.cn/problem/P6210


定义 \(g(n, m, c)\) 表示 \([1, m]\) 中存在多少个整数 \(i\) 满足 \(\gcd(i, n) \le c\)

则第二问的答案为 \(g(n, r, c) - g(n, l-1, c)\)

现在考虑如何求解 \(g(n, m, c)\)

\[g(n, m, c) = \sum_{i=1}^{m} [ \gcd(i, n) \le c ] \]

考虑枚举 \(d = \gcd(i, n)\),则上式转换为

\[g(n, m, c) = \sum_{d \mid n, 1 \le d \le c} \sum_{i=1}^m [\gcd(i,n) = d] \]

\(d\) 从 艾弗森表达式(中括号)中提取出来,得

\[g(n, m, c) = \sum_{d \mid n, 1 \le d \le c} \sum_{i=1}^{m/d} [\gcd(i,n) = 1] \]

接下来我们先化简一下第二个 sigma。

由 莫比乌斯反演的一个常用的性质(\([x = 1] = \sum_{e \mid x} \mu(e)\)

我们可以将上式的最后一个 sigma 转换一下

\[\sum_{i=1}^{m/d} [\gcd(i,n) = 1] \]

\[= \sum_{i=1}^{m/d} \sum_{e \mid \gcd(i, n)} \mu(e) \]

我们考虑枚举 \(e\),再考虑 \(i\)

由于 \(e \mid \gcd(i, n)\),所以 \(e \mid i\) \(\Rightarrow\) \(i\)\(e\) 的倍数;

同时 \(i \le m / d\)

所以,满足条件的 \(i\) 的个数为 \(\left\lfloor \frac{m/d}{e} \right\rfloor = m / (de)\)

所以,上面那个 sigma 转为

\[\sum_{e \mid n, e \le m/d} (m / (de)) \mu(e) \]

代入原式得

\[g(n, m, c) = \sum_{d \mid n}^{d \le c} \sum_{e \mid n}^{e \le m/d} \mu(e) \left\lfloor \frac{ m/d }{e} \right\rfloor \]

或者

\[g(n, m, c) = \sum_{d \mid n}^{d \le c} \sum_{e \mid n}^{e \le m/d} \mu(e) \left\lfloor \frac{ m }{de} \right\rfloor \]

我们可以 \(O(\sqrt n)\) 得到 \(n\) 的所有约数 \(d\)(约 \(\sim O( \log n)\) 个),然后对于每个 \(d\) 整除分块 \(O(\sqrt n)\) 求出第二个 sigma 的值。

http://www.jsqmd.com/news/823504/

相关文章:

  • 如何在10分钟内搭建AI与Figma双向通信系统:TalkToFigma MCP完整指南
  • 植物光合作用测定仪怎么样?农业科研人员关心的实测精度与选型指南 - 品牌推荐大师1
  • 如何用嘎嘎降AI处理MBA管理论文:案例分析密集的MBA毕业论文降AI完整操作流程教程
  • Trae如何把代码上方代码文件调为多行显示
  • 脑网络分析避坑指南:GLM模型中的三种编码方式(Dummy/Effect/Cell Means)到底怎么选?附R/Python代码对比
  • ZoneMinder开源监控系统:你的专业级安防解决方案终极指南
  • 3DMAX Quad Remesher插件避坑指南:参数没调对,你的四边面拓扑等于白做
  • 国产多模态新星:mPLUG-Owl全解析,从原理到落地
  • Ketcher:三步掌握开源化学绘图工具的完整使用指南
  • 主治医师考试课程推荐|4家高口碑机构实测,在职备考也能高效通关 - 医考机构品牌测评专家
  • 为什么92%的AI团队ArgoCD部署失败?DeepSeek官方认证架构师首次公开3个被忽略的CRD权限陷阱
  • 从lspci -xxx的十六进制输出里,我们能挖出什么硬件宝藏?
  • 一站式Steam Deck控制器配置方案:Windows平台完整游戏体验指南
  • 探寻压力传感器哪家好,广东犸力以核心技术引领产业发展 - 品牌速递
  • 弹球打砖块
  • Flash Attention 核心算法与 CUDA 实现精解:从 Tiling 到 Tensor Core 优化
  • 如何在Windows平台通过用户态驱动框架实现经典游戏外设的现代化兼容?
  • 巨头转身难的地方,我们的星辰大海:开发版机巢,为千行百业而生
  • DeepSeek等低价大模型实现低算力成本的5项核心技术‌与《论三生原理》思想技术同源?
  • 【maven内网依赖缺失解决办法】
  • py每日spider案例之某百du之登录接口密码参数逆向(rsa )
  • 如何基于 Git flow 工作流管理发布分支和热修复
  • 告别网盘下载烦恼:3步解锁9大网盘高效下载新体验
  • 2026年植物冠层图像分析仪厂家怎么选?从信誉、质量到售后服务一篇文章讲清楚 - 品牌推荐大师1
  • Installing the classic Jupyter Notebook interface
  • PPTAgent:当AI成为你的演示文稿架构师
  • 别再手动数脉冲了!用STM32定时器编码器模式搞定增量编码器(附CubeMX配置)
  • 做质量工程师:日常工作的五大核心模块 - 众智商学院职业教育
  • 2026年5月物联网水肥一体化智能灌溉系统实力厂家推荐榜,瑞华电子等品牌入选 - 品牌推荐大师1
  • 2026年|AI率90%怎么办?10款主流降ai率工具深度测评推荐,帮你搞定降aigc - 降AI实验室