P3291-[SCOI2016]妖怪【凸壳】
正题
题目链接:https://www.luogu.com.cn/problem/P3291
题目大意
给出 \(n\) 个数字对 \((atk,dnf)\),求一个\((a,b)\)。
对于每个数字对可以选择任意一个实数\(k\)让其变为\((atk+k\times a,dnf-k\times a)\),但是操作完之后两个数字都非负。记\(atk/dnf(a,b)\)表示在\((a,b)\)下\(atk/dnf\)的最大值。
然后要求最小化\(max\{atk_i(a,b),dnf_i(a,b)\}\)。
\(1\leq n\leq 10^6,1\leq atk,dnf\leq 10^8\)
解题思路
首先视\((atk,dnf)\)为一个点的话,那么对于任意一个\((a,b)\)答案肯定是在上凸壳上的。
然后考虑实际上我们并不需要用到\((a,b)\)只需考虑\(\frac{b}{a}\)的值,定义\(k=\frac{b}{a}\)
然后就是要求最小化(用\(a_i\)代\(atk_i\),\(d_i\)代\(dnf_i\))
\]
考虑这个点在\(k\)的哪些区间由它取到最大值,对于一个\(j\)需要满足
\]
化一下
\]
然后就是一个二次不等式,并且考虑到\(j\)只需考虑凸壳上\(i\)左右连接的两个点,解出来我们可以得到\(k\)的合法范围。
然后上面那个是一个对钩函数,现在只需在这个范围内求这个对钩函数的最小值就好了。
时间复杂度\(O(n\log n)\)
code
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=1e6+10;
struct node{
double x,y;
}p[N],s[N];
int n,top;double ans;
bool calc(double a,double b,double c,double &l,double &r){
double d=b*b-4.0*a*c;
if(d<0)return 0;d=sqrt(d);
double x0=(-b-d)/(2*a),x1=(-b+d)/(2*a);
if(x0>x1)swap(x0,x1);l=x0;r=x1;
return 1;
}
bool cmp(node x,node y)
{return (x.x==y.x)?(x.y>y.y):(x.x<y.x);}
double solpe(node x,node y)
{return (y.y-x.y)/(y.x-x.x);}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%lf%lf",&p[i].x,&p[i].y);
sort(p+1,p+1+n,cmp);
for(int i=1;i<=n;i++){
while(top>1&&solpe(s[top-1],s[top])<=solpe(s[top-1],p[i]))top--;
s[++top]=p[i];
}
ans=1e18;
for(int i=1;i<=top;i++){
double z=sqrt(s[i].y/s[i].x);
double l=0,r=1e18,L=1,R=1;bool flag=1;
if(i>1)calc(s[i].x-s[i-1].x,s[i].x-s[i-1].x+s[i].y-s[i-1].y,s[i].y-s[i-1].y,L,R);
if(i<top)flag&=calc(s[i].x-s[i+1].x,s[i].x-s[i+1].x+s[i].y-s[i+1].y,s[i].y-s[i+1].y,l,r);
if(!flag)continue;
if(L<l)l=max(R,l);if(R>r)r=min(r,L);
if(l>r||r<=0)continue;z=max(z,l);z=min(z,r);
if(z>L&&z<R){
if(L>=l)ans=min(ans,s[i].x+s[i].y+s[i].x*L+s[i].y/L);
if(R<=r)ans=min(ans,s[i].x+s[i].y+s[i].x*R+s[i].y/R);
}
else ans=min(ans,s[i].x+s[i].y+s[i].x*z+s[i].y/z);
}
printf("%.4lf\n",ans);
return 0;
}
P3291-[SCOI2016]妖怪【凸壳】的更多相关文章
- P3291 [SCOI2016]妖怪
传送门 我数学的确白学了--这种题目竟然一点思路都没有-- 首先可以把每个妖怪看成二维平面上的一个点,那么每一个环境\((a,b)\)就可以看成一条斜率\(k=-\frac{b}{a}\)的过该点的直 ...
- BZOJ 4570: [Scoi2016]妖怪
二次联通门 : BZOJ 4570: [Scoi2016]妖怪 二次联通门 : luogu P3291 [SCOI2016]妖怪 LibreOJ : LibreOJ #2015. 「SCOI2016 ...
- BZOJ4570:[SCOI2016]妖怪——题解
https://www.lydsy.com/JudgeOnline/problem.php?id=4570 邱老师是妖怪爱好者,他有n只妖怪,每只妖怪有攻击力atk和防御力dnf两种属性.邱老师立志成 ...
- 【BZOJ4570】 [Scoi2016]妖怪
Description 邱老师是妖怪爱好者,他有n只妖怪,每只妖怪有攻击力atk和防御力dnf两种属性.邱老师立志成为妖怪大师,于 是他从真新镇出发,踏上未知的旅途,见识不同的风景.环境对妖怪的战斗力 ...
- BZOJ4570: [Scoi2016]妖怪
题目传送门 4570: [Scoi2016]妖怪 Time Limit: 10 Sec Memory Limit: 64 MB Submit: 491 Solved: 125 [Submit][Sta ...
- BZOJ 3672 [Noi2014]购票 (熟练剖分+凸壳维护)
题目链接:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3672 题意:给出一棵有根树(1为根),边有长度.每个点u有三个属性(len[u], ...
- bzoj 3165: [Heoi2013]Segment 动态凸壳
3165: [Heoi2013]Segment Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 202 Solved: 89[Submit][Stat ...
- [CF1137E]Train Car Selection[维护凸壳]
题意 题目链接 分析 首先,如果加到了车头所有之前的车厢都不可能成为答案. 如果加到了车尾,容易发现对于 \(x_2<x_3\) 而言在某个时刻会出现 2 又比 3 优的情况. 具体来讲,如果存 ...
- BZOJ.1007.[HNOI2008]水平可见直线(凸壳 单调栈)
题目链接 可以看出我们是要维护一个下凸壳. 先对斜率从小到大排序.斜率最大.最小的直线是一定会保留的,因为这是凸壳最边上的两段. 维护一个单调栈,栈中为当前可见直线(按照斜率排序). 当加入一条直线l ...
随机推荐
- 【AIOT】智能感知--人
From: https://liudongdong1.github.io/ 1. 人体存在感知 目标:检测环境中的所有人体,标记出每个人体的坐标位置:不限人体数量,适应中低空斜拍.人体轻度遮挡.截断等 ...
- centos7上安装redis6-0-5
下载tar包 wget http://download.redis.io/releases/redis-6.0.5.tar.gz 解压tar包 tar -zxvf redis-6.0.5.tar.gz ...
- zabbix监控rabbitmq队列消费状态
使用rabbitmqctl 管理 mq -n 指定节点 [root@logging-master zabbix]# rabbitmqctl -n rabbit@localhost list_queue ...
- C# 获取计算机的硬件、操作系统信息
C# 获取计算机的硬件.操作系统信息 获得信息的命名空间是:System.Management 创建ManagementObjectSearcher对象获取相关信息 eg: using Syste ...
- C#调用C++ dll中返回值为字符串的函数问题
C#调用C++ dll函数,如果返回值为字符串,我们使用string去接收就会报错,因为C++返回的是char*,是个指针,所以c# 要用 IntPtr 来接收. C++: //预编译的标头 .h e ...
- 大数据Shell编程 之 常用正则表达式
使用man grep查看grep的帮助文档,有如下内容: 可以看出,正则表达式由三类,分别是 基本的正则表达式(Basic Regular Expression 又叫 Basic RegEx 简称 B ...
- MySQL常用权限操作
MySQL常用权限操作 ** ubuntu mysql 8.0.21修改root密码 ** 1.查看默认安装密码: sudo cat /etc/mysql/debian.cnf 2. 登录mysql ...
- mybaits源码分析--日志模块(四)
一.日志模块 首先日志在我们开发过程中占据了一个非常重要的地位,是开发和运维管理之间的桥梁,在Java中的日志框架也非常多,Log4j,Log4j2,Apache Commons Log,java.u ...
- Centos7部署Open-falcon 0.2.0
官方和github上都有教程,但是对于我来说有的部署内容较为陌生,有点错误官方也未在教程中说明,故在此记录方便以后快速部署,本文部署的时间是2018/10/10. 虽然open-falcon是采用了前 ...
- 一种封装Retrofit的方法,可以自动解析Gson,回避Method return type must not include a type variable or wildcard: retrofit2.Call<T>的问题
封装目的:屏蔽底层实现,提供统一接口,并支持Gson自动转化 最初封装: //请求方法 interface RequestListener { interface PostListener { @PO ...