矩阵链相乘问题

参考文章

5
30 35 15 5 10 20

输出样例:
在这里给出相应的输出。例如：

11875

AC代码：

#include <iostream>
#include <cstdio>
using namespace std;
const int maxn = 10000;
int p[maxn], n;
int m[maxn][maxn], s[maxn][maxn]; // 记录最小乘法运算次数,和最后一次运算的位置
int r; // r为链长
int MatrixChain(int P[], int n)
{
    for(int i = 1; i <= n; i++)
    {
        for(int j = 1; j <= n; j++)
        {
            m[i][j] = 0;
            s[i][j] = i;
        }
    }
    for(int r = 2; r <= n; r++)
    {
        for(int i = 1; i <= n - r + 1; i++) // 前边界为i
        {
            int j = i + r - 1; // 前边界i链长为r的后边界为j
            m[i][j] = m[i + 1][j] + p[i-1] * p[i] * p[j];
            s[i][j] = i; //记录分隔位置;
            for(int k = i + 1; k < j; k++)
            {
                int t = m[i][k] + m[k + 1][j] + p[i-1] * p[k] * p[j];
                if(t < m[i][j])
                {
                    m[i][j] = t;
                    s[i][j] = k;
                }
            }
        }
    }
    return m[1][n];
}
int main()
{
    cin >> n; //n个矩阵
    for(int  i = 0; i <= n; i++)
    {
        cin >> p[i];
    }
    cout << MatrixChain(p, n);
    return 0;
}

拓展代码：（递归法 动态规划迭代法 输出括号分隔形式 单输出分割点）

#include <iostream>
#include <cstdio>
#include <vector>
using namespace std;
const int maxn = 10000;
int p[maxn], n;
int m[maxn][maxn], s[maxn][maxn]; // 记录最小乘法运算次数,和最后一次运算的位置
int r, k; // r为链长
vector<int> result;
int MatrixChain(int P[], int n) //迭代法
{
    for(int i = 1; i <= n; i++)
    {
        for(int j = 1; j <= n; j++)
        {
            m[i][j] = 0;
            s[i][j] = i;
        }
    }
    for(int r = 2; r <= n; r++)
    {
        for(int i = 1; i <= n - r + 1; i++) // 前边界为i n - r + 1为最后一个以r为规模的矩阵链的前边界
        {
            int j = i + r - 1; // 前边界i链长为r的后边界为j
            m[i][j] = m[i + 1][j] + p[i-1] * p[i] * p[j]; //A1(A2A3)
            s[i][j] = i; //记录分隔位置;
            for(int k = i + 1; k < j; k++) //i<=k<j
            {
                int t = m[i][k] + m[k + 1][j] + p[i-1] * p[k] * p[j]; // 转移状态方程
                if(t < m[i][j])
                {
                    m[i][j] = t;
                    s[i][j] = k;
                }
            }
        }
    }
    return m[1][n];
}
int RecurMatrixChain(int p[], int i, int j) //递归
{
    if(i == j) //递归边界
    {
        m[i][j] = 0;
        s[i][j] = i;
        return m[i][j];
    }
    for(int k = i; k < j; k++) //i <= k < j
    {
        int q = RecurMatrixChain(p, i, k) + RecurMatrixChain(p, k + 1, j) + p[i-1] * p[k] * p[j];
        if(q < m[i][j])
        {
            m[i][j] = q;
            s[i][j] = k;
        }
//         m[i][j] = min(m[i][j], RecurMatrixChain(p, i, k) + RecurMatrixChain(p, k + 1, j) + p[i-1] * p[k] * p[j]); //如果不记录分割点直接找最小就好
    }
    return m[i][j];
}
void print(int i, int j) //输出带括号的分隔
{
	if(i == j)
		printf("A%d",i);
	else{
		printf("(");
		print(i, s[i][j]);
		print(s[i][j]+1, j);
		printf(")");
	}
}
vector<int> find() //记录保存分割点 开容器 直接追加至容器内做记录
{
    int i = 1;
    int j = n;
    while(1)
    {
        int t = s[i][j];
        result.push_back(t);
        j = t;
        if(i == j) break;
    }
    return result;
}
int main()
{
    cin >> n; //n个矩阵
    for(int  i = 0; i <= n; i++)
    {
        cin >> p[i];
    }
    cout << MatrixChain(p, n) << endl;  //迭代法
    cout << RecurMatrixChain(p, 1, n) << endl; //递归

    print(1, n); //输出
    cout << endl;
    vector<int> res = find();
    for(int i = 0; i < res.size(); i++)
    {
        cout << res[i] <<" ";
    }
    return 0;
}