[区间DP]石子合并

题目：

在一个圆形操场的四周摆放n堆石子(n≤500),现要将石子有次序地合并成一堆。规定每次只能选相邻的两堆合并成新的一堆，并将新的一堆的石子数，记为该次合并的得分。

编一程序，由文件读入堆数n及每堆的石子数；

⑴ 选择一种合并石子的方案，使得做n-1次合并，得分的总和最小；

⑵ 选择一种合并石子的方案，使得做n-1次合并，得分的总和最大。

输入文件：

第1行为n,代表n堆石子；

第2行为n个数，代表每堆石子的数量。

输出文件：

第1行为合并的得分总和最小值；

第2行为合并的得分总和最大值。

输入样例：

4

4 5 9 4(从最上面的一堆数起，顺时针数)

输出样例：

43(最小)

54(最大)

思路：

区间DP，状态转移方程：fmin[i][j] = min{f[i][k] + f[(i+k)%n][j-k] + sum[i][j]};

           fmax[i][j] = max{f[i][k] + f[(i+k)%n][j-k] + sum[i][j]};  

注意石子合并是成环的，所以 (i+k)%n;

以及注意初始化；

代码：

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <string.h>
#include <map>
#include <vector>
using namespace std;

typedef map<int, int> mii;
typedef vector<int> vi;
typedef pair<int, int> pii;
//f[i][j] = min{f[i][k] + f[(i+k)%n][j-k] + sum[i][j]};
const int N = 510;
#define inf 1<<29
int n;
int ans;
int a[N];//原始数据
int fMin[N][N];//第i堆起，顺时针合并j堆的总得分min
int fMax[N][N];//第i堆起，顺时针合并j堆的总得分max
/*
fmin与fmax出现error:'int fmin/fmax[510][510]' redeclared as different kind of symbol
原因：定义和声明与系统库冲突
*/
int sum[N][N];//代表这次合并的得分

void dp()
{
    int i, j, k;
    int t = 0;

    for(j = 2; j <= n; j++)//合并j堆
        for(i = 0; i < n; i++)//第i堆起
        {
            fMin[i][j] = inf;
            fMax[i][j] = 0;
            for(k = 1; k <= j-1; k++)
            {
                t = (i+k) % n;
                //就是此处出了问题！原来是t = (i+k-1) % n;
                if(fMin[i][k] + fMin[t][j-k] + sum[i][j] < fMin[i][j])
                    fMin[i][j] = fMin[i][k] + fMin[t][j-k] + sum[i][j];
                if(fMax[i][k] + fMax[t][j-k] + sum[i][j] > fMax[i][j])
                    fMax[i][j] = fMax[i][k] + fMax[t][j-k] + sum[i][j];
            }
        }
}
int main()
{
    int i, j;

    scanf("%d", &n);
    for(i = 0; i < n; i++)
    {
        scanf("%d", &a[i]);
        sum[i][1] = a[i];//初始化，第i堆起，合并一堆（即其本身）的总得分
        fMin[i][1] = 0;//初始化
        fMax[i][1] = 0;//初始化
    }
/*
若n: 1 to n 且 t=(i+1)%(n+1) 初始化有问题（如 sum[4][2] = a[4] + sum[0][1]）
所以，n: 0 to (n-1)
*/
    for(j = 2; j <= n; j++)
        for(i = 0; i < n; i++)
        {
            int t = (i+1) % n;
            sum[i][j] = a[i] + sum[t][j-1];
        }
    dp();
    int Min = inf;
    int Max = 0;
    for(i = 0; i < n; i++)
    {
        if(fMin[i][n] < Min)
            Min = fMin[i][n];
        if(fMax[i][n] > Max)
            Max = fMax[i][n];
    }
    printf("%d\n%d", Min, Max);

    return 0;
}
/*
4
4 5 9 4
------------
43
54
*/

反思：

我发现对我来说确定初始化很难，哪些状态应该提前确定下来？

经过练习，我认为那些“有确切数值的”“不会再改变的”状态应该被提前初始化；

还有一种，就是“会再改变的”，但“需要有数值”的，应该被提前初始化，但是这个度很难把握，需要多做题积累经验！