C - Maximize GCD(简单数论)

C - Maximize GCD

给定长度为n,(2≤3×105)n, (2 \le 3 \times 10 ^ 5)n,(23×105)的数组a,(1≤ai≤3×105)a, (1 \le a_i \le 3 \times 10 ^ 5)a,(1ai3×105),一个数字K,(1≤K≤1018)K, (1 \le K \le 10 ^{18})K,(1K1018)

我们可以对数组aaa进行最多kkk次操作,每次操作选定一个i,(1≤i≤n)i, (1 \le i \le n)i,(1in)使ai+1a_i + 1ai+1,为数组最大的可能gcdgcdgcd是多少。

我们设maxn=max({a1,…,an})maxn = max(\{a_1, \dots, a_n\})maxn=max({a1,,an}),如果最后的答案ansansans大于等于maxnmaxnmaxn,则一定有i∈[1,n],ai=ansi \in [1, n], a_i = ansi[1,n],ai=ans

否则我们假设答案为xxx,则判断是否合法:
∑i=1n⌈aix⌉×x−ai\sum_{i = 1} ^{n} \lceil\frac{a_i}{x}\rceil \times x - a_i\\ i=1nxai×xai
即上面那个式子求和是否小于等于kkk,这个式子可以前缀和,无穷级数的复杂度算一下就好了,整体复杂度最高maxnlog⁡maxnmaxn \log maxnmaxnlogmaxn

#include <bits/stdc++.h>using namespace std;const int N = 3e5 + 10;int a[N], cnt[N], n, maxn;long long sum[N], k, s;int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);cin >> n >> k;for (int i = 1; i <= n; i++) {cin >> a[i];s += a[i], cnt[a[i]]++, sum[a[i]] += a[i];maxn = max(maxn, a[i]);}if (1ll * maxn * n - s <= k) {// x * n - sum <= k -> x * n <= k + sum -> x = (k + sum) / n;cout << (k + s) / n << "\n";return 0;}for (int i = 1; i <= maxn; i++) {cnt[i] += cnt[i - 1];sum[i] += sum[i - 1];}int ans = 1;for (int i = 2; i < maxn; i++) {long long cur = 0;for (int j = i; j <= maxn; j += i) {// [j - i + 1, j] -> icur += 1ll * (cnt[j] - cnt[j - i]) * j - (sum[j] - sum[j - i]);}if (maxn % i) {int l = maxn / i * i, r = maxn;cur += 1ll * (cnt[r] - cnt[l]) * (maxn / i + 1) * i - (sum[r] - sum[l]);}if (cur <= k) {ans = i;}}cout << ans << "\n";return 0;
}

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

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

相关文章

推荐.neter常用优秀开源项目系列之二

.net社区有很多优秀的开源项目&#xff0c;我们今天再推荐12个开源项目&#xff1b;1. Domain-Driven-Design-ExampleDDD 示例 挺不错的。github https://github.com/zkavtaskin/Domain-Driven-Design-Example2.SmartStoreNET开源的电商项目github https://github.com/smartsto…

Zookeeper: Zookeeper架构及FastLeaderElection机制

本文转发自技术世界&#xff0c;原文链接 http://www.jasongj.com/zookeeper/fastleaderelection/ 一、Zookeeper是什么 Zookeeper是一个分布式协调服务&#xff0c;可用于服务发现&#xff0c;分布式锁&#xff0c;分布式领导选举&#xff0c;配置管理等。 这一切的基础&am…

I. Rise of Shadows(类欧几里得)

I. Rise of Shadows 一天有HHH个小时&#xff0c;MMM分钟&#xff0c;问&#xff0c;有多少个整数分钟&#xff0c;满足时针与分针的角度≤α\le \alpha≤α&#xff0c;α2πAHM\alpha \frac{2 \pi A}{HM}αHM2πA​。 ∑i0H−1∑j0M−1[∣2π(iMj)HM−2πjM∣≤2πAHM]∑i0…

与时俱进 | 博客现已运行在 .NET Core 3.0 及 Azure 上

点击上方蓝字关注“汪宇杰博客”导语9月23日&#xff0c;微软正式发布了 .NET Core 3.0&#xff0c;这个版本具有大量新功能和改进。我也在第一时间将自己的博客网站更新到了 .NET Core 3.0&#xff0c;并且仍然跑在微软智慧云 Azure 国际版的应用服务上。本文总结了我在博客迁…

Zookeeper:基于Zookeeper的分布式锁与领导选举

本文转发自技术世界&#xff0c;原文链接 http://www.jasongj.com/zookeeper/distributedlock/ 1、Zookeeper特点 1.1 Zookeeper节点类型 如上文《Zookeeper架构及FastLeaderElection机制》所述&#xff0c;Zookeeper 提供了一个类似于 Linux 文件系统的树形结构。该树形结构…

Asp.Net Core Mvc Razor之RazorPage

在AspNetCore.Mvc.Razor命名空间中的RazorPage继承RazorPageBase&#xff0c;并定义的属性为&#xff1a;HttpContext Context 表示当前请求执行的HttpContextRazorPageBase定义为抽象类&#xff0c;并继承了接口&#xff1a;IRazorPageIRazorPage接口定义属性如下&#xff1a;…

#138. 类欧几里得算法

#138. 类欧几里得算法 以下除法均为向下取整&#xff0c; 定义f(a,b,c,n,k1,k2)∑x0nxk1(axbc)k2f(a, b, c, n, k_1, k_2) \sum\limits_{x 0} ^{n} x ^{k_1} \left(\frac{a \times x b}{c}\right) ^ {k_2}f(a,b,c,n,k1​,k2​)x0∑n​xk1​(caxb​)k2​。 ∑x0nxk1(axbc)k2…

Spring Cloud——Consul——架构体系

我们知道&#xff0c;Eureka 2.X因遇到问题&#xff0c;已停止研发。Spring Cloud官方建议迁移到Consul或者Zookeeper等其他服务发现中间件。 下面是 Spring Cloud 支持的服务发现软件以及特性对比&#xff1a; 一、Consul 介绍 Consul 是 HashiCorp 公司推出的开源工具&…

ASP.NET Core 3.0 gRPC 双向流

目录ASP.NET Core 3.0 使用gRPCASP.NET Core 3.0 gRPC 双向流ASP.NET Core 3.0 gRPC 认证授权一.前言在前一文 《二. 什么是 gRPC 流gRPC 有四种服务类型&#xff0c;分别是&#xff1a;简单 RPC&#xff08;Unary RPC&#xff09;、服务端流式 RPC &#xff08;Server streami…

Spring cloud Loadbalance

feign原理介绍&#xff1a; https://segmentfault.com/a/1190000016776471?utm_sourcetag-newest http://techblog.ppdai.com/2018/05/28/20180528/ https://blog.csdn.net/xiaofeng10330111/article/details/85682513 负载均衡策略&#xff1a; 随机 (Random) 随机策略很…

类欧几里得(模板题推导)

类欧几里得 设三个函数f(a,b,c,n)∑i0naibc,g(a,b,c,n)∑i0niaibc,h(a,b,c,n)∑i0n(aibc)2f(a, b, c, n) \sum\limits_{i 0} ^{n} \frac{a \times i b}{c}, g(a, b, c, n) \sum\limits_{i 0} ^{n} i \times \frac{a \times i b}{c}, h(a, b, c, n) \sum\limits_{i 0} …

开源公司被云厂商“寄生”,咋整?

上周 OSS Capital 召集一些开源公司&#xff0c;组织了一场关于如何面对“云厂商给开源带来的危害”的会议。OSS Capital 是一家风险投资公司&#xff0c;该公司只投开源&#xff0c;其董事会合伙人之一是开源运动的先驱人物 Bruce Perens。网上有一个十分有名的“开源商业化独…

Spring Cloud Config——原理解析

springCloud config项目,用来为分布式的微服务系统中提供集成式外部配置支持,分为客户端和服务端 可以让你把配置放到远程服务器&#xff0c;目前支持本地存储、Git以及Subversion。 spring官方如下介绍: 简而言之: 通过配置服务(Config Server)来为所有的环境和应用提供外部配…

AWS加入.NET Foundation企业赞助商计划

.NET 走向开源&#xff0c;MIT许可协议。 微软为了推动.NET开源社区的发展&#xff0c;2014年联合社区成立了.NET基金会。.NET基金会是一个独立的组织&#xff0c;支持.NET社区和开源&#xff0c;旨在拓宽和加强.NET生态系统和社区。这可以通过多种方式完成&#xff0c;包括项目…

K. Easy Sigma(类欧几里得)

K. Easy Sigma ∑i1n(−1)⌊ik⌋,(n≤109,k≤104)\sum_{i 1} ^{n} (-1) ^{\lfloor i \times \sqrt k \rfloor}, (n \le 10 ^ 9, k \le 10 ^ 4)\\ i1∑n​(−1)⌊ik​⌋,(n≤109,k≤104) 考虑(−1)x1−2(xmod2)1−2(x−2x2)1−2x4⌊x2⌋(-1) ^{x} 1 - 2 \times (x \mod 2) 1 …

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中一样…