石子合并(GarsiaWachs算法)

对于石子合并问题，有一个最好的算法，那就是GarsiaWachs算法。时间复杂度为O(n^2)。

它的步骤如下：

设序列是stone[]，从左往右，找一个满足stone[k-1] <= stone[k+1]的k，找到后合并stone[k]和stone[k-1]，再从当前位置开始向左找最大的j，使其满足stone[j] > stone[k]+stone[k-1]，插到j的后面就行。一直重复，直到只剩下一堆石子就可以了。在这个过程中，可以假设stone[-1]和stone[n]是正无穷的。

举个例子：

186 64 35 32 103

因为35<103，所以最小的k是3，我们先把35和32删除，得到他们的和67，并向前寻找一个第一个超过67的数，把67插入到他后面，得到：186 67 64 103,现在由5个数变为4个数了，继续：186 131 103，现在k=2（别忘了，设A[-1]和A[n]等于正无穷大）234 186，最后得到420。最后的答案呢？就是各次合并的重量之和，即420+234+131+67=852。

基本思想是通过树的最优性得到一个节点间深度的约束，之后证明操作一次之后的解可以和原来的解一一对应，并保证节点移动之后他所在的深度不会改变。具体实现这个算法需要一点技巧，精髓在于不停快速寻找最小的k，即维护一个“2-递减序列”朴素的实现的时间复杂度是O(n*n)，但可以用一个平衡树来优化，使得最终复杂度为O(nlogn)。

题目：http://poj.org/problem?id=1738

1. #include <iostream>  
2. #include <string.h>  
3. #include <stdio.h>  
4.   
5. using namespace std;  
6. const int N = 50005;  
7.   
8. int stone[N];  
9. int n,t,ans;  
10.   
11. void combine(int k)  
12. {  
13.     int tmp = stone[k] + stone[k-1];  
14.     ans += tmp;  
15.     for(int i=k;i<t-1;i++)  
16.         stone[i] = stone[i+1];  
17.     t--;  
18.     int j = 0;  
19.     for(j=k-1;j>0 && stone[j-1] < tmp;j--)  
20.         stone[j] = stone[j-1];  
21.     stone[j] = tmp;  
22.     while(j >= 2 && stone[j] >= stone[j-2])  
23.     {  
24.         int d = t - j;  
25.         combine(j-1);  
26.         j = t - d;  
27.     }  
28. }  
29.   
30. int main()  
31. {  
32.     while(scanf("%d",&n)!=EOF)  
33.     {  
34.         if(n == 0) break;  
35.         for(int i=0;i<n;i++)  
36.             scanf("%d",stone+i);  
37.         t = 1;  
38.         ans = 0;  
39.         for(int i=1;i<n;i++)  
40.         {  
41.             stone[t++] = stone[i];  
42.             while(t >= 3 && stone[t-3] <= stone[t-1])  
43.                 combine(t-2);  
44.         }  
45.         while(t > 1) combine(t-1);  
46.         printf("%d\n",ans);  
47.     }  
48.     return 0;  
49. }