考虑分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
*/