矩阵链（DP思想）

引入

按顺序排列的的三个矩阵 M1,M2,M3

计算三个矩阵相乘结果，有两种乘法

（M1 M2) M3

M1(M2 M3)

但两种乘法计算次数不同

三个矩阵维度如下

4 * 5 5 * 6 6 * 7

第一种计算次数

4 * 5 * 6 + 4 * 6 * 7

第二种计算次数

5 * 6 * 7 + 4 * 5 * 7

计算次数会有不同

随着矩阵数量变大，计算的方式也会变得更多样，对于n个矩阵不同计算方式总数我们称之为卡特兰数。在此不多赘述。

寻找最优计算方式

暴力搜索

枚举所有矩阵计算方式按部就班的计算每一次记算次数，由于矩阵数量的增加，对应的卡特兰数（即计算方式种类也会急剧增加，不可取）

DP

引入一个二维矩阵辅助，用于存储中间变量

DP [i] [j] 表示从 i 到 j 个矩阵里最优的计算方式所对应的计算次数

而从 i 到 j 中，我们要将矩阵链分为两部分（前提为这两部分的计算方式同样为最优，这就需要先计算出这两部分再进行整个部分的最优方式求解）我们假设在k个位置来截取（k>=i && k <= j)

则有

c [i] [j] = max(c [i] [j] , c [i] [k-1] + c[k] [ j] + 矩阵i的行 * 矩阵k-1的列 * 矩阵j的列 )

因此在 i j 这两层循环 里还需要 一层循环来寻找 i 到 j 里分割为两部分最优的次数

我们需要初始化c [i] [i] = 0 ,然后寻找 c [i] [i+1] 接着 c [i] [i+2] 最后得到 c [1] [ n]

//
// Created by 29273 on 2022-03-23.
//
#include "bits/stdc++.h"

using namespace std;
#define inf 0x3f3f3f3f
const int maxN = 1e3 + 10;
int a[maxN];
int dp[maxN][maxN];
int N;

void solve() {
    for (int d = 1; d <= N - 1; ++d) {
        for (int i = 1; i <= N - d; ++i) {
            int j = i + d;
            dp[i][j] = inf;

            for (int k = i + 1; k <= j; ++k) {
                dp[i][j] = min(dp[i][j], dp[i][k - 1] + dp[k][j] + a[i] * a[k] * a[j + 1]);
            }
        }
    }
}

int main() {
    cin >> N;
    for (int i = 1; i <= N; ++i) {
        cin >> a[i];
    }
    solve();
    for (int i = 1; i < N; ++i) {
        for (int j = 1; j < N; ++j) {
            cout << dp[i][j] << " ";
        }
        cout << endl;
    }
    return 0;
}