bzoj4570 [Scoi2016]妖怪（凸包+对勾函数最值）

bzoj4570 [Scoi2016]妖怪

原题地址：http://www.lydsy.com/JudgeOnline/problem.php?id=4570

题意：
邱老师有n只妖怪，每只妖怪有攻击力atk和防御力dnf两种属性。
环境由a，b两个参数定义，在某种环境中，妖怪可以降低自己k×a点攻击力，提升k×b点防御力，或者，提升自己k×a点攻击力，降低k×b点防御力，a，b属于正实数，k为任意实数，但是atk和dnf必须始终非负。
妖怪在环境(a,b)中的战斗力为妖怪在该种环境中能达到的最大攻击力和最大防御力之和。
邱老师想知道在最为不利的情况下，他的n只妖怪能够达到的最强战斗力值，即存在一组正实数(a,b)使得n只妖怪在该环境下最强战斗力最低。

数据范围
1≤n≤10^6, 0＜atk,dnf≤10^8

题解：
好题。
首先，一个非常重要的转化，需要把原题抽象成：
假如我们认为atk为x，dnf为y，那么每个妖怪都是一个坐标为（x，y）的点。
给出一个环境变量（a,b），发现这个提升自己k×a点攻击力，降低k×b点防御力的过程，就是一条过（x，y），斜率$k=-\frac{b}{a}$的直线，
最大战斗力就是横纵截距之和。

问题转化为：
给出平面上n个点，求存在一个斜率k（k<0）使经过这些点的斜率为k的直线中的最大横纵截距和最小时，最大的横纵斜率和。

发现在斜率改变时，有些点永远不可能成为横纵截距和最大的点，能够成为答案的只有上凸壳上的点。
于是我们搞出这个上凸壳。
依次枚举上凸壳上的点，更新答案。
对于每个点（x，y），给出斜率k，横纵截距之和最小为：
$x+y-k*x-\frac{y}k$
在$k=-\sqrt{\frac{y}x}$时取到最值。
假如说对于枚举到的这个点（x，y），它的最小值斜率为k，它与它前后的两个点构成的斜率为k1，k2。
那么当k的直线不指向凸包内，就用k更新答案，（即k1<=k<=k2）
否则用k1，k2更新答案（因为函数的最小值取不到了，往两边走只会使答案变得更大，因此刚好卡着凸包是移动最小的）

（用k1，k2更新答案是必要的，因为虽然凸包上的点表示两个不同对勾函数的交点，但并不是交点处的最大值就是最小的）

代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
const int N=1000006;
const double EPS=1e-8;
struct Vec
{
    double x,y;
    Vec(){}
    Vec(double x,double y):x(x),y(y){}
}P[N],cvx[N];
Vec operator+(const Vec &A,const Vec &B) {return Vec(A.x+B.x,A.y+B.y);}
Vec operator-(const Vec &A,const Vec &B) {return Vec(A.x-B.x,A.y-B.y);}
Vec operator*(const Vec &A,double s) {return Vec(A.x*s,A.y*s);}
Vec operator/(const Vec &A,double s) {return Vec(A.x/s,A.y/s);}
double dot(const Vec &A,const Vec &B) {return A.x*B.x+A.y*B.y;}
double cross(const Vec &A,const Vec &B) {return A.x*B.y-A.y*B.x;}
int n;
int sign(double x)
{
    if(fabs(x)<=EPS) return 0;
    return x<0?-1:1;
}
bool cmp(const Vec &A,const Vec &B) {return (A.x!=B.x)?A.x<B.x:A.y<B.y;}
bool onleft(const Vec &A,const Vec &B,const Vec &C) {return sign(cross(B-A,C-A))>=0;}
int getcvx()
{
    int t=0; sort(P+1,P+n+1,cmp);
    for(int i=1;i<=n;i++)
    {
        while(t>=2&&onleft(cvx[t-1],cvx[t],P[i])) t--;
        cvx[++t]=P[i];
    }
    return t;
}
double min_slope(Vec A) {return -sqrt(A.x/A.y);}
double slope(Vec &A,Vec &B) {return (A.x-B.x==0)?1e60:(A.y-B.y)/(A.x-B.x);}
double cal(Vec &A,double k) 
{
    if(k>=0) return 1e20;
    return A.x+A.y-A.x*k-A.y/k;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) {double x,y; scanf("%lf%lf",&x,&y);P[i]=Vec(x,y);}
    n=getcvx(); double ans=1e20;
    if(n==1) printf("%0.4lf\n",cal(cvx[1],min_slope(cvx[1])));
    else
    {
        for(int i=2;i<=n;i++) 
        {
            double k1=slope(cvx[i],cvx[i-1]);
            double k2=(i==n)?-1e60:slope(cvx[i+1],cvx[i]);
            ans=min(ans,cal(cvx[i],k1));
            double k=min_slope(cvx[i]);
            if(k1>=k&&k>=k2) ans=min(ans,cal(cvx[i],k));
        }
        double k=min_slope(cvx[1]); 
        double k2=slope(cvx[1],cvx[2]);
        if(k>=k2) ans=min(ans,cal(cvx[1],k));
        printf("%0.4lf\n",ans);
    }
    return 0;
}