矩阵的乘法定义如下:设A是m×p的矩阵,B是p×n的矩阵,则A与B的乘积为m×n的矩阵,记作C=AB,其中,矩阵C中的第i行第j列元素cij可以表示为:
当多个矩阵相乘时,采用不同的计算顺序所需的乘法次数不相同。例如,A是50×10的矩阵,B是10×20的矩阵,C是20×5的矩阵,
计算ABC有两种方式:(AB)C和A(BC),前一种需要15000次乘法计算,后一种则只需3500次。
设A1,A2,…,An为矩阵序列,Ai是阶为Pi−1∗Pi的矩阵(1≤i≤n)。试确定矩阵的乘法顺序,使得计算A1A2…An过程中元素相乘的总次数最少。
输入格式
每个输入文件为一个测试用例,每个测试用例的第一行给出一个正整数n(1≤n≤100),表示一共有n个矩阵A1,A2…An第二行给出n+1个整数P0,P1…Pn. 以空格分隔,其中1≤Pi≤100(0≤i≤n),第i个矩阵Ai是阶为Pi-1 * Pi 的矩阵。
输出格式
获得上述矩阵的乘积,所需的最少乘法次数。
输入样例
在这里给出一组输入。例如:
5
30 35 15 5 10 20
输出样例
在这里给出相应的输出。例如:
11875
#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
int main()
{
int n;
cin >> n;
int Form[n + 1]{};
int solve[n][n];
for (int i = 0; i <= n; i++)
{
if (i != n)
solve[i][i] = 0; //自己不划分
cin >> Form[i]; //输入矩阵数据
}
for (int l = 2; l <= n; l++)
{
for (int i = 0; i <= n - l; i++)
{
int j = i + l - 1; //算上i,取l个元素还需l-1就行
solve[i][j] = -1;
for (int k = i; k < j; k++)
/*
for (int k = i; k <= j; k++)
不能取k=j,当k=j时,下面k+1会溢出从而出错
*/
{
/*
int temp = solve[i][k] + solve[k][j] + Form[i-1] * Form[k] * Form[j];
发生冲突,一开始打算solve表建立为0~n-1,但是当i = 1时,Form发生越界访问。
int temp = solve[i][k] + solve[k][j] + Form[i] * Form[k + 1] * Form[j + 1];
修改后仍然报错,应该是i~k k+1~j
*/
int temp = solve[i][k] + solve[k + 1][j] + Form[i] * Form[k + 1] * Form[j + 1];
if (solve[i][j] == -1 || temp < solve[i][j])
{
solve[i][j] = temp;
}
}
}
}
cout << solve[0][n - 1] << endl;
return {};
}