bzoj 3203: [Sdoi2013]保护出题人

题意

出题人在玩植物大战僵尸，卡关了。。
让你帮助他(雾)

题解

看了半天
才知道这是个超级豌豆射手
他的每一发子弹，如果打死了一个僵尸，剩下多余伤害可以给后面的僵尸承担。。
我就说怎么怪怪的

可以得到公式，第i轮的答案是

$$max((sum[i]-sum[j-1])/(x[i]+(i-j)*d)$$
这个形式十分像一个斜率的形式
可以把他看做斜率
又因为，sum的单调性
所以决策点肯定是在一个下凸的凸包上。。
维护一个凸包，然后在凸包上面三分就可以了
显然满足三分性
CODE:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long double LB;
typedef long long LL;
const LL N=1e5+5;
const LB eps=1e-15;
LL n,d;
LB a[N],x[N];
LB sum[N];
LL sta[N],top=0;
LB xielv (LL x,LL y)
{
    LB x1=d*x,y1=sum[x-1];
    LB x2=d*y,y2=sum[y-1];
    return (y1-y2)/(x1-x2);
}
LB calc (LL i,LL j){return (sum[i]-sum[j-1])/(x[i]+d*i-d*j);}
int main()
{
    scanf("%lld%lld",&n,&d);
    for (LL u=1;u<=n;u++)
    {
        scanf("%Lf%Lf",&a[u],&x[u]);
        sum[u]=sum[u-1]+a[u];
    }
    LB ans=0;
    for (LL u=1;u<=n;u++)
    {
        while (top>1&&xielv(u,sta[top])-xielv(sta[top],sta[top-1])<-eps) top--;
        sta[++top]=u;
        LL l=1,r=top;
        while (l+1<r)
        {
            LL mid=(r-l+1)/3;
            LL ll=l+mid-1,rr=ll+mid;
            LB lans=calc(u,sta[ll]);
            LB rans=calc(u,sta[rr]);
            if (lans-rans<-eps) l=ll+1;
            else r=rr-1;
        }
        ans=ans+max(calc(u,sta[l]),calc(u,sta[r]));
    }
    printf("%.0Lf\n",ans);
    return 0;
}