K. Easy Sigma(类欧几里得)

K. Easy Sigma

∑i=1n(−1)⌊i×k⌋,(n≤109,k≤104)\sum_{i = 1} ^{n} (-1) ^{\lfloor i \times \sqrt k \rfloor}, (n \le 10 ^ 9, k \le 10 ^ 4)\\ i=1n(1)i×k,(n109,k104)
考虑(−1)x=1−2×(xmod2)=1−2(x−2×x2)=1−2x+4×⌊x2⌋(-1) ^{x} = 1 - 2 \times (x \mod 2) = 1 - 2(x - 2 \times \frac{x}{2}) = 1 - 2x + 4 \times \lfloor \frac{x}{2} \rfloor(1)x=12×(xmod2)=12(x2×2x)=12x+4×2x
∑i=1n(−1)⌊i×k⌋n−2×∑i=1n⌊i×k⌋+4×∑i=1n⌊i×k2⌋\sum_{i = 1} ^{n} (-1) ^{\lfloor i \times \sqrt{k} \rfloor}\\ n - 2 \times \sum_{i = 1} ^{n} \lfloor i \times \sqrt k \rfloor + 4 \times \sum_{i = 1} ^{n} \lfloor \frac{ i \times \sqrt k}{2} \rfloor\\ i=1n(1)i×kn2×i=1ni×k+4×i=1n2i×k
我么假设f(a,b,c,n)=∑i=1n⌊(a×r+b)×ic⌋f(a, b, c, n) = \sum\limits_{i = 1} ^{n} \lfloor \frac{\left(a \times \sqrt r + b \right) \times i}{c} \rfloorf(a,b,c,n)=i=1nc(a×r+b)×i

如果a×r+bc≥1,d=⌊a×r+bc⌋\frac{a \times \sqrt r + b} {c} \ge 1, d = \lfloor \frac{a \times \sqrt r + b} {c} \rfloorca×r+b1,d=ca×r+b
d×n×(n+1)2+∑i=1n⌊(a×r+b−d×c)×ic⌋f(a,b,c,n)=d×n×(n+1)2+f(a,b−d×c,c,n)d \times \frac{n \times (n + 1)}{2} + \sum_{i = 1} ^{n} \lfloor \frac{\left(a \times \sqrt r + b - d \times c \right) \times i}{c} \rfloor \\ f(a, b, c, n) = d \times \frac{n \times (n + 1)}{2} + f(a, b - d \times c, c, n)\\ d×2n×(n+1)+i=1nc(a×r+bd×c)×if(a,b,c,n)=d×2n×(n+1)+f(a,bd×c,c,n)

否则a×r+bc<1\frac{a \times \sqrt r + b}{c} < 1ca×r+b<1

由类欧的一般套路,我们设m=⌊n×ar+bc⌋,t=a×r+bcm = \lfloor n \times \frac{a \sqrt r + b}{c} \rfloor, t = \frac{a \times \sqrt r + b}{c}m=n×car+b,t=ca×r+b,这里的ddd没有向下取整,是一个实数。
∑i=1n∑j=1m[j<t×i]∑i=1n∑j=1m[i>jt]∑j=1mn−⌊jt⌋∑i=1mn−⌊i×ca×r+b⌋n×m−∑i=1m⌊i×ca×r+b⌋\sum_{i = 1} ^{n} \sum_{j = 1} ^{m} [j < t \times i]\\ \sum_{i = 1} ^{n} \sum_{j = 1} ^{m} [i > \frac{j}{t}]\\ \sum_{j = 1} ^{m} n - \lfloor \frac{j}{t} \rfloor \\ \sum_{i = 1} ^{m} n - \lfloor \frac{i \times c}{a \times \sqrt r + b} \rfloor\\ n \times m - \sum_{i = 1} ^{m} \lfloor \frac{i \times c}{a \times \sqrt r + b} \rfloor \\ i=1nj=1m[j<t×i]i=1nj=1m[i>tj]j=1mntji=1mna×r+bi×cn×mi=1ma×r+bi×c
考虑将分母有理化,也就是乘上一个a×r−ba \times \sqrt r - ba×rb
n×m−∑i=1m⌊car−cba2r−b2i⌋n \times m - \sum_{i = 1} ^{m} \lfloor \frac{ca \sqrt r - cb}{a ^ 2 r - b ^ 2} i \rfloor\\ n×mi=1ma2rb2carcbi
综上,可以通过类欧来解决,递归出口就是n=0n = 0n=0

#include <bits/stdc++.h>using namespace std;double x;int r;long long f(long long a, long long b, long long c, long long n) {if (!n) {return 0;}long long t = (a * x + b) / c;if (t) {return n * (n + 1) / 2 * t + f(a, b - t * c, c, n);}long long A = c * a, B = -c * b, C = a * a * r - b * b;long long d = __gcd(A, __gcd(B, C)), m = n * (a * x + b) / c;A /= d, B /= d, C /= d;return n * m - f(A, B, C, m);
}int main() {// freopen("in.txt", "r", stdin);// freopen("ou.txt", "w", stdout);int T, n;scanf("%d", &T);while (T--) {scanf("%d %d", &n, &r);x = sqrt(r);int s = sqrt(r);if (s * s == r) {if (s % 2 == 0) {printf("%d\n", n);}else if (n % 2 == 0) {puts("0");}else {puts("-1");}}else {printf("%lld\n", n - 2 * f(1, 0, 1, n) + 4 * f(1, 0, 2, n));}}return 0;
}

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

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

相关文章

Spring cloud——Hystrix 原理解析

1、背景 分布式系统环境下&#xff0c;服务间类似依赖非常常见&#xff0c;一个业务调用通常依赖多个基础服务。如下图&#xff0c;对于同步调用&#xff0c;当库存服务不可用时&#xff0c;商品服务请求线程被阻塞&#xff0c;当有大批量请求调用库存服务时&#xff0c;最终可…

【B】替换 Quartz.net 默认使用的 MySql.Data 为 Mysqlconnector 的学习过程

文章转载授权级别&#xff1a;B无论是 Quartz.net 还是 MySql.Data 都是我们比较熟悉的库了&#xff0c;Quartz.net 如果配置为使用 MySql 数据库做持久化时&#xff0c;默认是硬编码了使用 MySql.Data 来操作 MySql 数据库的。下面是我的一些个人诉求和实践&#xff0c;和大家…

G. GCD Festival(莫比乌斯、欧拉函数)

G. GCD Festival ∑i1n∑j1ngcd⁡(ai,aj)gcd⁡(i,j)∑d1nd∑i1nd∑j1ndgcd⁡(aid,ajd)[gcd⁡(i,j)1]∑d1nd∑k1ndμ(k)∑i1nkd∑j1nkdgcd⁡(aikd,ajkd)Tkd∑T1n∑i1nT∑j1nTgcd⁡(aiT,ajT)∑d∣Tdμ(Td)∑T1nϕ(T)∑i1nT∑j1nTgcd⁡(aiT,ajT)\sum_{i 1} ^{n} \sum_{j 1} ^{n}…

APM(应用性能管理)与Dapper原理介绍

一、APM&#xff08;应用性能管理&#xff09; 1.1 什么是APM&#xff1f; APM (Application Performance Management) 即应用性能管理&#xff08;应用性能监控&#xff09; APM主要是针对企业 关键业务的IT应用性能和用户体验的监测、优化&#xff0c;提高企业IT应用的可靠…

asp.netcore3.0 使用 DbProviderFactories 连接数据库

在.netstandard2.0时 System.Data.Common 这个包里并没有加入DbProviderFactoriesDbProviderFactories类在.netframework中是非常重要的存在,依靠他可以适配各种数据库客户端&#xff08;sqlserver、mysql、sqllite等&#xff09;创建数据库连接。现在可以像.netframework中一样…

MIT 6.824 Lab 1 MapReduce

MapReduce 目标 根据论文所说明的&#xff0c;有MASTER和WORKER两类工作节点&#xff0c;以下实现大都按照论文所说的实现&#xff0c;但是在对MASTER的实现上有所改动&#xff1a; MASTER向WORKER发送心跳检测&#xff0c;这里改为了对分配出去的任务进行超时监控。 MASTER…

大家在寻找的高级程序员到底是什么样子的?

你好&#xff0c;我是Z哥。这篇文章主题很简单&#xff0c;就是一个很常见的话题“什么是高级程序员&#xff1f;”。文章稍微长了些&#xff0c;但是很容易阅读。我们的中国文化&#xff0c;对“面子”看的特别重&#xff0c;所以你会发现身边到处都是高级XXX&#xff0c;听着…

应用服务器——JBoss架构分析

JBoss是什么? JBoss是免费的&#xff0c;开放源代码J2EE的实现&#xff0c;它通过LGPL许可证进行发布。它提供了基本的EJB容器以及EJB(好像应该是J2EE)服务&#xff0c;例如&#xff1a;数据库访问JDBC、交易(JTA/JTS)、消息机制(JTS)、命名机制(JNDI)和管理支持(JMX)。目前的…

优秀的程序员是那种过单行线马路都要往两边看的人

最近一周帮我以前一个同事推荐工作&#xff0c;顺便了解下行情&#xff0c;我这个同事我感觉还行&#xff0c;技术不说有多好&#xff0c;但是往年绝对不至于简历筛选时被刷掉那种&#xff0c;最先开始推给了一个我比较信任的HR手里&#xff0c;她兼职猎头&#xff0c;推给这个…

应用服务器——tomcat架构分析

先mark&#xff0c;后续补充 https://blog.csdn.net/qq_38245537/article/details/79009448

【为自己相亲】单身小姐姐你在哪里,我是书豪,我在等你

笔者简介Introduction书豪&#xff1a;【人工智能爱好者社区】公众号负责人《R数据科学实战&#xff1a;工具详解与案例分析》书籍作者。 你没看错这是书豪在给自己寻觅良缘如果你有&#xff0c;或者身边的朋友有兴趣请与我联系基本信息 出生日期&#xff1a;1995年5月身高&am…

应用服务器——jetty架构分析

先mark&#xff0c;后续补充 PS&#xff1a;网上找的资料要么太旧&#xff08;2010年的&#xff09;&#xff0c;要么乱七八糟不知所云。百度真弱鸡

知道的越多,越感觉自己渺小

作者&#xff1a;猛哥&#xff0c;关注技术和人文发展的程序员&#xff0c;架构师社区合伙人芝诺说&#xff1a;“人的知识就像一个圆&#xff0c;圆圈外是未知的&#xff0c;圆圈内是已知的&#xff0c;你知道的越多&#xff0c;你的圆圈就会越大。圆的周长也就越大&#xff0…

Java8新特性解析

Java 8中的新特性有&#xff1a; 接口中默认方法&#xff0c;lambda 表达式,方法引用,重复注解&#xff0c;流、函数、接口、map扩展、日期中的新变化等&#xff0c;接下来一一介绍 1. Default Methods for Interfaces(接口中的默认方法) Java 8准许我们在接口中增加一个通过…

.NET Core 3.0 System.Text.Json 和 Newtonsoft.Json 行为不一致问题及解决办法

行为不一致.NET Core 3.0 新出了个内置的 JSON 库, 全名叫做尼古拉斯 System.Text.Json - 性能更高占用内存更少这都不是事...对我来说, 很多或大或小的项目能少个第三方依赖项, 还能规避多个依赖项的依赖 Newtonsoft.Json 版本不一致的问题, 是件极美的事情.但是, 结果总不是不…

Java9 新特性

在介绍 java9 之前&#xff0c;我们先来看看java成立到现在的所有版本。 1990年初&#xff0c;最初被命名为Oak&#xff1b;1995年5月23日&#xff0c;Java语言诞生&#xff1b;1996年1月&#xff0c;第一个JDK-JDK1.0诞生&#xff1b;1996年4月&#xff0c;10个最主要的操作系…

深入探究Kubernetes - 初识容器

♥2019年8月28星期三第47篇原创引言最近Kubernetes比较火&#xff0c;新技术快速火起来&#xff0c;一定有它强大的优势&#xff0c;Hr反馈&#xff0c;招聘时会Kubernetes的很少&#xff0c;风口上的Kubernetes一起学学&#xff1f;扫盲贴&#xff0c;参考《Kubernetes进阶实践…

Java10 新特性

作为当今使用最广泛的编程语言之一的 Java 在 2018 年 3 月 21 日发布了第十个大版本。为了更快地迭代、更好地跟进社区反馈&#xff0c;Java 语言版本发布周期调整为每隔 6 个月发布一次。Java 10 是这一新规则之后&#xff0c;采用新发布周期的第一个大版本。Java 10 版本带来…

WinForms项目升级.Net Core 3.0之后,没有WinForm设计器?

目录.NET Conf 2019Window Forms 设计器2019 9.23-9.25召开了 .NET Conf 2019 大会&#xff0c;大会宣布了 .Net Core 3.0 正式版。这两天我也开始试着将自己Github上的项目从 .Net Core 2.2升级到 .Net Core 3.0 。其中有一个项目&#xff0c;是fork下来的Windows Forms项目&a…

Java11 新特性

Java 11新特性的详细解释。JDK 11已经于 2018年9月25日正式发布&#xff0c;那么Java 11主要包含哪些新特性呢&#xff1f; JDK 11是Java SE 11平台版本11的开源参考实现&#xff0c;由JSR 384在Java Community Process中指定。 阿里巴巴是中国唯一的JCP委员会成员公司&#x…