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

[清华集训 2014] Sum

\(\sum\limits_{d=1}^{n} (-1)^{\lfloor d\sqrt{r} \rfloor}\) 的值。

\(T \leq 10^4\) 组数据,\(n \leq 10^9\)\(r \leq 10^4\)

首先特判 \(\sqrt{r}\) 为整数的情况,若 \(\sqrt{r}\) 为偶数则答案为 \(n\),否则根据 \(n\) 的奇偶性,\(n\) 为奇数则答案为 \(-1\),否则答案为 \(0\)

考虑等式 \((-1) ^ x = 1 - 2(x \bmod 2) = 1 - 2(x - 2\lfloor \dfrac{x}{2} \rfloor) = 1 - 2x + 4\lfloor \dfrac{x}{2} \rfloor\),此时原式可以化为 \(n- 2\sum\limits_{d=1}^{n} \lfloor d\sqrt{r} \rfloor + 4\sum\limits_{d=1}^{n} \lfloor \dfrac{d\sqrt{r}}{2} \rfloor\)。令 \(f(a,b,c,n)=\sum\limits_{i=1}^{n} \lfloor \dfrac{a\sqrt{r}+b}{c} \times i \rfloor\),则原题答案为 \(n-2 \times f(1,0,1,n)+4 \times f(1,0,2,n)\)

考虑如何求出 \(f(a,b,c,n)\)。我们不妨令 \(t=\dfrac{a\sqrt{r}+b}{c}\),则 \(f(a,b,c,n)=\sum\limits_{i=1}^{n} \lfloor i \times t \rfloor\),可以用类欧解决。

\(t \geq 1\) 时,记 \(p=\lfloor t \rfloor\),则有 \(f(a,b,c,n)=\sum\limits_{i=1}^{n} \lfloor i \times (t-p) + i \times p \rfloor=\sum\limits_{i=1}^{n} \lfloor i \times \dfrac{a\sqrt{r}+b-pc}{c} \rfloor + p \times \dfrac{n(n+1)}{2}=f(a,b-pc,c,n)+p \times \dfrac{n(n+1)}{2}\)

\(t < 1\) 时,考虑直线 \(y=tx\) 不会经过任何整点,则会产生贡献的最大点为 \((n,\lfloor nt \rfloor)\)。考虑容斥,可将答案化为 \(n \times \lfloor nt \rfloor - \sum\limits_{i=1}^{\lfloor nt \rfloor} \lfloor \dfrac{i}{t} \rfloor=n \times \lfloor nt \rfloor - \sum\limits_{i=1}^{\lfloor nt \rfloor} \lfloor i \times \dfrac{c}{a\sqrt{r}+b} \rfloor=n \times \lfloor nt \rfloor - \sum\limits_{i=1}^{\lfloor nt \rfloor} \lfloor i \times \dfrac{c \times (a\sqrt{r}-b)}{(a\sqrt{r}+b) \times (a\sqrt{r}-b)} \rfloor=n \times \lfloor nt \rfloor - \sum\limits_{i=1}^{\lfloor nt \rfloor} \lfloor i \times \dfrac{ac\sqrt{r}-bc}{a^2r-b^2} \rfloor=n \times \lfloor nt \rfloor - f(ac,-bc,a^r-b^2,\lfloor nt \rfloor)\)

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

相关文章:

  • 深入解析:HiTooler File Finder: macOS上速度碾压Spotlight,媲美「Everything」的文件搜索神器
  • P13552 鱼类考古学
  • P14134 【MX-X22-T5】「TPOI-4E」Get MiN? Get MeX!
  • 20231427田泽航ipsec协议验证
  • 29232428 2025-2026-1 《网络与系统攻防技术》实验六
  • 《道德经》第三十八章 - 教程
  • 2025年必收藏的8款AI论文写作神器!助你高效搞定学术写作
  • bfs dfs板子默写 真的好怕像上次一样这种题AC不了啊
  • 贪心题目
  • 【做题记录】HZOJ 多校-数论/多校-字符串/多校-图论Ⅲ
  • 2025软件工程L班
  • 2025-11-23
  • Chainlit+LlamaIndex 多模态 RAG 开发实战7:从系统架构到功能落地,搞定 PDF/PPT/ 图片全类型文件处理 - 详解
  • 使用Ansible批量安装JDK
  • 使用OpenZeppelin编写可升级智能合约(代理) - all-in
  • 实用指南:【逻辑回归】从线性模型到逻辑回归
  • vuepress2.x支持vue2吗?
  • 贪心专题 1 做题记录
  • static 静态变量
  • 【IO多路转接】IO 多路复用之 select:从接口解析到服务器实战 - 详解
  • java sql注入的危害有哪些
  • 单片机控制继电器及其原理
  • 2025-09-10-Wed-T-Milvus
  • 【Linux】 层层递进,抽丝剥茧:调度队列、命令行参数、环境变量 - 指南
  • 字符串大小写转换
  • vitepress如何支持vue2组件
  • 2025.11.23
  • 20231427田泽航第十周预习报告
  • java linux环境变量
  • java linux服务器