正题

题目链接: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\))

\[a_i+b_i+a_ik+b_i\frac{1}{k}
\]

考虑这个点在\(k\)的哪些区间由它取到最大值,对于一个\(j\)需要满足

\[a_i+b_i+a_ik+b_i\frac{1}{k}>a_j+b_j+a_jk+b_j\frac{1}{k}
\]

化一下

\[(a_i-a_j)k^2+(a_i-a_j+b_i-b_j)k+(b_i-b_j)>0
\]

然后就是一个二次不等式,并且考虑到\(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]妖怪【凸壳】的更多相关文章

  1. P3291 [SCOI2016]妖怪

    传送门 我数学的确白学了--这种题目竟然一点思路都没有-- 首先可以把每个妖怪看成二维平面上的一个点,那么每一个环境\((a,b)\)就可以看成一条斜率\(k=-\frac{b}{a}\)的过该点的直 ...

  2. BZOJ 4570: [Scoi2016]妖怪

    二次联通门 : BZOJ 4570: [Scoi2016]妖怪 二次联通门 : luogu P3291 [SCOI2016]妖怪 LibreOJ : LibreOJ  #2015. 「SCOI2016 ...

  3. BZOJ4570:[SCOI2016]妖怪——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=4570 邱老师是妖怪爱好者,他有n只妖怪,每只妖怪有攻击力atk和防御力dnf两种属性.邱老师立志成 ...

  4. 【BZOJ4570】 [Scoi2016]妖怪

    Description 邱老师是妖怪爱好者,他有n只妖怪,每只妖怪有攻击力atk和防御力dnf两种属性.邱老师立志成为妖怪大师,于 是他从真新镇出发,踏上未知的旅途,见识不同的风景.环境对妖怪的战斗力 ...

  5. BZOJ4570: [Scoi2016]妖怪

    题目传送门 4570: [Scoi2016]妖怪 Time Limit: 10 Sec Memory Limit: 64 MB Submit: 491 Solved: 125 [Submit][Sta ...

  6. BZOJ 3672 [Noi2014]购票 (熟练剖分+凸壳维护)

    题目链接:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3672 题意:给出一棵有根树(1为根),边有长度.每个点u有三个属性(len[u], ...

  7. bzoj 3165: [Heoi2013]Segment 动态凸壳

    3165: [Heoi2013]Segment Time Limit: 40 Sec  Memory Limit: 256 MBSubmit: 202  Solved: 89[Submit][Stat ...

  8. [CF1137E]Train Car Selection[维护凸壳]

    题意 题目链接 分析 首先,如果加到了车头所有之前的车厢都不可能成为答案. 如果加到了车尾,容易发现对于 \(x_2<x_3\) 而言在某个时刻会出现 2 又比 3 优的情况. 具体来讲,如果存 ...

  9. BZOJ.1007.[HNOI2008]水平可见直线(凸壳 单调栈)

    题目链接 可以看出我们是要维护一个下凸壳. 先对斜率从小到大排序.斜率最大.最小的直线是一定会保留的,因为这是凸壳最边上的两段. 维护一个单调栈,栈中为当前可见直线(按照斜率排序). 当加入一条直线l ...

随机推荐

  1. 【AIOT】智能感知--人

    From: https://liudongdong1.github.io/ 1. 人体存在感知 目标:检测环境中的所有人体,标记出每个人体的坐标位置:不限人体数量,适应中低空斜拍.人体轻度遮挡.截断等 ...

  2. 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 ...

  3. zabbix监控rabbitmq队列消费状态

    使用rabbitmqctl 管理 mq -n 指定节点 [root@logging-master zabbix]# rabbitmqctl -n rabbit@localhost list_queue ...

  4. C# 获取计算机的硬件、操作系统信息

    C# 获取计算机的硬件.操作系统信息   获得信息的命名空间是:System.Management 创建ManagementObjectSearcher对象获取相关信息 eg: using Syste ...

  5. C#调用C++ dll中返回值为字符串的函数问题

    C#调用C++ dll函数,如果返回值为字符串,我们使用string去接收就会报错,因为C++返回的是char*,是个指针,所以c# 要用 IntPtr 来接收. C++: //预编译的标头 .h e ...

  6. 大数据Shell编程 之 常用正则表达式

    使用man grep查看grep的帮助文档,有如下内容: 可以看出,正则表达式由三类,分别是 基本的正则表达式(Basic Regular Expression 又叫 Basic RegEx 简称 B ...

  7. MySQL常用权限操作

    MySQL常用权限操作 ** ubuntu mysql 8.0.21修改root密码 ** 1.查看默认安装密码: sudo cat /etc/mysql/debian.cnf 2. 登录mysql ...

  8. mybaits源码分析--日志模块(四)

    一.日志模块 首先日志在我们开发过程中占据了一个非常重要的地位,是开发和运维管理之间的桥梁,在Java中的日志框架也非常多,Log4j,Log4j2,Apache Commons Log,java.u ...

  9. Centos7部署Open-falcon 0.2.0

    官方和github上都有教程,但是对于我来说有的部署内容较为陌生,有点错误官方也未在教程中说明,故在此记录方便以后快速部署,本文部署的时间是2018/10/10. 虽然open-falcon是采用了前 ...

  10. 一种封装Retrofit的方法,可以自动解析Gson,回避Method return type must not include a type variable or wildcard: retrofit2.Call<T>的问题

    封装目的:屏蔽底层实现,提供统一接口,并支持Gson自动转化 最初封装: //请求方法 interface RequestListener { interface PostListener { @PO ...