题解:[COCI2011-2012#5] BLOKOVI

Description

PDF : https://hsin.hr/coci/archive/2011_2012/contest5_tasks.pdf

题意:

有 \(N\) 个已知质量的矩形,长都是 2 ,高都是 \(h\) 。把它们放到一个平面直角坐标系里,满足:

  1. 各矩形的边与坐标轴平行;
  2. 各矩形下底边的纵坐标各不相同,分别为 \(0,h,2h,\dots,(N-1)h\) ;
  3. 最底端的矩形左下角坐标恒为 \((-2,0)\) ,右下角为原点 \((0,0)\) 。

我们记一个矩形的 \(Xcentre\) 为其下底边的中点,

一个或多个矩形的 \(Xbarycentre\) 是这些矩形的 \(Xcentre\) 的加权平均数,公式为:

\[Xbarycentre = \frac{\sum_i m_i \cdot Xcentre(i)}{\sum_im_i}
\]

其中 \(m_i\) 是各个矩形的质量。

如果 \(N\) 的矩形的摆放方案是稳定的,它应当满足对于任意的矩形 \(A\) ,在 \(A\) 上方的各矩形的 \(X-barycentre\) 符合 \(|Xbarycentre - Xcentre(A)| \leq 1\) 。

左图的摆放方案就是不稳定的,因为最上方两个矩形的 \(Xcentre\) 距离落在的下面那个矩形的右侧。

右图的的拜访方案则是稳定的。

现在我们给出从下至上各个矩形的质量,在保证摆放方案稳定的情况下,请你找出所有矩形右下角顶点的横坐标最大值可能是多少。

注意,你不能改变矩形的摆放顺序和第一个矩形的位置。

\(2\leq N \leq 3e5\)

至今,我仍然无法用自然流畅的语言翻译英文题面……大家姑且一看,建议还是读原题。

Algorithm

太神了,我只能说不愧是克罗地亚题。

首先, \(N\leq3e5\) 的规模要求了一个线性的算法,这似乎是比较困难的。

为什么呢?我们自底而上地考虑,考虑放置某个矩形的时候,未来放置的其它矩形很可能会影响它的放置。

换言之,我们需要「回头」。

说到「回头」,我们很容易想到动态规划中「后效性」的概念。

而这道题正是一个 多阶段决策最优化问题 。我们不妨就试着消除后效性,使用动态规划来解决它。

试着消除本题的后效性是一件并不困难的问题。只需要回顾刚刚我们需要「回头」的情形,很容易就能发现:

任何一个矩形的放置方法只可能影响在它之下的矩形放置。

当然上面矩形的绝对坐标也会受到影响,但是相对位置并不会发生变化。

因此,我们当然可以自顶向下地考虑各个矩形的放置,这样就自然可以消除后效性了。

当然,题目还要求最底下的那个矩形的位置固定,因此我们需要将自顶向下考虑得到的相对坐标转化为答案要求的绝对坐标。

既然后效性没了,那么接下来,我们考虑动态规划的状态转移。

我们在这里仍然存在一个小小的问题:

我们是在一个笛卡尔坐标系上放矩形,可以放置的位置可是有不可数无穷多个啊。

直觉敏感的同学或者勇敢莽撞的同学一定发现了:我们只要往最左边或者最右边放就行了。

这点是容易感性地证明的:

考虑不是把某个矩形放在最左或最右的情况,将其与其下的矩形在不改变其相对位置的情况下、整体地左移或右移,一定可以使得答案不变得更劣。

这句话稍稍有点复杂,不过其意涵是比较容易理解的。

每个矩形能放置的极限位置也是容易计算的。我们先计算向右放最远能放多少。

\[\begin{align*}
&d_k = Xcentre(k) - Xbarycentre(k + 1\sim n) \leq 1
\end{align*}
\]

由定义又有:

\[Xbarycentre(k\sim n) =\frac{ (\sum_{i = k + 1}^n m _i) \cdot Xbarycentre(k+1\sim n) + m_k \cdot Xcenter(k)}{\sum_{i = k}^n m_i}
\]

其中,

\[Xcentre(k) = d_k + Xbarycentre(k + 1 \sim n)
\]

然后一个带入,

\[Xbarycentre(k\sim n) = Xbarycentre(k+1 \sim n) + d_k \cdot \frac{m_k}{\sum_{i=k}^n m_i}
\]

按照上述分析,直接取 \(d_k = 1\) ,那么放置导致的重心向右偏移量即为 \(\Large \frac{m_k}{\sum_{i = k}^n m_i}\) 了。

由对称性,向左放最大偏移量即为 \(2 - \Large \frac{m_k}{\sum_{i = k}^n m_i}\) 。

接下来我们只要顺次考虑某个矩形是放在最左还是最右即可了,这是一个类似 01背包的模型。

代码很简单:

#include<bits/stdc++.h>
using namespace std; template<class T>
inline void read(T &x)
{
char c = getchar(); x = 0;
while(c < '0' || '9' < c) c = getchar();
while('0' <= c && c <= '9')
{
x = (x << 1) + (x << 3) + c - 48;
c = getchar();
}
} const int N = 3e5 + 10;
int n, a[N]; int main()
{
read(n);
for(int i = 1; i <= n; read(a[i++])); double ans = 0, sum = 0;
for(int i = n; i > 1; --i)
{
sum += a[i];
double del = a[i] / sum;
ans = max(ans, max(ans + del, 2 - del));
} printf("%.8f\n", ans); return 0;
}

题解:[COCI2011-2012#5] BLOKOVI的更多相关文章

  1. NOIP 2012 Day2T2 借教室题解

    NOIP 2012 Day2T2 借教室题解 题目传送门:http://codevs.cn/problem/1217/ 题目描述 Description 在大学期间,经常需要租借教室.大到院系举办活动 ...

  2. 【题解】Ples [COCI2011]

    [题解]Ples [COCI2011] 依旧是没有传送门,只有提供了数据的官网. [题目描述] \(N\) 个汉子和 \(N\) 个妹纸一起参加舞会,跳舞时只能是一个汉子一个妹纸配对,现在给出每个人的 ...

  3. 【题解】Dvoniz [COCI2011]

    [题解]Dvoniz [COCI2011] 没有传送门,只有提供了数据的官网. [题目描述] 对于一个长度为 \(2*K\) 的序列,如果它的前 \(K\) 个元素之和小于等于 \(S\) 且后 \( ...

  4. 网易游戏QA工程师笔试回忆-2012.9【个人题解】

    ========================转帖======================== 网易游戏QA工程师笔试回忆-2012.9 刚刚从武大回来,趁热回忆下题目,给以后的XDJMs参考. ...

  5. 【2012天津区域赛】部分题解 hdu4431—4441

    1001: 题意:给你13张麻将牌,问可以胡哪些张 思路: 枚举可能接到的牌,然后dfs判断能否胡 1002: 题意: 已知n,m 求 n的所有约数在m进制下的平方和 做法:队长用java高精度写的 ...

  6. 【2012长春区域赛】部分题解 hdu4420—4430

    这场比赛特点在于两个简单题太坑,严重影响了心情..导致最后只做出两题....当然也反映出心理素质的重要性 1002: 题意:一个矩阵b[n][n]通过数组 a[n]由以下规则构成,现在已知b[n][n ...

  7. ACM-ICPC Dhaka Regional 2012 题解

    B: Uva: 12582 - Wedding of Sultan 给定一个字符串(仅由大写字母构成)一个字母表示一个地点,经过这个点或离开这个点都输出这个地点的字母) 问: 每一个地点经过的次数(维 ...

  8. 【USACO 2012 Open】奶牛赛跑_题解

    奶牛赛跑 目录 奶牛赛跑 题目描述 输入格式 输出格式 样例 样例输入#1 样例输出#1 题解 代码 题目描述 约翰有头奶牛,他为这些奶牛准备了一个周长为的环形跑牛场.所有奶牛从起点同时起跑,奶牛在比 ...

  9. Luogu P1082 同余方程(NOIP 2012) 题解报告

    题目传送门 [题目大意] 求关于x的同余方程 ax≡1(mod b)的最小整数解. [思路分析] 由同余方程的有关知识可得,ax≡1(mod b)可以化为ax+by=1,此方程有解当且仅当gcd(a, ...

随机推荐

  1. 索引对单列极值查询的显著性影响(百万级别表单列最值查询 Cost由1405变成3)

    表结构: create table hy_emp( id number(7,0) primary key, name nvarchar2(20) not null, salary number(5,0 ...

  2. Centos7安装后进不去,死活就要填licence,该怎么办?

    遇到这个问题不要麻爪,跟着我做: 1 回车 2 回车 c 回车 c 回车 然后就进入系统了. 要使它联网,点右上角的开关按钮,将PCI Ethernet选择为connect状态. 我的centos7是 ...

  3. Promise核心实现

    核心 构造函数核心 维护状态变量,只能由pending变为resolve或者reject 维护一个存储结果的变量 维护一个回调数组,执行到then,如果我们传入的立即执行函数没有立即执行resolve ...

  4. 软件测试----xml文件介绍

    软件测试 目录 软件测试 一.什么是XML?: 二.XML和HTML的差异: 三.XML的特点 1.XML可以自定义标签 2.XML必须包含根元素 如上所示, 3.XML标签对大小写敏感 4.XML ...

  5. Jmeter测试工具

    jmeter的简单应用 目录 jmeter的简单应用 1.Jmeter 的基本概念 2.我们 为什么 使用 Jmeter 3.Jmeter的作用 4.Jmeter怎么用 5.安装JAVA环境 6.Jm ...

  6. 关于px、pt、em、rem四个单位的解释

    写在前面 最近在群里突然看到一个问题,就是px pt em rem 三者的区别,这个问题看起来非常基础,也非常容易被忽略,however,面试会问到~,那我就解释一下 px px的英文是pixel,翻 ...

  7. 【Flutter 实战】文件系统目录

    老孟导读:Flutter 中获取文件路径,我们都知道使用 path_provider,但对其目录对含义不是很清楚,此文介绍 Android.iOS 系统的文件目录,不同场景下建议使用的目录. 不同的平 ...

  8. Linux实战(8):查看文件大小

    du -sh *

  9. leetcode1546题解【前缀和+贪心】

    leetcode1546.和为目标值的最大数目不重叠非空子数组数目 题目链接 算法 前缀和+贪心 时间复杂度O(n). 1.对nums数组求前缀和: 2.在求前缀和过程中将前缀和sum插入到set集合 ...

  10. golang 协程学习

    协程数据传递问题 func TestGoroutineData(t *testing.T) { var wg sync.WaitGroup wg.Add(1) i := 0 go func(j int ...