bzoj2216: [Poi2011]Lightning Conductor 决策单调性

这是一个非常经典的模型了，1D1D的动态规划
这类问题最基础的是斜率优化。难一点的就很麻烦了……什么用平衡树维护凸包之类的。不过我完全不会
这题是我见到的第一道非斜率优化的1D1D了……

p=F[i]=max{a[j]+sqrt(abs(i-j))}-a[i]

注意到函数f(x)=sqrt(x)-sqrt(x-1)单调递减，所以若k<j<i且对于i而言j比k优，则k没用了
然后可以yy出很多东西
比如对于i，最优解是i的必然是一段区间。
所以我们维护一个栈，栈里每一个元素都代表一段区间，表示这段区间的最优解都来自这个元素p。所以一个元素是个三元组(p,l,r)
要保证对于i最优决策是队首区间
然后插入i，如果他能改变队列的情况，那么i对应的区间必然是[l,n]
所以步骤如下：
    如果对于n，当前的i比队尾更优，就不断的将队尾区间弹出栈。将区间(p,l,r)弹出来的条件是对于l,i比p更优。
    现在我们的队尾(p,l,r)满足，对于l，p比i优；对于r，i比p优。我们可以通过二分查找，找到一个临界点q,然后将队尾变成(p,l,q-1)并将元素(i,q+1,n)插入队尾
    这只是大概思路，实现过程中要注意队列为空

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n;
int a[500005];
double f[500005],g[500005];
struct data
{
    int l,r,id;
}q[500005];
double cal(int j,int i)
{
    return a[j]+sqrt(abs(i-j))-a[i];
}
int find(data t,int x)
{
    int l=t.l,r=t.r;
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(cal(t.id,mid)>cal(x,mid))
            l=mid+1;
        else r=mid-1;
    }
    return l;
}
void dp(double *F)
{
    int head=1,tail=0;
    for(int i=1;i<=n;i++)
    {
        q[head].l++;
        if(head<=tail&&q[head].r<q[head].l)head++;
        if(head>tail||cal(i,n)>cal(q[tail].id,n))
        {
            while(head<=tail&&cal(i,q[tail].l)>cal(q[tail].id,q[tail].l))tail--;
            if(head>tail) q[++tail]=(data){i,n,i};
            else
            {
                int t=find(q[tail],i);
                q[tail].r=t-1;
                q[++tail]=(data){t,n,i};
            }
        }
        F[i]=cal(q[head].id,i);
    }
}
int main()
{
     n=read();
    for(int i=1;i<=n;i++)a[i]=read();
    dp(f);
    for(int i=1;i<=n/2;i++)swap(a[i],a[n-i+1]);
    dp(g);
    for(int i=1;i<=n;i++)
        printf("%d\n",max(0,(int)ceil(max(f[i],g[n-i+1]))));
    return 0;
}