矩阵链乘

矩阵链乘问题，是个动态规划大展身手的地方。如果您还不了解矩阵链乘，建议您先到其他博主的博文上学习，这里目的是要给出矩阵链乘精简的C++代码。
参考：《算法导论》

测试数据

有若干个矩阵（行列给出）

矩阵链乘目的

寻找乘法次数最少的矩阵运算方案

原理

矩阵相乘的结合律

方法

动态规划

思路

对于A1,A2,A3…An的矩阵，肯定先用“一刀”将矩阵链分成两部分，这样，大问题就能化解为子问题，并且，该问题具有最有子结构。那这一刀实际上会有若干种可能，于是，我们需要对这每一种可能都做遍历，更新最小值。

变量准备

定义dp[i,j]为矩阵链Ai,Ai+1,…,Aj的最小相乘次数
定义s[i,j]为矩阵链Ai,Ai+1,…,Aj得到相乘最少次数时，“一刀”砍在的位置（为了方便重构解）
例如，dp[1,3]=6，s[1,3]=2，表示A1A2A3矩阵相乘的最少次数为6，为了使相乘次数为6，这一刀应该这样砍：(A1A2)A3
那么，这样一来，最终次数就应该是：dp[1,n]

输入

9
50 10
10 20
20 5
5 10
10 20
20 25
25 30
30 35
35 15

C++源代码

#include <iostream>
#include <algorithm>
#define INF 0x7fffffff
#define N 12
using namespace std;

int n; // 矩阵个数
struct matrix {
    int row; //行数
    int col; //列数
}a[N];
int dp[N][N]; // dp[i,j]为矩阵链Ai,Ai+1,...,Aj的最小相乘次数
int s[N][N]; // s[i,j]为矩阵链Ai,Ai+1,...,Aj得到相乘最少次数时，“一刀”砍在的位置（为了方便重构解）

void print(int i, int j) {
    if (i == j) {
        cout << "A" << i;
    } else {
        cout << "(";
        print(i,s[i][j]);
        print(s[i][j]+1,j);
        cout << ")";
    }
}

int main() {
    cin >> n;
    for (int i = 1; i <= n; ++i) {
        cin >> a[i].row >> a[i].col;
    }
    for (int l = 2; l <= n; ++l) { // 枚举矩阵链的长度，肯定要从短矩阵链开始枚举
        for (int i = 1; i+l-1 <= n; ++i) { // 枚举矩阵链的起点
            int j = i+l-1; // 枚举矩阵链的终点
            dp[i][j] = INF;
            for (int k = i; k <= j-1; ++k) { // 枚举此段所有可能切割点
                int tmp = dp[i][k]+dp[k+1][j]+a[i].row*a[k].col*a[j].col;
                if (tmp < dp[i][j]) {
                    dp[i][j] = tmp;
                    s[i][j] = k;
                }
            }
        }
    }
    cout << dp[1][n] << endl;
    print(1,n);
    return 0;
}

如有错误，欢迎指正。