[bzoj3203][SDOI2013]保护出题人

3203: [Sdoi2013]保护出题人

Time Limit: 3 Sec Memory Limit: 128 MB
Submit: 389 Solved: 218
[Submit][Status][Discuss]
Description

Input

第一行两个空格隔开的正整数n和d，分别表示关数和相邻僵尸间的距离。接下来n行每行两个空格隔开的正整数，第i + 1行为Ai和 Xi，分别表示相比上一关在僵尸队列排头增加血量为Ai 点的僵尸，排头僵尸从距离房子Xi米处开始接近。

Output

一个数，n关植物攻击力的最小总和 ，保留到整数。

Sample Input

5 2

3 3

1 1

10 8

4 8

2 3

Sample Output

7
HINT

第一关：距离房子3米处有一只血量3点的僵尸，植物最小攻击力为1.00000； 第二关：距离房子1米处有一只血量1点的僵尸、3米处有血量3点的僵尸，植物最小攻击力为1.33333； 第三关：距离房子8米处有一只血量10点的僵尸、10米处有血量1点的僵尸、12米处有血量3点的僵尸，植物最小攻击力为1.25000； 第四关：距离房子8米处有一只血量4点的僵尸、10米处有血量10点的僵尸、12米处有血量1点的僵尸、14米处有血量3点的僵尸，植物最小攻击力为1.40000； 第五关：距离房子3米处有一只血量2点的僵尸、5米处有血量4点的僵尸、7米处有 血量10点的僵尸、9米处有血量1点的僵尸、11米处有血量3点的僵尸，植物最小攻击力 为2.28571。 植物攻击力的最小总和为7.26905。

对于100%的数据， 1≤n≤10^5，1≤d≤10^12,1≤x≤ 10^12,1≤a≤10^12

感觉这道题还是比较好的。
首先对于第i关来说，对应的yi应该怎样求呢？？
题目中说一个僵尸死了之后，植物会接着攻击后一个僵尸而且僵尸的移动速度是一样的，也就是只有打死了之前的所有僵尸，才能打到最后一个僵尸。这样我们就可以把之前所有僵尸的血全部放到最后一个身上，再除以最后一个到终点的距离，就是速度。
可是这样不能保证前面的僵尸在被打死之前就吃到了植物。所以我们需要将i之前的全部都处理一遍。
yi=max{(sum[i]-sum[j-1])/(x[i]+(i-j)*d)} 1<=j<=i
这样我们处理出了每一关的y。
这样每次做都是n^2的。我们观察一下上面的式子，发现这是个斜率的形式，因为k=(y1-y2)/(x1-x2)。
所以我们设x1=x[i]+i * d,x2=j * d,y1=sum[i],y2=sum[j-1]。
这样就变成了一个球斜率的最大值。
我们维护一个下凸包，每次新加入一个点的时候，因为上面的那个函数是个单峰的，所以在凸包上面三分就行了。

#include
#include
#include
using namespace std;
#define eps 1e-6
#define LL long long
#define mid1 l+(r-l)/3
#define mid2 l+(r-l)/3*2
const int N=100010;
int n,siz;
double ans,d,sum[N],xi[N];
struct Point{double x,y;}q[N];
Point operator - (Point x,Point y){return (Point){x.x-y.x,x.y-y.y};}
inline double Cross(Point x,Point y){return x.x*y.y-x.y*y.x;}
inline double calc(Point x){
    int l=1,r=siz,i;
    double maxn=0,k1,k2;
    while(l+2x.y-q[mid1].y)/(x.x-q[mid1].x);
        k2=(x.y-q[mid2].y)/(x.x-q[mid2].x); 
        maxn=max(maxn,max(k1,k2));
        if(k1else r=mid2;
    }
    for(i=l;i<=r;++i) maxn=max(maxn,(x.y-q[i].y)/(x.x-q[i].x));
    return maxn;
}
int main(){
    int i,j;
    scanf("%d%lfd",&n,&d);
    for(i=1;i<=n;++i){
        scanf("%lf%lf",&sum[i],&xi[i]);
        sum[i]+=sum[i-1];
    }
    ans=sum[1]/xi[1];q[siz=1]=(Point){d,sum[0]};
    for(i=2;i<=n;++i){
        Point x=(Point){i*d,sum[i-1]};
        while(siz>=2&&Cross(q[siz]-q[siz-1],x-q[siz-1])<=0) --siz;
        q[++siz]=x;
        x=(Point){xi[i]+i*d,sum[i]};
        ans+=calc(x);
    }
    printf("%.0f\n",ans);
}