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‘]+\sqrt{i-j‘}$  ①

现要证$ f[j]+\sqrt{i-j+1} \geqslant f[j‘]+\sqrt{i-j‘+1}$ ②

即证$ \sqrt{i-j‘}+\sqrt{i-j+1}\geqslant \sqrt{i-j}+\sqrt{i-j‘+1}$(即想法求出一个不等式与①相加消项可得②)

即证$ \sqrt{i-j‘+1}-\sqrt{i-j‘}\geqslant \sqrt{i-j+1}+\sqrt{i-j}$(上式移项)

那么把它看成关于$j$的函数看单调性,设$g(j)=\sqrt{i-j+1}+\sqrt{i-j}$

对其求导,$g‘(j)=[(i-j+1)^{\frac{1}{2}}]‘ - [(i-j)^{\frac{1}{2}}]‘=-\frac{1}{2} (i-j+1)^{-\frac{1}{2}} + \frac{1}{2} (i-j)^{-\frac{1}{2}}=\frac{1}{2} (\frac{1}{\sqrt{i-j}}-\frac{1}{\sqrt{i-j+1}})$

由$ i-j<i-j+1$知$\frac{1}{2} (\frac{1}{\sqrt{i-j}}-\frac{1}{\sqrt{i-j+1}}) > 0$则函数$g(j)$单调增,则上不等式成立,满足单调性。

希望自己以后就按照这种方法想,虽然并卵。

证完决策单调性优化即可。



错误记录:第二次写的时候line37写成l<r了,,丢人。

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 typedef long long ll;
 4 typedef double db;
 5 template<typename T>inline char MIN(T&A,T B){return A>B?A=B,1:0;}
 6 template<typename T>inline char MAX(T&A,T B){return A<B?A=B,1:0;}
 7 template<typename T>inline T _min(T A,T B){return A<B?A:B;}
 8 template<typename T>inline T _max(T A,T B){return A>B?A:B;}
 9 template<typename T>inline T read(T&x){
10     x=0;int f=0;char c;while(!isdigit(c=getchar()))if(c==‘-‘)f=1;
11     while(isdigit(c))x=x*10+(c&15),c=getchar();return f?x=-x:x;
12 }
13 const int N=500000+7;
14 struct kochiya_sanae{
15     int l,r,pos;
16     kochiya_sanae(int l0=0,int r0=0,int pos0=0):l(l0),r(r0),pos(pos0){}
17 }q[N];
18 db sq[N],f[N],h[N];
19 int n,l,r;
20 inline void preprocess(){for(register int i=1;i<=n;++i)sq[i]=sqrt((db)i);}
21 inline db calc(int j,int i){return (db)h[j]+sq[i-j];}
22 inline int find_pos(int L,int R,int j,int i){
23     ++R;int mid;
24     while(L<R){
25         mid=L+R>>1;
26         if(calc(j,mid)<=calc(i,mid))R=mid;
27         else L=mid+1;
28     }
29     return R;
30 }
31 inline void dp(){
32     q[l=r=1]=kochiya_sanae(0,n,0);
33     for(register int i=1;i<=n;++i){
34         if(q[l].r<i)++l;else ++q[l].l;
35         MAX(f[i],calc(q[l].pos,i)-h[i]);
36         while(l<=r&&calc(q[r].pos,q[r].l)<=calc(i,q[r].l))--r;
37         if(r<l)q[r=l]=kochiya_sanae(i,n,i);
38         else{
39             int k;
40             if(calc(q[r].pos,q[r].r)>calc(i,q[r].r))k=q[r].r+1;
41             else k=find_pos(q[r].l,q[r].r,q[r].pos,i);
42             if(k<=n)q[r].r=k-1,q[++r]=kochiya_sanae(k,n,i);
43         }
44     }
45 }
46
47 int main(){//freopen("tmp.in","r",stdin);freopen("tmp.out","w",stdout);
48     read(n);for(register int i=1;i<=n;++i)read(h[i]);h[0]=-3;
49     preprocess();dp();reverse(h+1,h+n+1);reverse(f+1,f+n+1);dp();
50     for(register int i=n;i;--i)printf("%d\n",(int)ceil(f[i]));
51     return 0;
52 }

原文地址:https://www.cnblogs.com/saigyouji-yuyuko/p/10468742.html

时间: 2024-10-05 12:42:43

P3515 [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>

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来说, 它的答案来

【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

决策单调性优化dp

决策单调性: 对于一些dp方程,经过一系列的猜想和证明,可以得出,所有取的最优解的转移点(即决策点)位置是单调递增的. 即:假设f[i]=min(f[j]+b[j]) (j<i) 并且,对于任意f[i]的决策点g[i],总有f[i+1]的决策点g[i+1]>=g[i](或者<=g[i]) 那么,这个方程就具备决策单调性. 这个有什么用吗? 不懂具体优化方法的话确实也没有什么用.可能还是n^2的.只不过范围可能少了一些. 经典入门例题: Description: [POI2011]Ligh

决策单调性优化dp 专题练习

决策单调性优化dp 专题练习 优化方法总结 一.斜率优化 对于形如 \(dp[i]=dp[j]+(i-j)*(i-j)\)类型的转移方程,维护一个上凸包或者下凸包,找到切点快速求解 技法: 1.单调队列 : 在保证插入和查询的x坐标均具有单调性时可以使用 2.单调栈+二分:保证插入有单调性,不保证查询有单调性 3.分治+ 1 或 2:在每次分治时将\([l,mid]\)这段区间排序后插入,然后更新右区间\([mid+1,r]\)的答案 二.分治.单调队列维护有单调性的转移 (甚至还有分治套分治)

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[

Bzoj 1563: [NOI2009]诗人小G(决策单调性优化)

原题面 带有详细证明的转这里 题意:每一个线段有一个长度,有一个标准长,现在要把这些线段按照顺序分行,每行的不和谐值等于标准长和该行线段总长的差的绝对值的p次方.现在要求最小的不和谐值之和. 开始的时候完全读错题了,以为p==2 for ever.真是太天真.后来看数据范围才发现.我真是面向数据编程? n^2的dp是一眼秒的.然后如果是p=2,那就是一个非常像玩具装箱的斜率优化dp.对于4.5的数据范围.可以想到用贪心优化dp.因为每一行的长度不会通过拼接线段(线段条数>=2)达到2*标准长,这

题解——[NOI2009]诗人小G 决策单调性优化DP

第一次写这种二分来优化决策单调性的问题.... 调了好久,,,各种细节问题 显然有DP方程: f[i]=min(f[j] + qpow(abs(sum[i] - sum[j] - L - 1))); 其中f[i]代表到了第i个句子的最小答案 qpow用于处理^p sum为前缀和 (同时为了处理句子之间的空格问题,我们在统计前缀和的时候就默认在句子后面加一个空格, 然后在计算的时候,由于每一行只有最后一个不用加空格,直接减掉这个多加的空格即可获得正确长度) 首先我们可以打表发现是满足决策单调性的,