BZOJ 3203: [Sdoi2013]保护出题人

题意：

考虑暴力，ans[i]=max((sum[i]-sum[j-1])/(x[i]+i*d-j*d))(其中sum记录僵尸血量的前缀和）。接下来有一个巧妙的想法，我们观察式子后发现可以抽象成平面点对的斜率，对于一个i,可以表示成（sum[i-1],i*d),因为取最大值，所以就像做DP的斜率优化一样维护单调的凸壳，然后用单调队列进行删除（发现taii和i的交点小于tail和tail-1的交点，则tail无用，弹出），之后我们要在凸壳上找最大斜率，根据sum[i-1],i*d的递增性质（不会证= =），可以用三分解决。

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
using namespace std;
const int maxn=100000+10;
long long x[maxn],a[maxn],sum[maxn],d;
int st[maxn],top,n;
double ans;
double calc(int i,int j)
{
  double x=(double)d*i-(double)d*j;
  double y=sum[i-1]-sum[j-1];
  return y/x;
}
double work(int i,int j)
{
  double res=(double)(sum[i]-sum[j-1])/(double)(x[i]+d*i-d*j);
  return res;
}
int main()
{
  //freopen("3203.in","r",stdin);
  //freopen("3203.out","w",stdout);
  scanf("%d%lld",&n,&d);
  for(int i=1;i<=n;i++)
  {
    scanf("%lld%lld",&a[i],&x[i]);
    sum[i]=sum[i-1]+a[i];
  }
  for(int i=1;i<=n;i++)
  {
    while(top>1&&calc(i,st[top])<calc(st[top],st[top-1])) top--;
    st[++top]=i;
    int l=1,r=top;
    while(l!=r)
    {
      int t=(r-l+1)/3;  
      int x1=l-1+t;  
      int x2=l-1+2*t;  
      if(x1==x2) break;
      if(work(i,st[x1])<work(i,st[x2])) l=x1+1;
      else r=x2-1;
    }
    ans+=max(work(i,st[l]),work(i,st[r]));
  }
  printf("%.0lf\n",ans);
  return 0;
}