# 【洛谷3515】[POI2011] Lightning Conductor（决策单调性）

### 代码

#include<bits/stdc++.h>
#define Tp template<typename Ty>
#define Ts template<typename Ty,typename... Ar>
#define Reg register
#define RI Reg int
#define Con const
#define CI Con int&
#define I inline
#define W while
#define N 500000
#define swap(x,y) (x^=y^=x^=y)
using namespace std;
int n,a[N+5];double s1[N+5],s2[N+5];

class FastIO
{
private:
#define FS 100000
#define pc(c) (C==E&&(clear(),0),*C++=c)
#define tn (x<<3)+(x<<1)
#define D isdigit(c=tc())
int T;char c,*A,*B,*C,*E,FI[FS],FO[FS],S[FS];
public:
I FastIO() {A=B=FI,C=FO,E=FO+FS;}
Tp I void read(Ty& x) {x=0;W(!D);W(x=tn+(c&15),D);}
Tp I void write(Ty x) {W(S[++T]=x%10+48,x/=10);W(T) pc(S[T--]);}
Tp I void writeln(Con Ty& x) {write(x),pc('\n');}
I void clear() {fwrite(FO,1,C-FO,stdout),C=FO;}
}F;
I void Solve(CI l,CI r,CI tl,CI tr,double *s)//分治处理决策单调性
{
if(l>r) return;RI mid=l+r>>1,p=tl;
for(RI i=tl;i<=tr;++i) if(s[mid]<a[i]-a[mid]+sqrt(mid-i)) s[mid]=a[p=i]-a[mid]+sqrt(mid-i);//找到最远能影响到mid的位置，同时更新mid的答案
Solve(l,mid-1,tl,p,s),Solve(mid+1,r,p,tr,s);//递归处理子区间
}
int main()
{