石子归并

http://blog.csdn.net/abcjennifer/article/details/5805330

1.链式归并

 

问题描述

设有N堆沙子排成一排，其编号为1,2,3,…,N(N<=100)。每堆沙子有一定的数量。现要将N堆沙子并成为一堆。归并的过程只能每次将相邻的两堆沙子堆成一堆，这样经过N-1次归并后成为一堆。找出一种合理的归并方法，使总的代价最小。

【输入格式】

    输入由若干行组成，第一行有一个整数，n（1≤n≤100）；表示沙子堆数。第2至m+1行是每堆沙子的数量。 

【输出格式】

    一个整数，归并的最小代价。

【输入样例】

输入文件名：shizi.in

7
13
7
8
16
21
4
18

【输出样例】

输出文件名：shizi.out

239

---

令f[i,j]表示归并第i个数到第j数的最小代价,sum[i,j]表示第i个数到第j个数的和，这个可以事先计算出来。sum[i,j]可以在O(1)的时间内算出.
容易的到以下的动态转移方程：

阶段：以归并石子的长度为阶段，一共有n-1个阶段。
状态：每个阶段有多少堆石子要归并，当归并长度为2时，有n-1个状态；
当归并长度为3时，有n-2个状态；
当归并长度为n时，有1个状态。
决策：当归并长度为2时，有1个决策；当归并长度为3时，有2个决策；
当归并长度为n时，有n-1个决策。

 

#include <iostream>   
using namespace std;
#define M 101
#define INF 1000000000
int n,f[M][M],sum[M][M],stone[M];
int main()
{
	int i,j,k,t;
	cin>>n;
	for(i=1;i<=n;i++)
		scanf("%d",&stone[i]);

	for(i=1;i<=n;i++)
	{
		f[i][i]=0;
		sum[i][i]=stone[i];
		for(j=i+1;j<=n;j++)
			sum[i][j]=sum[i][j-1]+stone[j];
	}

	for(int len=2;len<=n;len++)//归并的石子长度
	{
		for(i=1;i<=n-len+1;i++)//i为起点，j为终点
		{
			j=i+len-1;
			f[i][j]=INF;
			for(k=i;k<=j-1;k++)
			{
				if(f[i][j]>f[i][k]+f[k+1][j]+sum[i][j])
					f[i][j]=f[i][k]+f[k+1][j]+sum[i][j];
			}
		}
	}
	printf("%d/n",f[1][n]);  
	return 0;
}

 

[算法优化]

 见 http://it.dgzx.net/drkt/oszt/zltk/yxlw/dongtai3.htm

 

优化后，无非加了一个s[i,j] 记录i...j之间f[i,j]为最优状态的分割点

#include <iostream>   
using namespace std;
#define M 101
#define INF 1000000000
int n,f[M][M],sum[M][M],stone[M],s[M][M];
int main()
{
    int i,j,k,t;
    cin>>n;
    for(i=1;i<=n;i++)
        scanf("%d",&stone[i]);

    for(i=1;i<=n;i++)
    {
        f[i][i]=0;
        s[i][i]=i;
        sum[i][i]=stone[i];
        for(j=i+1;j<=n;j++)
            sum[i][j]=sum[i][j-1]+stone[j];
    }

    for(int len=2;len<=n;len++)//归并的石子长度
    {
        for(i=1;i<=n-len+1;i++)//i为起点，j为终点
        {
            j=i+len-1;
            f[i][j]=INF;
            for(k=s[i][j-1];k<=s[i+1][j];k++)
            {
                if(f[i][j]>f[i][k]+f[k+1][j]+sum[i][j])
                {
                    f[i][j]=f[i][k]+f[k+1][j]+sum[i][j];
                    s[i][j]=k;
                }
            }
        }
    }
    printf("%d/n",f[1][n]);
    return 0;
}

2.

环型石子合并

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

输入文件
输入文件stone.in包含两行，第1行是正整数n(1≤n≤100)，表示有n堆石子。第2行有n个整数，分别表示每堆石子的个数。

输出文件
输出文件stone.out 包含两行，第1行中的数是最小得分；第2行中的数是最大得分。

输入样例
4
4 4 5 9

输出样例
43
54

#include<cstdio>
using namespace std;

typedef int Arr[210][210];
int N;
Arr M1,M2,S1,S2;
int main()
{
    scanf("%d",&N);
    for(int i=1;i<=N;++i)
	{
        int a; scanf("%d",&a);
        M1[i][i]=M2[i][i]=M1[i+N][i+N]=M2[i+N][i+N]=a;
        S1[i][i]=S2[i][i]=S1[i+N][i+N]=S2[i+N][i+N]=0;
    }
    for(int k=2;k<=N;++k)
	{
        for(int i=1;i<=2*N-k+1;++i)
		{
            int j=i+k-1;
            M1[i][j]=S1[i][j]=0xFFFFFFF;
            M2[i][j]=S2[i][j]=0;
            for(int t=j-1;t>=i;--t)
			{
                if(M1[i][j]>M1[i][t]+M1[t+1][j])
				{
                    M1[i][j]=M1[i][t]+M1[t+1][j];
                    if(S1[i][j]>M1[i][j]+S1[i][t]+S1[t+1][j])
                        S1[i][j]=M1[i][j]+S1[i][t]+S1[t+1][j];
                }
                else if(M1[i][j]==M1[i][t]+M1[t+1][j])
				{
                    if(S1[i][j]>M1[i][j]+S1[i][t]+S1[t+1][j])
                        S1[i][j]=M1[i][j]+S1[i][t]+S1[t+1][j];
                }
                if(M2[i][j]<=M2[i][t]+M2[t+1][j])
				{
                    M2[i][j]=M2[i][t]+M2[t+1][j];
                    if(S2[i][j]<M2[i][j]+S2[i][t]+S2[t+1][j])
                        S2[i][j]=M2[i][j]+S2[i][t]+S2[t+1][j];
                }
                else if(M2[i][j]==M2[i][t]+M2[t+1][j])
				{
                    if(S2[i][j]<M2[i][j]+S2[i][t]+S2[t+1][j])
                        S2[i][j]=M2[i][j]+S2[i][t]+S2[t+1][j];
                }
            }
        }
    }
    int ans=0xFFFFFFF;
    for(int i=1;i<=N;++i)
        if(S1[i][i+N-1]<ans)
            ans=S1[i][i+N-1];
    printf("%d/n",ans);
    ans=0;
    for(int i=1;i<=N;++i)
        if(S2[i][i+N-1]>ans)
            ans=S2[i][i+N-1];
    printf("%d/n",ans);
    return 0;
}