HYSBZ 1010 玩具装箱toy （决策单调DP）

题意：

　　有n个玩具，要将它们分为若干组，玩具长度C可能不同。给出n个玩具的摆放顺序，连续的任意多个玩具都可以成为一组。区间[i,j]成为一组的费用是cost=(j-i+Sigma(C_k)-L)²且i<=k<=j。给定n和L和每个玩具的长度，问分组后费用总和是多少? (n<=5*10⁴)。

思路：

　　转移方程：dp[i]=min( dp[j]+(sum[i]-sum[j]+i-j+1-L)²  )。sum[i]表示前i件玩具长度的总和，0<j<i，（i-j+1）表示与i同组的玩具个数。

　　根据方程是可以推出这题是满足决策单调性的。以下是抄来的证明，稍微修改：

　　令f[i]=sum[i]+i, c=1+L，则dp[i]=min( dp[j]+(f[i]-f[j]-c)²  )

　　1.证明决策单调性

　　假设在状态i处的k决策优于j决策，且j<k，那么 dp[k]+(f[i]-f[k]-c)²<=dp[j]+(f[i]-dp[j]-c)²。

　　而对于i后面的某个状态t，设f[t]=f[i]+v，先不管v是多少。

　　要证明：dp[k]+(f[t]-f[k]-c)²<=dp[j]+(f[t]-f[j]-c)²

　　只要证(将f[t]=f[i]+v代入)：dp[k]+(f[i]+v-f[k]-c)²<=dp[j]+(f[i]+v-f[j]-c)²

　　只要证dp[k]+(f[i]-f[k]-c)²+2v*(f[i]-f[k]-c)+v²  <=  dp[j]+(f[i]-f[j]-c)²+2v*(f[i]-f[j]-c)+v²。

　　由于假设，所以只要证： 2v*(f[i]-f[k]-c)<=2v*(f[i]-f[j]-c)。

　　即证：f[k]>=f[j]（显然）

　　证明完毕

　　思路很明确，一直卡在二分上面，噗。

　　用一个队列来维护这些区间段，由于区间段必定是连在一起的，所以只需要记录左端点L以及更新这个区间的决策k。如果队列为空，则后面全部由i来更新得到，若非空，那么判断队尾的L，是否由i来更新会更优，若是，则pop掉队尾，继续同样的动作，直到队列为空或者i作为决策不如队尾的L更好，那么i可以更新的就是[L,n]之中的尾部区间[r,n]，而r可以用二分查找的方式。细节上很容易写挫，比如i决策可能完全都可以用武之地，不用二分去找了，否则会错；二分时必定要保证r由i来更新更佳，且有可能会出现等于的情况。复杂度O(nlogn)，斜率优化等再写。

 //#include <bits/stdc++.h>
 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <cmath>
 #include <set>
 #include <deque>
 #include <map>
 #include <algorithm>
 #include <vector>
 #include <iostream>
 #define pii pair<int,int>
 #define back que[rear-1]
 #define INF 0x7f7f7f7f
 #define LL long long
 #define ULL unsigned long long
 using namespace std;
 const double PI  = acos(-1.0);
 const int N=;
 
 LL len[N], dp[N], L;
 int q[N], d[N], n, l, r;   //区间以及决策
 LL cost(int j,int i)        //用j来更新i的费用
 {
     return  dp[j]+(len[i]-len[j]-L)*(len[i]-len[j]-L);
 }
 
 int find(int i,int k,int st)
 {
     int ll=st, rr=n;
     while(ll<rr)
     {
         int mid=rr-(rr-ll+)/;
         if( cost(i,mid)<cost(k,mid)) rr=mid;
         else    ll=mid+;
     }
     return rr;
 }
 LL cal()
 {
     l=r=;
     d[]=;q[]=;  //初始时，0可以更新[1,n]
     for(int i=; i<=n; i++)
     {
         dp[i]=cost(d[l], q[l]++);   //q[l]永远等于i
         if( l<r && q[l]==q[l+] )  l++;
 
         while( l<=r && cost(i,q[r])<cost(d[r],q[r]) )   r--;
         if(l>r) //只能用i来更新
         {
             q[++r]=i+;
             d[r]=i;
         }
         else if( cost(i,n)<cost(d[r],n))
         {
             int tmp=find(i, d[r], q[r]);
             q[++r]=tmp;
             d[r]=i;
         }
     }
     return dp[n];
 }
 
 int main()
 {
     //freopen("input.txt","r",stdin);
     while(~scanf("%d%lld",&n,&L))
     {
         L++;len[]=;
         for(int i=; i<=n; i++)
         {
             scanf("%lld",&len[i]);
             len[i]+=len[i-];
         }
         for(int i=; i<=n; i++) len[i]+=i;
         printf("%lld\n", cal() );
     }
     return ;
 }

AC代码