【BZOJ2216】Lightning Conductor（动态规划）

最新推荐文章于 2025-08-09 11:46:48 发布

weixin_30767921

最新推荐文章于 2025-08-09 11:46:48 发布

阅读量56

点赞数

CC 4.0 BY-SA版权

文章标签：数据结构与算法

原文链接：http://www.cnblogs.com/cjyyb/p/9308498.html

本文介绍了解决BZOJ2216问题的一种动态规划方法，通过二分搜索和单调队列实现高效求解。文章详细解释了算法原理及其实现过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

【BZOJ2216】Lightning Conductor（动态规划）

题面

BZOJ，然而是权限题
 洛谷

题解

\(\sqrt {|i-j|}\)似乎没什么意义，只需要从前往后做一次再从后往前做一次就好了。
只考虑从前往后，把给定的式子移项，可以得到
\(p\ge a[j]-a[i]+\sqrt{i-j}\)
而\(a[i]\)是当前的枚举的位置\(i\)的值，这个是不会变化的。
所以要求的就是\(max(a[j]-\sqrt{i-j})\)
画出\(\sqrt x\)的函数图像，是一个增长率越来越慢的函数。
那么，如果当前转移\(i\)的时候，\(j\)优于\(k\)(\(k\lt 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\n",(int)(ceil(f[i])));
    return 0;
}

转载于:https://www.cnblogs.com/cjyyb/p/9308498.html