【算法实验四】【动态规划】【1048】计算矩阵连乘积

时限：1000ms 内存限制：10000K  总时限：3000ms

描述

在科学计算中经常要计算矩阵的乘积。矩阵A和B可乘的条件是矩阵A的列数等于矩阵B的行数。若A是一个p×q的矩阵，B是一个q×r的矩阵，则其乘积C=AB是一个p×r的矩阵。计算C=AB总共需要p×q×r次乘法。
现在的问题是，给定n个矩阵{A1,A2,…,An}。其中Ai与Ai+1是可乘的，i=1,2,…,n-1。
要求计算出这n个矩阵的连乘积A1A2…An最少需要多少次乘法。

 

输入

输入数据的第一行是一个整树n（0 < n <= 10），表示矩阵的个数。
接下来的n行每行两个整数p,q( 0 < p,q < 100)，分别表示一个矩阵的行数和列数。

 

输出

输出一个整数：计算连乘积最少需要乘法的次数。

 

输入样例

10
1 2
2 3
3 4
4 5
5 6
6 7
7 8
8 9
9 10
10 11

 

输出样例

438

#include<stdio.h>
#include<iostream>
using namespace std;

int matrixchain();
void readdata();

int p[101];//存放矩阵信息 
int m[101][101];//m[i][j] 矩阵Mi到Mj相乘最小计算次数
int s[101][101];//s[i][j] 矩阵Mi到Mj相乘 分开位置
int n;//矩阵数目
int r[100], c[100];

int main() {
	readdata();
	cout << matrixchain() << endl;
	return 0;
}

int matrixchain() {
	for (int i = 1; i < n; i++)
		m[i][i] = 0;
	for(int r=1;r<n;r++)//r为i与j的差
		for (int i = 1; i < n; i++) {
			int j = i + r;
			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];
}

void readdata() {
	cin >> n;
	for (int i = 0; i < n; i++) {
		cin >> r[i] >> c[i];
	}
	for (int i = 0; i < n; i++)
		p[i] = r[i];
	p[n] = c[n - 1];
	
}