【NOI2009】诗人小G【决策单调性dp】
傳送門
設(shè)句子加上空格長度前綴和為sss
轉(zhuǎn)移方程為
fi=min?1≤j<i{fj+∣si?sj?L?1∣P}f_i=\min_{1\leq j<i}\{f_j+|s_i-s_j-L-1|^P\}fi?=1≤j<imin?{fj?+∣si??sj??L?1∣P}
不難發(fā)(cai)現(xiàn)(chu)決策具有單調(diào)性
本文重點(diǎn)講決策單調(diào)dp的實(shí)現(xiàn)
由于決策單調(diào),所以每個(gè)決策有效的都是一個(gè)區(qū)間
用隊(duì)列按順序維護(hù)所有區(qū)間
開始時(shí)加入“決策000的生效區(qū)間為整個(gè)序列”
每到一個(gè)位置,彈出所有已經(jīng)失效的區(qū)間,隊(duì)頭就是當(dāng)前最優(yōu)決策,算出答案
然后用當(dāng)前位置更新隊(duì)列
先將隊(duì)尾無用的彈掉,因?yàn)閱握{(diào),所以判斷如果隊(duì)尾決策對(duì)區(qū)間最左邊沒有當(dāng)前位置作為決策優(yōu),彈隊(duì)尾。
如果隊(duì)列中仍有元素,二分出當(dāng)前砍掉的位置。
最后把挖掉的部分都用當(dāng)前位置填上。注意如果沒挖直接忽略
本題中由于主要是大小比較,可以犧牲精度追求值域,用long double存
#include <iostream> #include <cstdio> #include <cstring> #include <cctype> #include <cmath> #define MAXN 100005 using namespace std; typedef long double ld; struct key{int c,p;}q[MAXN]; int head,tail; char str[MAXN][35]; int s[MAXN],n,l,p,ans[MAXN]; ld dp[MAXN]; inline ld abso(const ld& x){return x<0? -x:x;} inline ld qpow(ld a,int p) {ld ans=1;while (p){if (p&1) ans=ans*a;a=a*a;p>>=1;}return ans; } inline ld calc(int n,int i){return dp[i]+qpow(abso(s[n]-s[i]-l-1),p);} void print(int n) {if (!n) return;print(ans[n]);for (int i=ans[n]+1;i<=n;i++){printf("%s",str[i]);if (i<n) putchar(' ');}puts(""); } int main() {int T;scanf("%d",&T);while (T--){head=1,tail=0;scanf("%d%d%d",&n,&l,&p);for (int i=1;i<=n;i++) scanf("%s",str[i]),s[i]=s[i-1]+strlen(str[i])+1;q[++tail]=(key){0,n};for (int i=1;i<=n;i++){while (head<tail&&i>q[head].p) ++head;dp[i]=calc(i,ans[i]=q[head].c);while (head<=tail&&calc(q[tail-1].p+1,i)<calc(q[tail-1].p+1,q[tail].c)) --tail;if (head<=tail){int l=q[tail-1].p+1,r=q[tail].p,mid;while (l<r){mid=(l+r+1)>>1;if (calc(mid,q[tail].c)<calc(mid,i)) l=mid;else r=mid-1;}q[tail].p=l;}if (q[tail].p<n) q[++tail]=(key){i,n};}if (dp[n]>1e18) puts("Too hard to arrange");else{printf("%lld\n",(long long)dp[n]);print(n);}puts("--------------------");}return 0; }總結(jié)
以上是生活随笔為你收集整理的【NOI2009】诗人小G【决策单调性dp】的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 风云决的歌词 风云决歌简介
- 下一篇: 【NOI2012】迷失游乐园【概率期望】