【BZOJ2216】Lightning Conductor（动态规划）
题面
题解
(sqrt {|i-j|})似乎没什么意义，只需要从前往后做一次再从后往前做一次就好了。
只考虑从前往后，把给定的式子移项，可以得到
(pge a[j]-a[i]+sqrt{i-j})
而(a[i])是当前的枚举的位置(i)的值，这个是不会变化的。
所以要求的就是(max(a[j]-sqrt{i-j}))
画出(sqrt x)的函数图像，是一个增长率越来越慢的函数。
那么，如果当前转移(i)的时候，(j)优于(k)((klt j))，
那么接下来转移(i+1)的时候同样(j)更优。
既然具有了决策单调性，直接二分就好了。
维护一个单调队列，存下当前位置(j)转移到哪些区间是更优的，
不难发现这个区间一定是([l,n])，当然也可能是个空区间。
插入的时候也二分修改一下就好了。有点类似于诗人小G那题。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define ll long long
#define RG register
#define MAX 500500
inline int read()
{
RG int x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
struct Node{int i,l,r;}Q[MAX];
int h,t;
int n,a[MAX];
double f[MAX];
double Trans(int j,int i){return a[j]+sqrt(i-j);}
void Work()
{
Q[h=t=1]=(Node){1,2,n};
for(int i=2;i<=n;++i)
{
while(h<=t&&Q[h].r<i)++h;Q[h].l=i;
f[i]=max(f[i],Trans(Q[h].i,i)-a[i]);
while(h<=t&&Trans(Q[t].i,Q[t].l)<Trans(i,Q[t].l))--t;
if(h>t){Q[++t]=(Node){i,i,n};continue;}
int l=Q[t].l,r=Q[t].r,ret=Q[t].r+1;
while(l<=r)
{
int mid=(l+r)>>1;
if(Trans(Q[t].i,mid)<Trans(i,mid))ret=mid,r=mid-1;
else l=mid+1;
}
if(ret!=Q[t].l)Q[t].r=ret-1;else --t;
if(ret<=n)Q[++t]=(Node){i,ret,n};
}
}
int main()
{
n=read();
for(int i=1;i<=n;++i)a[i]=read();
Work();reverse(&a[1],&a[n+1]);reverse(&f[1],&f[n+1]);
Work();reverse(&f[1],&f[n+1]);
for(int i=1;i<=n;++i)printf("%d
",(int)(ceil(f[i])));
return 0;
}