[POI2011]Lightning Conductor

题面在这里

description

已知一个长度为\(n\)的序列\(a_1,a_2,...,a_n\)。

对于每个\(1\le i\le n\),找到最小的非负整数\(p\),

满足对于任意的\(1\le j\le n\),\(a_j\le a_i+p-\sqrt{|i-j|}\)

data range

\[n\le 5\times 10^5,a_i\le 10^9\]

solution

绝对值怎么办?

我们先从左到右\(DP\ j< i\)的部分(此时有\(|i-j|=i-j\)),

再右到左\(DP\ j> i\)的部分(此时有\(|i-j|=j-i\));

那么我们现在考虑这个式子:

\[f[i]=max_{j=1}^{i-1}(a[j]-a[i]+\sqrt{i-j})\]

对于每一个相同的\(i\),

记\(w[j]=\sqrt{i-j}\),那么我们很容易发现这个函数是凸函数

对于两个决策\(j<k\),当\((j-k)\)一定时,我们发现\(w[j]-w[k]\)随着\(i\)的递增而递减

而对于这两个决策,\(j\)比\(k\)优当且仅当对于\(j<k<i\le n\),

\[a[j]+\sqrt{i-j}\ge a[k]+\sqrt{i-k}\]

\[g[j,k]=a[j]-a[k]+w[j]-w[k]\ge 0\]

当\(i\)递增时,\(g[j,k]\)这个函数单调递减,那么存在一个分界点\(i‘\),

使得当\(i\le i‘\)时,\(g[j,k]\ge 0\)(\(j\)更优),当\(i\ge i‘\)时,\(g[j,k]\le 0\)(\(k\)更优);

那么我们知道这个\(DP\)方程具有决策单调性

具体如何实现?

考虑相邻的两个决策\(j,k\),

由于上面的性质,我们可以二分出这两个决策的分界点\(t[j,k]\);

对于相邻的三个决策\(x,y,z\),如果\(t[x,y]\ge t[y,z]\)那么决策\(y\)是没有用的可以扔掉

考虑维护一个单调队列,存放所有合法的决策;

由于我们在维护的同时\(i\)是变化的,

因此我们需要求出每个决策此时转移最优的区间

即单调队列存放的元素为\((pos,l_{pos},r_{pos})\)(注意\(pos < l_{pos}\le r_{pos}\))

当加入新决策\(j\)时,考虑队尾的决策\((k,l_k,r_k)\);

如果\(j\)总是比\(k\)要劣(在这里指对于\(i=n\)仍有\(a[j]+\sqrt{i-j}\le a[k]+\sqrt{i-k}\)),

则不考虑决策\(j\);

否则(此时\(r_j=n\)),如果对于\(i=l_k\)有\(a[j]+\sqrt{i-j}\ge a[k]+\sqrt{i-k}\),

一直弹掉队尾决策并考虑队尾的前一个决策,如果队列为空则压入\((j,j,n)\);

否则(此时存在一个决策\(k\)无法被弹出),考虑二分决策\(k\)和\(j\)的分界点\(mid\)

(本人的定义是,对于\(i\le mid\),有\(a[k]+\sqrt{i-k}> a[j]+\sqrt{i-j}\)),

之后将\(r_k\)修改为\(mid\),压入决策\((j,mid,n)\)。

il dd calc(int j,int i){return a[j]-a[i]+sqrt(i-j);}
//计算j转移至i时对应的DP值
il int divide(int k,int j,int l,int r){
    RG int mid,ans=l;
    while(l<=r){
        mid=(l+r)>>1;
        if(calc(k,mid)>calc(j,mid))
            ans=mid,l=mid+1;
        else r=mid-1;
    }
    return ans;
}
//二分分界点
il void push(int i){
    if(calc(Qx[R],n)>calc(i,n))return;
    while(calc(Qx[R],Ql[R])<calc(i,Ql[R])&&L<=R)R--;
    if(L>R){Ql[++R]=i;Qr[R]=n;Qx[R]=i;return;}
    Qr[R]=divide(Qx[R],i,Ql[R],Qr[R]);
    R++;Ql[R]=Qr[R-1]+1;Qr[R]=n;Qx[R]=i;
}
//插入新决策

当要查询位置\(pos\)时,直接从队头开始查找,

如果队首决策\((k,l_k,r_k)\)中\(r_k<pos\)或队首决策劣于队首后决策则将其弹出;

否则,将其作为此次查询的答案

il int pop(int i){
    while(Qr[L]<i||(L<R&&calc(Qx[L],i)<calc(Qx[L+1],i)))L++;
    return Qx[L];
}
//查询答案

这里是总的\(DP\)模板

il void DP(int *a,int *p,int n){
    L=0;R=-1;
    for(RG int i=1;i<=n;i++){
        if(i!=1)p[i]=ceil(calc(pop(i),i));Ql[L]=i;push(i);
        //答案为非负整数,这里做了上取整
    }
}
//DP序列

code

#include<bits/stdc++.h>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<iomanip>
#include<cstring>
#include<complex>
#include<vector>
#include<cstdio>
#include<string>
#include<bitset>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define mp make_pair
#define pub push_back
#define puf push_front
#define pob pop_back
#define pof pop_front
#define RG register
#define il inline
using namespace std;
typedef unsigned long long ull;
typedef vector<int>VI;
typedef long long ll;
typedef double dd;
const dd eps=1e-10;
const int N=1000000;
const int inf=1e9+7;
il ll read(){
    RG ll data=0,w=1;RG char ch=getchar();
    while(ch!=‘-‘&&(ch<‘0‘||ch>‘9‘))ch=getchar();
    if(ch==‘-‘)w=-1,ch=getchar();
    while(ch<=‘9‘&&ch>=‘0‘)data=data*10+ch-48,ch=getchar();
    return data*w;
}

il void file(){
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
}

int n,a[N],p[N],s[N];
int Qx[N],Ql[N],Qr[N],L,R;

il dd calc(int j,int i){return a[j]-a[i]+sqrt(i-j);}
//计算j转移至i时对应的DP值
il int divide(int k,int j,int l,int r){
    RG int mid,ans=l;
    while(l<=r){
        mid=(l+r)>>1;
        if(calc(k,mid)>calc(j,mid))
            ans=mid,l=mid+1;
        else r=mid-1;
    }
    return ans;
}
//二分分界点
il void push(int i){
    if(calc(Qx[R],n)>calc(i,n))return;
    while(calc(Qx[R],Ql[R])<calc(i,Ql[R])&&L<=R)R--;
    if(L>R){Ql[++R]=i;Qr[R]=n;Qx[R]=i;return;}
    Qr[R]=divide(Qx[R],i,Ql[R],Qr[R]);
    R++;Ql[R]=Qr[R-1]+1;Qr[R]=n;Qx[R]=i;
}
//插入新决策
il int pop(int i){
    while(Qr[L]<i||(L<R&&calc(Qx[L],i)<calc(Qx[L+1],i)))L++;
    return Qx[L];
}
//查询答案
il void DP(int *a,int *p,int n){
    L=0;R=-1;
    for(RG int i=1;i<=n;i++){
        if(i!=1)p[i]=ceil(calc(pop(i),i));Ql[L]=i;push(i);
    }
}
//DP序列
int main()
{
    n=read();
    for(RG int i=1;i<=n;i++)a[i]=read();
    DP(a,p,n);reverse(a+1,a+1+n);DP(a,s,n);
    for(RG int i=1;i<=n;i++)
        printf("%d\n",max(0,max(p[i],s[n-i+1])));
    //记得答案是非负整数
    return 0;
}

原文地址:https://www.cnblogs.com/cjfdf/p/8870067.html

时间: 2024-11-06 03:52:39

[POI2011]Lightning Conductor的相关文章

【BZOJ2216】[Poi2011]Lightning Conductor 决策单调性

[BZOJ2216][Poi2011]Lightning Conductor Description 已知一个长度为n的序列a1,a2,...,an.对于每个1<=i<=n,找到最小的非负整数p满足 对于任意的j, aj < = ai + p - sqrt(abs(i-j)) Input 第一行n,(1<=n<=500000)下面每行一个整数,其中第i行是ai.(0<=ai<=1000000000) Output n行,第i行表示对于i,得到的p Sample I

P3515 [POI2011]Lightning Conductor(决策单调性分治)

P3515 [POI2011]Lightning Conductor 式子可转化为:$p>=a_j-a_i+sqrt(i-j) (j<i)$ $j>i$的情况,把上式翻转即可得到 下面给一张图证明这是满足决策单调性的 把$a_j+sqrt(i-j)$表示在坐标系上 显然$sqrt(i-j)$的增长速度趋缓 曲线$a$被曲线$b$超过后是无法翻身的 对两个方向进行决策单调性分治,取$max$即可 #include<iostream> #include<cstdio>

BZOJ2216 : [Poi2011]Lightning Conductor

$f[i]=\max(a[j]+\lceil\sqrt{|i-j|}\rceil)$, 拆开绝对值,考虑j<i,则决策具有单调性,j>i同理, 所以可以用分治$O(nlogn)$解决. #include<cstdio> #include<cmath> #define N 500010 int n,i,l,r,mid,a[N],b[N],f[N],g[N]; inline void read(int&a){char c;while(!(((c=getchar())

bzoj2216: [Poi2011]Lightning Conductor(分治决策单调性优化)

每个pi要求 这个只需要正反DP(?)一次就行了,可以发现这个是有决策单调性的,用分治优化 #include<iostream> #include<cstring> #include<cstdlib> #include<cstdio> #include<cmath> #include<algorithm> using namespace std; const int maxn=500010,inf=1e9; int n; int a[

【bzoj2216】[Poi2011]Lightning Conductor 1D1D动态规划优化

Description 已知一个长度为n的序列a1,a2,…,an.对于每个1<=i<=n,找到最小的非负整数p满足 对于任意的j, aj < = ai + p – sqrt(abs(i-j)) Input 第一行n,(1<=n<=500000)下面每行一个整数,其中第i行是ai.(0<=ai<=1000000000) Output n行,第i行表示对于i,得到的p Sample Input 6532424 Sample Output 235354 题解 http

@bzoj - [email&#160;protected] [Poi2011]Lightning Conductor

目录 @[email protected] @[email protected] @part - [email protected] @part - [email protected] @part - [email protected] @accepted [email protected] @version - [email protected] @version - [email protected] @[email protected] @[email protected] 已知一个长度为

P3515 [POI2011]Lightning Conductor[决策单调性优化]

给定一序列,求对于每一个$a_i$的最小非负整数$p_i$,使得$\forall j \neq i $有$ p_i>=a_j-a_i+ \sqrt{|i-j|}$. 绝对值很烦 ,先分左右情况单独做.现在假设j都在i左边,则$p_i=max{a_j-a_i+ \sqrt{i-j}}=max{a_j+ \sqrt{i-j} }-a_i$.带根号,不易斜率优化,考虑证决策单调性. 假设最优决策为j,j之前的任意决策称之为$j'$,则有 $f[j]+\sqrt{i-j} \geqslant f[j']

P3515 [POI2011]Lightning Conductor

首先进行一步转化 $a_j \leq a_i + q - sqrt(abs(i - j))$ $a_i + q \geq a_j + sqrt(abs(i-j))$ 即 $q = max (a_j + sqrt(abs(i-j))) - a_i $ 我们对$i \geq j 和 j > i$ 分类讨论, 其实解决一种情况后将序列翻转再做一遍即可 有一种O($n^2$)的dp暴力应该不难想到 那么我们现在思考如何以比较优秀的时间复杂度解决 这里涉及到决策单调性 简单的说, 对于i来说, 它的答案来

BZOJ 2216 Lightning Conductor

决策单调.整体二分. #include<iostream> #include<cstdio> #include<cmath> #include<cstring> #include<algorithm> #define maxn 500500 using namespace std; int n,a[maxn],f[maxn],g[maxn]; void dp1(int left,int right,int l,int r) { if (left&