[bzoj4850][Jsoi2016]灯塔

来自FallDream的博客，未经允许，请勿转载，谢谢。

---

JSOI的国境线上有N一座连续的山峰,其中第ii座的高度是hi??.为了简单起见,我们认为这N座山峰排成了连续一条

直线.如果在第ii座山峰上建立一座高度为p(p≥0)的灯塔,JYY发现,这座灯塔能够照亮第jj座山峰,当且仅当满足如

下不等式:hj≤hi+p-（是减号)sqrt(|i-j|)JSOI国王希望对于每一座山峰,JYY都能提供建造一座能够照亮全部其他山峰的灯

塔所需要的最小高度.你能帮助JYY么?

1< N ≤ 10^5

0 < hi ≤ 10^9

 

话说这个出题人非常的不走心 同bzoj2216 5年前的题 样例都不改就扔了一个假的公式出来。

暴力比较显然，可以预处理rmq，然后根号枚举。

考虑两个点j，k，那么显然一个点分别从两个点转移的情况都是一段区间，可以二分求出最小的满足从j转移比从k转移优的i。

然后开一个单调队列，维护队列中两两相邻元素，前一个比后一个优的最小的i 单调递增。

复杂度nlogn

#include<iostream>
#include<cstdio>
#include<cmath>
#define getchar() (*S++)
#define MN 500000
#define INF 2000000000
char B[<<],*S=B;
using namespace std;
inline int read()
{
    int x =  , f = ; char ch = getchar();
    while(ch < '' || ch > ''){ if(ch == '-') f = -;  ch = getchar();}
    while(ch >= '' && ch <= ''){x = x *  + ch - '';ch = getchar();}
    return x * f;
}

int n,a[MN+],top,tail,q[MN+];
double F[MN+],G[MN+],sq[MN+];
inline int My_abs(int x){return x<?-x:x;}
double Get(int x,int y){return a[x]+sq[My_abs(y-x)];}

int Calc(int x,int y)
{
    int l=y,r=n,mid,ans=INF;
    while(l<=r)
    {
        mid=l+r>>;
        if(Get(y,mid)>=Get(x,mid)) ans=mid,r=mid-;
        else l=mid+;
    }
    return ans;
}

void Solve(double*f)
{
    top=;tail=;
    for(register int i=;i<=n;++i)
    {
        if(top<tail||a[i]>a[q[top]])
        {
            while(top>tail&&Calc(q[top],i)<=Calc(q[top-],q[top])) --top;
            q[++top]=i;
        }
        while(top>tail&&Calc(q[tail],q[tail+])<=i) ++tail;
        f[i]=Get(q[tail],i)-a[i];
    }
}

int main()
{
    fread(B,,<<,stdin);
    n=read();
    for(int i=;i<=n;++i) sq[i]=sqrt(i);
    for(int i=;i<=n;++i) a[i]=read();
    Solve(F);
    for(int i=;i<=n>>;++i) swap(a[i],a[n+-i]);
    Solve(G);
    for(int i=;i<=n;++i) printf("%d\n",max(,(int)ceil(max(F[i],G[n+-i]))));
    return ;
}