【SSL】2021-12-04 1596.矩阵链相乘

原题网址

由于某些原因，这个网址会进不去…

题目描述

假设我们要用标准的矩阵乘法来计算$M_1,M_2,M_3$三个矩阵的乘积$M_1{M}_2{M}_3$，这三个矩阵的维度分别是$2\times 10,10\times 2,2\times 10$，如果把$M_1,M_2$相乘，然后再与$M_3$相乘，那么要$2\times 10\times 2+2\times 2\times 10=80$，如果代之以用$M_2,M_3$相乘的结果去成$M_1$，那么矩阵乘法的次数变成了$10\times 2\times 10+2\times 10\times 10=400$，执行$M_1(M_2{M}_3)$耗费的时间是执行乘法$(M_1{M}_2)M_3$的$5$倍。

格式

输入格式

$n$表示矩阵的个数$(\le 100)$
$n+1$个数,表示矩阵$(\le 100)$

输出格式

最小的乘法次数

样例

输入样例

5
5 10 4 6 10 2

输出样例

348

解题思路

一般来说，顺序数等于这个$n$个矩阵每一种可能的途径放置括号数的方法数。设$f(n)$是求$n$个知阵乘积的所有放置括号的方法数，假定要进行下面的乘法：
$(A_1{A}_2...A_k)\times (A_{k+1}{A}_{k+2}...A_n)$
那么，对于前$k$个矩阵有$f(k)$对于中的每一种方法，可对余下的$f(n-k)$个矩阵放置括号的方法，总共有$f(k)\times f(n-k)$种方法。由于可以假设$k$是$1$到$n-1$中任意值，对于$n$个矩阵放置括号弧的所有方法数由
$f(n)={\sum }_{k=1}^{n-1}f(n)f(n-k)$
给出。显然，两个矩阵相乘只有一种方法，三个矩阵相乘有两种方法。因此$f(2)=1,f(3)=2$，为了使递推有意义，令$f(1)=1$。我们可以证明
$f(n)=C_{2n-2}^{n-1}/n$
这个递推式就产生的CATALAN数（即卡特兰数）
百度链接
$1,1,2,5,14,42,132,429,1430,4862,16796...$
因此，$10$个矩阵相乘存在$4862$种方法。
矩阵$A$和$B$是可乘的条件是矩阵$A$的列数等于矩阵$B$的行数。若$A$是一个$p\times q$的矩阵，$B$是一个$q\times r$的矩阵，则其乘积$C=AB$是一个$p\times r$矩阵。在计算$C$的标准算法中，主要计算量在三重循环，总共需要$p\times q\times r$次数成。
$3\times 22\times 23\times 2$

$122181$
$21\ast 30=73$
$0130$

将矩阵连乘积$A_i{A}_{i+1}\dots A_j$简记为$A[i:j]$。设最优次序在矩阵$A_k$和$A_{k+1}$A之间将矩阵链断开，$1\le k<n$，则完全加括号的方式为$((A_1\dots A_{k-1})(A_k\dots A_n))$。计算$F[1][n]$的一个最优次序所包含的计算矩阵子链$F[1][k-1]$和$F[k][n]$的次序也是最优的。因此，矩阵连乘积计算次序问题的最优解包含着其子问题的最优解。
f [ i ] [ j ] = m i n { f [ i ] [ k − 1 ] + f [ k ] [ j ] + r i r k r j + 1 } ( i < k ≤ j ) f[i][j]=min\{f[i][k-1]+f[k][j]+r_ir_kr_{j+1}\}(i<k\le j) f[i][j]=min{f[i][k−1]+f[k][j]+ri​rk​rj+1​}(i<k≤j)
而本问题的问题求解则为：
f [ 1 ] [ n ] = m i n { f [ 1 ] [ k − 1 ] + f [ k ] [ n ] + r 1 r k r n + 1 } f[1][n]=min\{f[1][k-1]+f[k][n]+r_1r_kr_{n+1}\} f[1][n]=min{f[1][k−1]+f[k][n]+r1​rk​rn+1​}

Code

#include<iostream>
#include<istream>
#include<ostream>
#include<sstream>
#include<fstream>
#include<iomanip>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<string>
#include<cstring>
#include<map>
#include<vector>
using namespace std;
const int maxn=100,inf=2147483647;
int n,m[maxn*2],f[maxn*2][maxn*2];
void init() {
    ios::sync_with_stdio(false); // 快读
    cin.tie(); // 快读
    cin>>n;
    for(int i=1;i<=n+1;i++) cin>>m[i];
}
int main() {
    init();
    for(int i=1;i<=n;i++) f[i][i]=0;
    for(int i=1;i<n;i++) {
        for(int j=1;j<=n-i;j++) {
            int k=i+j;
            f[j][k]=inf;
            for(int l=j+1;l<=k;l++) f[j][k]=min(f[j][k],f[j][l-1]+f[l][k]+m[j]*m[k+1]*m[l]);
        }
    }
    cout<<f[1][n];
    return 0;
}