BZOJ4570: [Scoi2016]妖怪
4570: [Scoi2016]妖怪
Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 491 Solved: 125
[Submit][Status][Discuss]
Description
邱老师是妖怪爱好者,他有n只妖怪,每只妖怪有攻击力atk和防御力dnf两种属性。邱老师立志成为妖怪大师,于
是他从真新镇出发,踏上未知的旅途,见识不同的风景。环境对妖怪的战斗力有很大影响,在某种环境中,妖怪可
以降低自己k×a点攻击力,提升k×b点防御力或者,提升自己k×a点攻击力,降低k×b点防御力,a,b属于正实数
,k为任意实数,但是atk和dnf必须始终非负。妖怪在环境(a,b)中的战斗力为妖怪在该种环境中能达到的最大攻击
力和最大防御力之和。strength(a,b)=max(atk(a,b))+max(dnf(a,b))环境由a,b两个参数定义,a,b的含义见前
文描述。比如当前环境a=3,b=2,那么攻击力为6,防御力为2的妖怪,能达到的最大攻击力为9,最大防御力为6。
所以该妖怪在a=3,b=2的环境下战斗力为15。因此,在不同的环境,战斗力最强的妖怪可能发生变化。作为一名优
秀的妖怪训练师,邱老师想发掘每一只妖怪的最大潜力,他想知道在最为不利的情况下,他的n只妖怪能够达到的
最强战斗力值,即存在一组正实数(a,b)使得n只妖怪在该环境下最强战斗力最低。
Input
第一行一个n,表示有n只妖怪。接下来n行,每行两个整数atk和dnf,表示妖怪的攻击力和防御力。
1≤n≤10^6, 0<atk,dnf≤10^8
Output
输出在最不利情况下最强妖怪的战斗力值,保留4位小数。
Sample Input
3
1 1
1 2
2 2
Sample Output
8.0000
题解
首先我们可以得到 $strength= atk + a / b * dnf + dnf + b / a * atk $ ,进而得到 \(strength=(a+b)/a*atk+(a+b)/b * dnf\)
所以我们将atk看为x,dnf看为y,则strength就是一条斜率为$ -b/a $ 的直线,而此时最大战力就是就是最外面的那条直线。所以我们可以将这些直线跑一个半凸包。其实按照我的想法就该三分了,然出题人告诉我这样会挂。。。
正解是这些点些点对应的斜率都有一个范围,然后将 \(strength\) 变为一个双钩函数,之后讨论最值,跟新ans。
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long LL;
LL Ina; bool InSign; char Inc;
inline LL geti() {
InSign = false;
while ((Inc = getchar()) < '0' || Inc > '9') InSign |= Inc == '-';
Ina = Inc - '0';
while ((Inc = getchar()) >= '0' && Inc <= '9') Ina = (Ina << 3) + (Ina << 1) + Inc - '0';
return InSign ? -Ina : Ina;
}
const int N = 1e6 + 5;
const double OO = 1e9;
struct P{
LL x, y;
P(LL a = 0, LL b = 0) : x(a), y(b) {}
P operator - (const P &a) const { return P(x - a.x, y - a.y); }
LL operator * (const P &a) const { return x * a.y - y * a.x; }
bool operator < (const P &a) const { return (x < a.x) || (x == a.x && y < a.y); }
void Read() { x = geti(), y = geti(); }
}p[N], Stack[N];
//k = -b/a
inline double get(const P &a) { return -sqrt((double)a.y / a.x); }
inline double getk(const P &a, const P &b) {
return (a.x ^ b.x) ? ((double)(a.y - b.y) / (double)(a.x - b.x)) : OO;
}
inline double Cal(const P &a, const double &k) {
return k >= 0 ? OO : (double)a.x + a.y + - k * a.x - a.y / k;
}
template <class T>
inline void SelfMin(T &a, const T &b) { if (b < a) a = b; }
int main() {
int n = geti(), i, j, top = 0;
for (i = 1; i <= n; ++i) p[i].Read();
sort(p + 1, p + n + 1);
Stack[++top] = p[1];
for (i = 2; i <= n; ++i) {
while (top >= 2 && (Stack[top]-Stack[top-1])*(p[i]-Stack[top-1]) >= 0)
--top;
Stack[++top] = p[i];
}
if (top < 2) return printf("%.4lf\n", Cal(Stack[1], get(Stack[1]))), 0;
double k, k1, k2, ans = OO;
k2 = getk(Stack[1], Stack[2]); k = get(Stack[1]);
if (k >= k2) SelfMin(ans, Cal(Stack[1], k));
k1 = getk(Stack[top - 1], Stack[top]); k = get(Stack[top]);
if (k <= k1) SelfMin(ans, Cal(Stack[top], k));
SelfMin(ans, Cal(Stack[top], k1));
for (i = 2; i < top; ++i) {
k1 = getk(Stack[i - 1], Stack[i]), k2 = getk(Stack[i], Stack[i + 1]);
k = get(Stack[i]);
SelfMin(ans, Cal(Stack[i], k1));
if ((k <= k1) && (k >= k2)) SelfMin(ans, Cal(Stack[i], k));
}
return printf("%.4lf\n", ans), 0;
}
BZOJ4570: [Scoi2016]妖怪的更多相关文章
- [Bzoj4570][Scoi2016]妖怪(右上凸包)
4570: [Scoi2016]妖怪 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 1110 Solved: 336[Submit][Status][ ...
- [BZOJ4570][SCOI2016]妖怪(凸包)
两种做法,前一种会TLE. 第一种是高一数学题做法,设一个妖怪的atk和dnf分别为x和y,则它在(a,b)环境下的战斗力为x+y/a*b+y+x/a*b. 设t为b/a,则战斗力即$f(x,y,t) ...
- BZOJ4570 SCOI2016妖怪(三分)
strength=atk*(1+b/a)+dnf*(1+a/b).设a/b=x,可以发现这是一个关于x的对勾函数.开口向上的一堆凸函数取max还是凸函数,三分即可. 然而无良出题人既卡精度又卡时间.众 ...
- 2018.10.15 bzoj4570: [Scoi2016]妖怪(凸包)
传送门 不得不说这题有点东西啊. 看到题第一眼二分,用二次函数求范围来进行checkcheckcheck,20分滚粗了233. 于是开始思考正解. 发现可以把每只怪物的二元组属性看成二维坐标. 这时对 ...
- BZOJ 4570: [Scoi2016]妖怪
二次联通门 : BZOJ 4570: [Scoi2016]妖怪 二次联通门 : luogu P3291 [SCOI2016]妖怪 LibreOJ : LibreOJ #2015. 「SCOI2016 ...
- 【bzoj4570 scoi2016】妖怪
题目描述 邱老师是妖怪爱好者,他有n只妖怪,每只妖怪有攻击力atk和防御力dnf两种属性.邱老师立志成为妖怪大师,于是他从真新镇出发,踏上未知的旅途,见识不同的风景. 环境对妖怪的战斗力有很大影响,在 ...
- BZOJ4570:[SCOI2016]妖怪——题解
https://www.lydsy.com/JudgeOnline/problem.php?id=4570 邱老师是妖怪爱好者,他有n只妖怪,每只妖怪有攻击力atk和防御力dnf两种属性.邱老师立志成 ...
- 【BZOJ4570】 [Scoi2016]妖怪
Description 邱老师是妖怪爱好者,他有n只妖怪,每只妖怪有攻击力atk和防御力dnf两种属性.邱老师立志成为妖怪大师,于 是他从真新镇出发,踏上未知的旅途,见识不同的风景.环境对妖怪的战斗力 ...
- P3291 [SCOI2016]妖怪
传送门 我数学的确白学了--这种题目竟然一点思路都没有-- 首先可以把每个妖怪看成二维平面上的一个点,那么每一个环境\((a,b)\)就可以看成一条斜率\(k=-\frac{b}{a}\)的过该点的直 ...
随机推荐
- ORA-01012: not logged on
关于ORA-01012这个错误,惜分飞的博客ORA-01012: not logged on里面已经做了一些介绍,原因就不多说了,看看他的描述说明: 现象说明: 1)终于发现了ORA-01012错误, ...
- Snapshot Instance 操作详解 - 每天5分钟玩转 OpenStack(36)
本节我们通过日志详细讨论 instance 的 snapshot 操作. 有时候操作系统损坏得很严重,通过 Rescue 操作无法修复,那么我们就得考虑通过备份恢复了.当然前提是我们之前对instan ...
- 实战:rsync+inotify实现数据实时同步
Linux 内核从 2.6.13 版本开始提供了 inotify 通知接口,用来监控文件系统的各种变化情况,如文件存取.删除.移动等.利用这一机制,可以非常方便地实现文件异动告警.增量备份,并针对目录 ...
- Redmine 插件安装
将对应的插件都复制进redmine的plugins 安装对应所需要的GEMS bundle install --without development test rmagick 执行插件合并 bund ...
- 将Apache手动安装成Windows的服务
将Apache手动安装成Windows的服务 可以选择在安装Apache时自动将其安装为一个服务.如果选择"for all users",那么Apache将会被安装为服务. 如果选 ...
- 【转帖】分享一个迅为4412开发板OTG烧录批处理文件
平台:iTOP-4412开发板 Bat 功能: 1.可以分条的执行烧录,不需要每次烧录都去复制命令 2.可以批量烧录 开发板系统烧录批处理文件,请将此文件放置在fastboot程序同目录下,下载地址: ...
- UI自动化,你值得拥有
去年春节联欢晚会,为了那张“敬业福”,全家都卯足了劲儿“咻一咻”,连节目都顾不上看了.当时我就想,要是能自动化该多好,不停点击屏幕,屏幕不疼手还疼呢,何况还不好分心,生怕错过了“敬业福”.玩“咻一咻” ...
- Java面向对象 第一章 面向对象开发方法概述
一.软件开发经历的生命周期: ①软件分析 ②软件设计 ③软件编码 ④ 软件测试 ⑤ 软件部署 ⑥软件维护 二.为了提高软件开发效率,降低软件开发成本,一个优良的软件系统应该具备以下特点: ① 可重用性 ...
- Winform(DataGridView)控件及通过此控件中实现增删改查
DataGridView:显示数据表,通过此控件中可以实现连接数据库,实现数据的增删改查 一.后台数据绑定: List<xxx> list = new List<xxx> ...
- JAVA object
1.toString 把对象变成字符串 对于一个引用型变量toString写不写都是一样的 2.equals() 比较两个对象的引用是否一样. 3. public class EqualsTest ...