NC15052 求最值

题目

题目描述

给你一个长为 \(n\) 的序列 \(a\)

定义 \(f(i,j)=(i-j)^2+g(i,j)^2\)

\(g\) 是这样的一个函数

求最小的 \(f(i,j)\) 的值,\(i!=j\)

输入描述

第一行一个数 \(n\)

之后一行 \(n\) 个数表示序列 \(a\)

输出描述

输出一行一个数表示答案

示例1

输入

4
1 0 0 -1

输出

1

备注

对于 \(100\%\) 的数据,\(2 \leq n \leq 100000 , |ai| \leq 10000\)

题解

思路

方法一

知识点:分治,排序。

观察到, \(g(i,j) = (i-j)^2 + (sum(j)-sum(i))^2\) ,其中 \(sum(i)\) 指 \([1,i]\) 的前缀和。

将其转换为平面上的点 \((i,sum(i)),(j,sum(j))\) 则函数 \(g\) 为两点之间距离的平方,问题转换为平面最近点对问题。

这里一开始处理完已经是按 \(x\) 排序好的,直接开始分治。设区间左右端点为 \(l,r\) ,若 \(l == r\) ,则直接返回无穷大;若 \(l+1 == r\) ,则返回两点距离;若其他情况则先求得对半左右两个区间的临时最小值 \(ans\) ,用这个 \(ans\) 作为阈值找到 \(x\) 轴坐标与区间分界坐标 \(mid\) 差的平方小于 \(ans\) 的点对,然后用朴素算法找到他们之间的最小值,过程中不断更新 \(ans\) 。在最后一步,有个关键优化,先按 \(y\) 轴排序,这样对于一个定点找另一点时,若 \(y\) 轴差值平方大于等于目前 \(ans\) 可以直接跳过,可以证明大部分情况可以将 \(O(n^2)\) 降到 \(O(n)\) 。

时间复杂度 平均: \(O(n\log^2n)\) 最差:\(O(n^2 \log n)\)

空间复杂度 \(O(n)\)

方法二

知识点:分治,排序。

方法二实际上是方法一的优化。通过 \([l,mid]\) 和 \([mid+1,r]\) 的 \(d^2\) 最小值,返回合并区间的 \([l,r]\) 中 \(d^2\) 最小值,在合并时用归并把当前区间的区间按 \(y\) 排序,方便查找合并区间的最小值,是 \(O(n\log n)\) 。如果每次不排序,直接取点 \(sort\) 按 \(y\) 排序,会重复排序,最差会 \(O(n^2)\) 。

要注意的是,如果采用这种方法我们就必须是先排序好再确定可行下标,但是确定下标需要原来 \(mid\) 处的 \(x\) 坐标,因此需要额外设置变量保存。

时间复杂度 \(O(n\log^2n)\)

空间复杂度 \(O(n)\)

代码

方法一

#include <bits/stdc++.h>
#define ll long long using namespace std; struct Point {
ll x, y;
}a[100007];
int c[100007]; ll sqr(ll x) { return x * x; } ll dist2(Point a, Point b) { return sqr(a.x - b.x) + sqr(a.y - b.y); } ll solve(int l, int r) {
if (l == r) return ~(1LL << 63);
if (l + 1 == r) return dist2(a[l], a[r]);
int mid = l + r >> 1;
ll ans = min(solve(l, mid), solve(mid + 1, r));
int cnt = 0;
for (int i = l;i <= r;i++) if (sqr(a[i].x - a[mid].x) < ans) c[cnt++] = i;
sort(c, c + cnt, [&](int p, int q) {return a[p].y < a[q].y;});
for (int i = 0;i < cnt;i++) {
for (int j = i + 1;j < cnt;j++) {
if (sqr(a[c[i]].y - a[c[j]].y) >= ans) break;
ans = min(ans, dist2(a[i], a[j]));
}
}
return ans;
} int main() {
std::ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
int n;
cin >> n;
for (int i = 1;i <= n;i++) a[i].x = i, cin >> a[i].y, a[i].y += a[i - 1].y;
cout << solve(1, n) << '\n';
return 0;
}

方法二

#include <bits/stdc++.h>
#define ll long long using namespace std; struct Point {
ll x, y;
}a[100007], b[100007];
ll c[100007]; ll sqr(ll x) { return x * x; } ll dist2(Point a, Point b) { return sqr(a.x - b.x) + sqr(a.y - b.y); } ll solve(int l, int r) {
if (l == r) return ~(1LL << 63);
if (l + 1 == r) return dist2(a[l], a[r]);
int mid = l + r >> 1;
ll midx = a[mid].x;
ll ans = min(solve(l, mid), solve(mid + 1, r));
merge(a + l, a + mid + 1, a + mid + 1, a + r + 1, b + l, [&](Point a, Point b) {return a.y < b.y;});
for (int i = l;i <= r;i++) a[i] = b[i];
int cnt = 0;
for (int i = l;i <= r;i++) if (sqr(a[i].x - midx) < ans) c[cnt++] = i;
for (int i = 0;i < cnt;i++) {
for (int j = i + 1;j < cnt;j++) {
if (sqr(a[i].y - a[j].y) >= ans) break;
ans = min(ans, dist2(a[i], a[j]));
}
}
return ans;
} int main() {
std::ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
int n;
cin >> n;
for (int i = 1;i <= n;i++) a[i].x = i, cin >> a[i].y, a[i].y += a[i - 1].y;
cout << solve(1, n) << '\n';
return 0;
}

NC15052 求最值的更多相关文章

  1. ACM3 求最值

    /*2*2014.11.18*求最值*描述:给定N个整数(1<=N<=100),求出这N个数中的最大值,最小值.*输入:多组数据,第一行为一个整数N,第二行为N个不超过100的正整数,用空 ...

  2. [NOI2005]维修数列 Splay tree 区间反转,修改,求和,求最值

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1500 Description Input 输入文件的第1行包含两个数N和M,N表示初始时数 ...

  3. hdu4521-小明系列问题——小明序列(线段树区间求最值)

    题意:求最长上升序列的长度(LIS),但是要求相邻的两个数距离至少为d,数据范围较大,普通dp肯定TLE.线段树搞之就可以了,或者优化后的nlogn的dp. 代码为  线段树解法. #include ...

  4. javascript之求最值

    求最值: var selections = $("#deliveryGridSalesOrGoods").datagrid('getRows'); var costPrice = ...

  5. poj3264(线段树区间求最值)

    题目连接:http://poj.org/problem?id=3264 题意:给定Q(1<=Q<=200000)个数A1,A2,```,AQ,多次求任一区间Ai-Aj中最大数和最小数的差. ...

  6. Sql示例说明如何分组后求中间值--【叶子】

    原文:Sql示例说明如何分组后求中间值--[叶子] 这里所谓的分组后求中间值是个什么概念呢? 我举个例子来说明一下: 假设我们现在有下面这样一个表: type        name price -- ...

  7. hdu 1754 I Hate It(树状数组区间求最值)2007省赛集训队练习赛(6)_linle专场

    题意: 输入一行数字,查询第i个数到第j个数之间的最大值.可以修改其中的某个数的值. 输入: 包含多组输入数据. 每组输入首行两个整数n,m.表示共有n个数,m次操作. 接下来一行包含n个整数. 接下 ...

  8. C语言 · 求arccos值

    算法提高 7-2求arccos值   时间限制:10.0s   内存限制:256.0MB      问题描述 利用标准库中的cos(x)和fabs(x)函数实现arccos(x)函数,x取值范围是[- ...

  9. 位运算求最值 学习笔记 (待补充QAQ)

    没有什么前言?直接进入正题qwq 俩俩异或 求最值: 建trie树 O(n)枚举每个数找这个数的最值,每次反走就成,还可以剪枝一波(如果在某位已经小于ans显然可以直接return? void Ins ...

随机推荐

  1. 共读《redis设计与实现》-数据结构篇

    准备将之前攒下的书先看一遍,主要是有个大概的了解,以后用的时候也知道在哪里找.所以准备开几篇共读的帖子,激励自己多看一些书. Redis 基于 简单动态字符串(SDS).双端链表.字典.压缩列表.整数 ...

  2. MKL与VS2019配置方法

    VS2019配置oneAPI并调用MKL库 oneAPI oneAPI是一个跨架构的编程工具,旨在简化跨GPU.CPU.FPGA和AI加速器之间的编程,可以与英特尔自身设备,或其他厂商的芯片配合使用, ...

  3. 解决go-micro与其它gRPC框架之间的通信问题

    在之前的文章中分别介绍了使用gRPC官方插件和go-micro插件开发gRPC应用程序的方式,都能正常走通.不过当两者混合使用的时候,互相访问就成了问题.比如使用go-micro插件生成的gRPC客户 ...

  4. 关于fiddler抓包一键生成python脚本

    本人贡献一篇关于抓包转换成脚本的文章 步骤一 打开fiddler,抓到包之后,保存成txt文件 步骤二 脚本里str_filename改成保存的文件名 步骤三 执行脚本一键转换 附上脚本,感谢关注~ ...

  5. Linux和kali Linux 介绍

    常用的渗透测试平台 CTFTools kali (近亲 Ubuntu) Parrot Security OS PentestBox --由印度人开发,运行在Windows下的渗透测试环境 kali L ...

  6. shell脚本实现MySQL全量备份+异地备份

    一.知识储备工作: Mysql导出数据库语法: mysqldump -u用户名 -p密码 数据库名 > 数据库名.sql shell脚本for循环及if条件判断基本语法 gzip压缩文件用法 r ...

  7. golang get process name by pid

    一个很好的问题:How golang to get process name by process id (pid)? 目前看来go api并没有提供通过pid获取进程名称的方法,可以通过 /proc ...

  8. Spring 源码(10)Spring Bean 的创建过程(1)

    Spring Bean的创建刚开始进行了一些准备工作,比如转换服务的初始化,占位符解析器的初始化,BeanDefinition元数据的冻结等操作,都是为了在创建Bean的过程中保证Bean的正确的创建 ...

  9. zabbix5.0报错PHP时区未设置(配置参数"date.timezone")

    解决办法 : #1.编辑文件/etc/opt/rh/rh-php72/php-fpm.d/zabbix.conf,取消注释并设置为所在地时区 vim /etc/opt/rh/rh-php72/php- ...

  10. vue中blob文件下载及其它下载方式

    一.Blob对象的了解 1:blob表示一个不可变.原始数据的类文件对象.Blob()构造函数返回一个新的blob对象:blob对象的内容由参数给出的值串联组成: 2:new Blob(array, ...