【NOI2009诗人小G】动态规划、1D1D经典优化

本题其实是个大裸题

f[i]=max{f[j]+w(j,i)}的方程基本上是摆明了的，这个方程显然是O(n^2)，虽然有30s的时限，但也是过不了的

优化，势在必行

这个方程最麻烦的一点就在于w(j,i)的形式过于飘逸，连p次方都来了，以至于诸如斜率优化，单调队列等一般方法在这里用不上

那怎么办？只剩神法决策单调性了，这个方程有决策单调性吗？怎么证？ p.s 决策单调性就是对于i<j,i的最优决策<=j的最优决策

本人数学不才，但可以隐约感觉到应该是有的，可以这样想，如果没有，岂不没什么方法了？

上句话是废话，你可以略过

证明过程byvoid博客讲的非常清楚，此处略过（其实是自己不会）

光有结论还不够，具体怎么做还得学

这里我就再偷一懒，网上有一篇论文讲得非常好，名字好像是1D1D优化初步什么的，我就不讲了

大概就是每次二分一个决策的“优势”区间，然后用个双端队列维护决策集合就可以了

写的有点丑，比盾哥和mt都慢了1、2倍，好像是int64和extended慢了，反正感觉写的不怎么漂亮，但还是贴出来吧

program syj; {dp} type real=extended; var n,len,p,task,i,st,ed,ll,rr,mid:longint; ok:boolean; t:real; l,r,q:array[0..100005]of longint; ss:ansistring; f,s:array[0..100005]of int64; function calc(j,i:longint):real; begin if s[i]-s[j]-1-len=0 then calc:=f[j] else calc:=f[j]+exp(p*ln(abs(s[i]-s[j]-1-len))); end; begin assign(input,'poet.in');reset(input); assign(output,'poet.out');rewrite(output); readln(task); for task:=1 to task do begin readln(n,len,p); for i:=1 to n do begin readln(ss);s[i]:=length(ss); inc(s[i],s[i-1]); end; for i:=1 to n do s[i]:=s[i]+i; st:=1;ed:=1;q[1]:=0;l[0]:=1;r[0]:=n; ok:=true; for i:=1 to n do begin if r[q[st]]<i then inc(st); t:=calc(q[st],i); if t>1e18 then begin t:=1.5e18; end; f[i]:=round(t); l[q[st]]:=i+1; if l[q[st]]>r[q[st]] then inc(st); while (st<=ed)and(calc(i,l[q[ed]])<calc(q[ed],l[q[ed]])+1e-8) do dec(ed); inc(ed);q[ed]:=i; if st<ed then begin ll:=l[q[ed-1]]+1;rr:=n; while ll<rr do begin mid:=(ll+rr) div 2; if calc(i,mid)<calc(q[ed-1],mid)+1e-8 then rr:=mid else ll:=mid+1; end; if calc(i,ll)<calc(q[ed-1],ll)+1e-8 then begin l[i]:=ll;r[i]:=n;r[q[ed-1]]:=ll-1; end else dec(ed); end else begin l[i]:=i+1;r[i]:=n; end; end; if f[n]>1e18 then writeln('Too hard to arrange') else writeln(f[n]); end; close(input); close(output); end.