BZOJ3203 SDOI2013保护出题人（三分）

　　给a做一个前缀和，那么现在每次所查询的就是(s_n-s_k)/(b_n+nd-(k+1)d)的最大值。这个式子可以看成是(b_n+nd,s_n)和((k+1)d,s_k)所成直线的斜率。

　　脑补一条直线不断减小斜率，容易发现可能成为最大值的点一定在下凸壳上。并且凸壳上的点和该点间的斜率变化情况是一个凸函数。于是维护出凸壳，在凸壳上三分即可。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 100010
#define ll long long
ll read()
{
    ll x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
int n,q[N],tail;
double ans=;
ll d,a[N],b[N];
double calc(int x,int y)
{
    return 1.0*(a[y]-a[x])/((y-x)*d);
}
double getans(int n,int k)
{
    return 1.0*(a[n]-a[k])/(b[n]+(n-k-)*d);
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj3203.in","r",stdin);
    freopen("bzoj3203.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),d=read();
    for (int i=;i<=n;i++) a[i]=read(),b[i]=read();
    q[tail=]=;
    for (int i=;i<=n;i++)
    {
        a[i]+=a[i-];
        int l=,r=tail;
        while (l+<r)
        {
            int mid=l+r>>;
            double x=getans(i,q[mid]),y=getans(i,q[mid+]);
            if (x>y) r=mid+;else l=mid;
        }
        int mx=l;
        if (l+<=r&&getans(i,q[l+])>getans(i,q[mx])) mx=l+;
        if (l+<=r&&getans(i,q[l+])>getans(i,q[mx])) mx=l+;
        ans+=getans(i,q[mx]);
        while (tail>=&&calc(q[tail-],q[tail])>calc(q[tail],i)) tail--;
        q[++tail]=i;
    }
    printf("%0.lf",ans);
    return ;
}