“玲珑杯”ACM比赛 Round #21-A- 西方制度 （区间DP+平行四边形优化）

A -- 西方制度

Time Limit：2s Memory Limit：256MByte

Submissions：162Solved：44

DESCRIPTION

法力浮·鳝AK为了开启异世界的大门，踏上了收集五行之灵。众所周知，五行之灵分布在世界的五个方向。因此，在成功获得了土之灵之后，鳝前往西方寻找金之灵。
经过重重难关，鳝终于到达了亚瑟王（Arthur）的安息之地——阿瓦隆（Avalon）。然而，大贤者梅林（Merlin）告知鳝：想要得到金之灵必须说出亚瑟王密码，而密码正是不列颠王国制度的编码，具体是这样的：
有n条制度，每条制度有ai个单位长。给每一条制度一个唯一的0-1串，并且这种0-1串可以保证任何一种0-1串都不是别的0-1串的前缀。（否则负责审判的法官会疯掉的）
并且，为了减少制法经费以便给亚瑟买吃的建设王国，所以她定义了一个码长L，L=a1编码长度a1 + a2编码长度a2 + …… + an的编码长度*an。
例如现在有五条制度ABCDE，分别有1,3,4,2,5长，此时给它们一个编号：
A:000,B:001,C:01,D:10,E:11，此时对应的码长就是L=34。
这种编码方式显然太naïve显示不出不列颠王国的伟大。因此亚瑟有定义了一种新的规则，就是在满足上述规则的情况下，必须使编码升序排列。
比如存在一种方案:
A:000 B:01 C:10 D:001 E:11
但是不满足升序，所以就是错误的。
我们姑且把这种编码方案称为亚瑟编码。
那么现在的问题就是，请你对给定的制度序列，找到一个最短的合法的编码长度。

题意：N组元素，每次合并相邻的两组元素，合并代价是两组元素个数之和，求合并代价最小。

题解：

设dp[i][j]表示第i到第j合并的最优值，sum[i][j]表示第i到第j的总数量。那么就有状态转移公式：

1、dp[i][j]=0 (i==j)
2、dp[i][j]=min(dp[i][k]+dp[k][j])+sum[i][j] (i!=j)
此时算法复杂为O(n^3)，这里可以利用平行四边形优化降为O(n^2)：
由上面的方程式可知我们每次求dp[i][j]的关键是找到合适的k值，设p[i][j]为dp[i][j]的这个合适的k值，根据平行四边形规则有以下不等式：p[i][j-1]<=p[i][j]<=p[i+1][j]。
那么求解dp[i]i+L的复杂度就为：
(p[2,L+1]-p[1,L])+(p[3,L+2]-p[2,L+1])…+(p[n-L+1,n]-p[n-L,n-1])=p[n-L+1,n]-p[1,L]≤n。
复杂度为O(n)。然后L从1循环至n，总复杂度就为O(n^2)。

这里有关于四边形优化的证明：http://www.cnblogs.com/jiu0821/p/4493497.html

#include<set>  
#include<map>     
#include<stack>            
#include<queue>            
#include<vector>    
#include<string> 
#include<time.h>
#include<math.h>            
#include<stdio.h>            
#include<iostream>            
#include<string.h>            
#include<stdlib.h>    
#include<algorithm>   
#include<functional>    
using namespace std;            
#define ll long long       
#define inf  1000000000       
#define mod 1000000007             
#define maxn 1005
#define lowbit(x) (x&-x)            
#define eps 1e-9
const int INF = 0x3f3f3f3f;
int n, a[maxn];
int p[maxn][maxn];
ll dp[maxn][maxn], sum[maxn];
int main(void) 
{
    while(scanf("%d", &n)!=EOF) 
	{
        for(int i = 1; i <= n; i++)
			scanf("%d", &a[i]);
        memset(dp, 0x3f, sizeof(dp));
        sum[0] = 0;
        for(int i = 1; i <= n; i++)
		{
            dp[i][i] = 0;
            p[i][i] = i;
            sum[i] = sum[i - 1] + a[i];
        }
        for(int len = 2; len <= n; len++) 
		{
            for(int i = 1; i + len - 1 <= n; i++)
			{
                int j = i + len - 1;
                for(int k = p[i][j - 1]; k <= p[i + 1][j]; k++)
				{
                    ll tmp = dp[i][k] + dp[k + 1][j] + (sum[j] - sum[i - 1]);
                    if(tmp < dp[i][j])
					{
                        dp[i][j]= tmp;
                        p[i][j] = k;
                    }
                }
            }
        }
        printf("%lld\n", dp[1][n]);
    }
    return 0;
}