BZOJ2216 [Poi2011]Lightning Conductor

考虑分j在i前边和后边分别算然后取max

前边后边都是一样的，考虑前边的情况则f[i]=max(1<=j<i)(a[j]+sqrt(i-j))

因为根号x函数差分之后单调递减，所以有严格意义上的决策单调性，即对于i的任意一个最优决策点k，j的任意一个最优决策点l，i<j，都有k<=l

注意是根号函数，double类型的，不是根号上取整函数！要最后再取整，过程中取整没法搞的，因为要是根号上取整的话虽然也是决策单调性，但是如果i有多个决策点都可行的话不一定所有决策点都满足单调性，然而我们无法判断哪个是有单调性的，所以必须用double

然后分治优化就好了

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<map>
#include<set>
#include<bitset>
#include<queue>
#include<stack>
using namespace std;
#define MAXN 500010
#define MAXM 1010
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
#define ll long long
int n;
int a[MAXN];
double f1[MAXN],f2[MAXN];
void dp(int l,int r,int L1,int R1,int L2,int R2){
	if(l>r){
		return ;
	}
	int i;
	int mid=l+r>>1;
	int k1,k2;
	f1[mid]=f2[mid]=-INF;
	for(i=L1;i<=R1&&i<=mid;i++){
		if(a[i]+sqrt(mid-i)>=f1[mid]){
			f1[mid]=a[i]+sqrt(mid-i);
			k1=i;
		}
	}
	for(i=R2;i>=L2&&i>=mid;i--){
		if(a[i]+sqrt(i-mid)>=f2[mid]){
			f2[mid]=a[i]+sqrt(i-mid);
			k2=i;
		}
	}
	dp(l,mid-1,L1,k1,L2,k2);
	dp(mid+1,r,k1,R1,k2,R2);
}
int main(){
	int i;
	scanf("%d",&n);
	for(i=1;i<=n;i++){
		scanf("%d",&a[i]);
	}
	dp(1,n,1,n,1,n);
	for(i=1;i<=n;i++){
		printf("%d\n",int(ceil(max(f1[i],f2[i])))-a[i]);
	}
	return 0;
}

/*
9
7 9 7 10 10 4 9 6 1 

*/