石子合并问题

一、问题分析

(一)问题描述：

在一个圆形操场的四周摆放着n堆石子。现要将石子有次序地合并成一堆。规定每次只能选相邻的2 堆石子合并成新的一堆，并将新的一堆石子数记为该次合并的得分。试设计一个算法，计算出将n堆石子合并成一堆的最小得分和最大得分。

(二)要求：

对于任意给定的n堆石子,计算合并成一堆的最小得分和最大得分。

二、问题的解决方案/算法选择/设计思路

环形区间DP本质上还是求解链形区间DP，只是一个有n个节点的环可以拆分成n个不同的链。如果分别计算每个链的结果，那么本题一定会超时。优化策略：将一个链复制一份并连接在其尾部，形成一个2n个节点的“长链”。上述n条不同的链都能在该“长链”中找到。枚举len的时候只需要枚举到len = n即可。

三、算法设计

(一)算法时间复杂度的分析

此问题的求解过程需要枚举区间长度，然后枚举起点和终点(因为长度确定，所以终点也就确定)，最后枚举断点，这样算下来是O(n^3)的时间复杂度。对于正常计算机来说1秒内可以进行10^7~10^8操作 ，所以n的极限应该是在400多。

(二)初始化

假设圆形石子长度为n，基于动态规划的问题求解,可以采用集合的方式求解，区间DP问题一般采用二维空间计算。对于问题需要求最大最小两个值我们可以开两个数组f(最大值),s(最小值)。因为需要将链再复制一次，所以数组开空间的时候要开到n*2，即f[n<<1][n<<1],s[n<<1][n<<1]。除此之外我们还要开两个数组，一个数组e用来存储我们输入每一堆石子的，即e[n<<1]，数组c用来存储前缀和,便于计算两堆石子合并时的花费。因为求最大值所以数组s里的每个值都初始化为0x3f3f3f3f(一个极大值，视情况而定)，由于石子合并的花费不会为负数，所以f数组不用管。

(三)状态转移方程计算(仅以f数组做例子)

状态转移方程的计算我们可以用到分集合的方式,对于f[i][j](j>=i)来讲，它代表的是合并从i到j这好几堆石子的最大花费。不管从i到j有多少堆石子，合并到最后一定是左边剩一堆，右边剩一堆，然后再合并成一堆。所以我们就需要枚举断点。所以我们的转态转移方程就有了f[i][j]=max(f[i][k]+f[k+1][j]+c[j]-c[i-1],f[i][j]),k是一层循环，从i到j-1。如果len==1，则f[i][i]=0。f[i][i]的状态是合理的，自己和自己合并的花费是0。如果len！=0则进行状态转移。

(四)结果输出

对于f求最大值，我们只要定义一个极小值a=0来保留最大值即可，对于这道题目我们取最小值的话应该写成a=max(a,f[i][i+n-1])。这样写的目的是我们要求的是一个圆形展开的链，长度是不应该变的，f[i][i+n-1]便是从i开始到i+n-1这n堆石子合并的最大花费。

#include<bits/stdc++.h>

using namespace std;

const int N=410;

int n;

int f[N][N],s[N][N];
int e[N],c[N];

int main(){
    
    scanf("%d",&n);
    
    for(int i=1;i<=n;++i) scanf("%d",&e[i]),e[i+n]=e[i];
    
    for(int i=1;i<=n<<1;++i) c[i]=c[i-1]+e[i];
    
    memset(s,0x3f,sizeof s);
    
    for(int len=1;len<=n;++len){
        for(int i=1;i+len-1<=n<<1;++i){
            int j=i+len-1;
            if(len==1) f[i][i]=s[i][i]=0;
            else{
                for(int k=i;k<j;++k){
                    f[i][j]=max(f[i][k]+f[k+1][j]+c[j]-c[i-1],f[i][j]);				
                    s[i][j]=min(s[i][k]+s[k+1][j]+c[j]-c[i-1],s[i][j]);
                }
            }
        }
    }
    
    int a=0x3f3f3f3f,b=0;
    
    for(int i=1;i<=n;++i){
        a=min(a,s[i][i+n-1]);
        b=max(b,f[i][i+n-1]);
    }
    
    printf("%d\n%d\n",a,b);
    
    return 0;
    
    
}

#include<bits/stdc++.h>
using namespace std;
const int N=2010;

int n;

int f[N][N],s[N][N];//f鏈€澶у€?锛宻鏈€灏忓€?
int e[N],c[N];
int d[N][N];//s鐨勬柇鐐?
 
int main(){
	
	scanf("%d",&n);	
	for(int i=1;i<=n;++i) scanf("%d",&e[i]),e[i+n]=e[i];	
	for(int i=1;i<=n<<1;++i) d[i][i]=i,c[i]=c[i-1]+e[i];
	
	memset(s,0x3f,sizeof s);
	
	for(int len=1;len<=n;++len){
		for(int i=1;i+len-1<=n<<1;++i){
			int j=i+len-1;
			if(len==1) f[i][j]=s[i][j]=0;
			else{
			    f[i][j]=max(f[i][j-1],f[i+1][j])+c[j]-c[i-1];
					
				for(int k=d[i][j-1];k<=d[i+1][j];k++){
					int tmp=s[i][k]+s[k+1][j]+c[j]-c[i-1];
							
					if(tmp<s[i][j]){
					s[i][j]=tmp;
					d[i][j]=k;
					}
				}	
			}
		}
	}
	
	int a=0x3f3f3f3f,b=0;
	for(int i=1;i<=n;++i){
		a=min(a,s[i][i+n-1]);
		b=max(b,f[i][i+n-1]);
	}
	printf("%d\n%d\n" ,a,b);
	
	return 0;
	 
}