51 nod 1023 石子归并 V3(GarsiaWachs算法)

1023 石子归并 V3
基准时间限制：2 秒 空间限制：131072 KB 分值: 320 难度：7级算法题

N堆石子摆成一条线。现要将石子有次序地合并成一堆。规定每次只能选相邻的2堆石子合并成新的一堆，并将新的一堆石子数记为该次合并的代价。计算将N堆石子合并成一堆的最小代价。

 

例如： 1 2 3 4，有不少合并方法

1 2 3 4 => 3 3 4(3) => 6 4(9) => 10(19)

1 2 3 4 => 1 5 4(5) => 1 9(14) => 10(24)

1 2 3 4 => 1 2 7(7) => 3 7(10) => 10(20)

 

括号里面为总代价可以看出，第一种方法的代价最低，现在给出n堆石子的数量，计算最小合并代价。

 

Input

第1行：N（2 <= N <= 50000)
第2 - N + 1：N堆石子的数量（1 <= A[i] <= 10000）

Output

输出最小合并代价

Input示例

4
1
2
3
4

Output示例

19

/*
51 nod 1023 石子归并 V3(GarsiaWachs算法)

problem:
给你n个石碓,相邻两个可以合并代价是它们的和. 求总体的最小代价

参考:
http://blog.sina.com.cn/s/blog_a825ada90101no1m.html
这博客解释不错,但仍不是很懂这个算法的原理 T T.

hhh-2016/09/05-21:14:18
*/
#pragma comment(linker,"/STACK:124000000,124000000")
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <vector>
#include <math.h>
#include <queue>
#include <set>
#include <map>
#define lson  i<<1
#define rson  i<<1|1
#define ll long long
#define clr(a,b) memset(a,b,sizeof(a))
#define scanfi(a) scanf("%d",&a)
#define scanfs(a) scanf("%s",a)
#define scanfl(a) scanf("%I64d",&a)
#define scanfd(a) scanf("%lf",&a)
#define key_val ch[ch[root][1]][0]
#define eps 1e-7
#define inf 0x3f3f3f3f3f3f3f3f
using namespace std;
const ll mod = 1000000007;
const int maxn = 50050;
const double PI = acos(-1.0);

ll t[maxn];
ll ans;
int num;
void dfs(int cur)
{
    ll tval = t[cur-1] + t[cur];
    ans = ans+(ll)tval;
    for(int i = cur;i < num-1;i++)
        t[i] = t[i+1];
    int k;
    num --;
    for(k = cur-1;k >= 1 && t[k-1] < tval;k--)
    {
        t[k] = t[k-1];
    }
    t[k] = tval;

    while(k >= 2 && t[k] >= t[k-2])
    {
        int len = num-k;
        dfs(k-1);
        k = num - len;

    }
}

int main(){
    int n;
    while(scanfi(n)!=EOF){
//    clr(t,inf);
    for(int i =0 ;i < n;i++)
        scanfi(t[i]);
    ans = 0;
    num = 1;
    for(int i =1;i < n;i++)
    {
        t[num ++ ] = t[i];
        while(num >= 3 &&  t[num-3] <= t[num-1])
            dfs(num-2);
    }
    while(num > 1) dfs(num-1);
    printf("%I64d\n",ans);
    }
    return 0;
}