具体数学 - 底和顶的和式 Floor and Ceiling Sums

news/2025/2/24 14:34:50/文章来源:https://www.cnblogs.com/wenbinteng/p/18734004

通常在含有取整的和式中,求封闭形式解的行之有效的技巧是通过引入一个新的变量来规避底或顶。

例 1:求 \(\sum_{0 \leq k \leq n} \lfloor \sqrt{k} \rfloor\)

\(m=\lfloor \sqrt{k} \rfloor\),则

\[\begin{align*} \sum_{0 \leq k \leq n} \lfloor \sqrt{k} \rfloor &= \sum_{k,m \geq 0} m[k<n][m=\lfloor \sqrt{k} \rfloor] \\ &= \sum_{k,m \geq 0} m[k<n][m \leq \sqrt{k} < m+1] \\ &= \sum_{k,m \geq 0} m[k<n][m^2 \leq k < (m+1)^2] \\ &= \sum_{k,m \geq 0} m[m^2 \leq k < (m+1)^2 \leq n] + \sum_{k,m \geq 0} m[m^2 \leq k < n < (m+1)^2]. \end{align*} \]

我们假设 \(n=a^2\),则加号右侧的式子条件为空,求和结果为零。加号左侧的式子为

\[\begin{align*} \sum_{k,m \geq 0} m[m^2 \leq k < (m+1)^2 \leq n] &= \sum_{m \geq 0} m((m+1)^2-m^2)[m+1 \leq a] \\ &= \sum_{m \geq 0} m(2m+1)[m<a] \\ &= \sum_{m \geq 0} (2m^{\underline{2}}+3m^{\underline{1}})[m<a] \\ &= \sum_0^a (2m^{\underline{2}}+3m^{\underline{1}}) \delta m \\ &= \frac{2}{3} a(a-1)(a-2)+\frac{3}{2} a(a-1) \\ &= \frac{1}{6}(4 a+1) a(a-1). \end{align*} \]

一般情形下,我们假设 \(a = \lfloor \sqrt{n} \rfloor\),加号左侧的式子求解过程同上。加号右侧的式子为

\[\sum_{k,m \geq 0} m[m^2 \leq k < n < (m+1)^2] = (n-a^2)a. \]

则所求封闭形式的解为

\[\sum_{0 \leq k \leq n} \lfloor \sqrt{k} \rfloor = na-\frac{1}{3}a^3-\frac{1}{2}a^2-\frac{1}{6}a,a = \lfloor \sqrt{n} \rfloor. \]

例 2:求 \(\sum_{0 \leq k < m} \lfloor \frac{nk+x}{m} \rfloor\)

通过枚举 \(m=1,2,\dots\),可以发现和式满足形式

\[a \lfloor \frac{x}{a} \rfloor + bn + c. \]

我们对式子做以下变换

\[\lfloor \frac{nk+x}{m} \rfloor = \lfloor \frac{(nk+x) \bmod m}{m} \rfloor + \frac{kn}{m} - \frac{kn \bmod m}{m} \]

对第一项,通过枚举后发现其值是周期性的,其周期重复次数是 \(d=\gcd(m,n)\),则第一项求和的值为

\[\begin{align*} \sum_{0 \leq k < m} \lfloor \frac{(nk+x) \bmod m}{m} \rfloor &= \lfloor \frac{x \bmod m}{m} \rfloor + \lfloor \frac{(n+x) \bmod m}{m} \rfloor + \cdots + \lfloor \frac{(n(m-1)+x) \bmod m}{m} \rfloor \\ &= d(\overbrace{\lfloor \frac{x}{m} \rfloor + \lfloor \frac{(x+d)}{m} \rfloor + \cdots + \lfloor \frac{(x+m-d)}{m} \rfloor}^m) \\ &= d(\overbrace{\lfloor \frac{x/d}{m/d} \rfloor + \lfloor \frac{x/d+1}{m/d} \rfloor + \cdots + \lfloor \frac{x/d+m/d-1}{m/d} \rfloor}^m) \\ &= d \lfloor \frac{x}{d} \rfloor. \end{align*} \]

对第二项求和有

\[\sum_{0 \leq k < m}\frac{kn}{m} = \frac{(m-1)n}{2}. \]

对第三项,其值同样有周期性,其求和的值为

\[\sum_{0 \leq k < m} \frac{kn \bmod m}{m} = d(\frac{0}{m} + \frac{d}{m} + \cdots + \frac{m-d}{m}) = \frac{m-d}{2}. \]

则所求

\[\sum_{0 \leq k < m} \lfloor \frac{nk+x}{m} \rfloor = d \lfloor \frac{x}{d} \rfloor + \frac{m-1}{2}n + \frac{d-m}{2}. \]

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/889017.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

ios SDK AB 开关切换

在数据库的这个服务器 然后再ctest1数据库新建编辑器然后查询select* fromapp_config ac whereaccess_no = 12100186 //这个是应用IDand module = abSwitchand param_name = export_otel_ab查到后,把param_value改为B,或者A,然后回车,然后点击图中的保存 保存后等两分钟,…

CS Course Learning

【李宏毅】2024大语言模型课程 课程学习课程链接:https://speech.ee.ntu.edu.tw/~hylee/genai/2024-spring.php Bilibili相关视频链接:https://www.bilibili.com/video/BV1XS411w7qrGPT: Autoregressive model In-context LearningChain of Thoughts (CoT) Tree of Thoughts …

跟着狂神学markdown作业01天

markdown学习 标题 一共可以做六级标题 格式为#+空格+标题 几级标题就打几个空格 字体 粗体:hello,world 两边各加两个*号 斜体:hello,world 两边各加一个*号 粗体+斜体:hello,world 两边各加三个***号 删除效果:hello,world 引用选择狂神说java,走向人生巅峰(用>…

java知识面试day4

1.常见的关键字有哪些static:静态变量,静态变量被所有对象共享,在内存中只有一个副本。具有静态变量,静态方法块,静态代码块(在类加载时候被指执行一次),静态内部类:非静态内部类需要依赖外部实列,但静态内部类不需要。final 基本数据类型用final修饰不能修改,引用对象被…

[QOJ 8366] 火车旅行

毒瘤边化点,有人说非排列只需要加一些细节,但是这个题毒瘤在于非排列。 statement 给定一个长度为 \(n\) 的序列 \(a_i\)。 对于位置 \(x\) 和 \(y\):若 \(y < x\) 且 \(max_{y < i < x} a_i < min(a_x, a_y)\) 则位于 \(x\) 的棋子可以花费 \(L_x\) 的代价跳到…

uipath更新到最新版本2025.0.161出现严重问题

uipath更新到最新版本2025.0.161出现严重问题:1. 打开既有项目,会报CS0246错误2. 无法创建新项目,一直报无权限访问尝试办法:1. 重新安装uipath,未解决2. 删除项目重新添加,未解决3. 给账户添加最高权限,未解决 workaround:把项目从默认文件夹复制到其他盘(除了C盘外…

Python正则表达式之re.compile函数

​在Python编程语言中,re.compile函数是正则表达式模块(re)中的一个核心组件,它负责将文本形式的正则表达式编译成一个正则表达式对象。这个对象随后可以被用来执行高效的模式匹配操作,如查找、替换或者分割字符串等。理解并有效利用 re.compile对于编写高效且可维护的正则表…

Unity Addresable打包总结第二弹

前言 前文介绍了Addressables在本地打包是怎么使用,这里介绍下怎么打远程包,并且怎么做到打增量包,Lets Go! 远程包新建一个Group,将它的 BUild & Load Paths 改为Remote,并将RemoteRes资源文件夹塞入Remote Group,其中包含一个Capsule.prefab资源:在Addressabvles …

BUUCTF-RE-[2019红帽杯]easyRE

这道题很难,但是并不难在他的解题要用到的方法和技巧上,而是难在它的题目设计。做的过程中真的有一种闯关的感觉,非常有趣 首先我们通过对字符的定位我们可以来到sub_4009C6函数 __int64 sub_4009C6() {__int64 result; // raxint i; // [rsp+Ch] [rbp-114h]__int64 v2; // …

2025年免费项目管理软件哪家强?5款零成本工具实测报告

在项目管理领域,众多团队尤其是初创企业和小型项目组,都渴望找到功能实用且零成本的软件来助力项目推进。2025 年,有 5 款免费项目管理软件表现突出,它们分别是禅道、Trello、Asana、Redmine 以及国内新兴的钉钉项目管理相关功能,下面将为大家带来详细的实测报告。一、禅道…

BT151-ASEMI电机控制专用BT151

BT151-ASEMI电机控制专用BT151编辑:ll 在当今科技飞速发展的时代,电子元件宛如繁星般点缀着各个领域,而 BT151 单向可控硅无疑是其中一颗极为耀眼的明星。 BT151-ASEMI电机控制专用BT151 型号:BT151 品牌:ASEMI 封装:TO-220F 正向电流:12A 反向电压:600V~800V 引脚数量…

销售必备60招

脸皮厚、嘴巴甜、手脚勤。先放下面子努力拿下客户,提高收入,往后在亲友面前抢着买单时,面子自然就有了。干销售得有狼性、野性和匪性,还要有胆量。推销是常识的运用,但只有将经实践验证的理念用于积极的客户身上,才会产生效果。想要取得惊人成绩,必须先做好枯燥乏味的准…