PTA——矩阵链相乘问题

矩阵链相乘问题

矩阵的乘法定义如下：设A是m×p的矩阵，B是p×n的矩阵，则A与B的乘积为m×n的矩阵，记作C=AB，其中，矩阵C中的第i行第j列元素c​ij
​​ 可以表示为：

当多个矩阵相乘时，采用不同的计算顺序所需的乘法次数不相同。例如，A是50×10的矩阵，B是10×20的矩阵，C是20×5的矩阵， 计算ABC有两种方式：(AB)C和A(BC)，前一种需要15000次乘法计算，后一种则只需3500次。
设A1 ,A2,…,An 为矩阵序列，Ai是阶为P​i−1 ∗Pi的矩阵(1≤i≤n)。试确定矩阵的乘法顺序，使得计算A1A2…A​n过程中元素相乘的总次数最少。
输入格式:
每个输入文件为一个测试用例，每个测试用例的第一行给出一个正整数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>
using namespace std;
 
const int size=1500;
int p[size];
int m[size][size];
int n;
 
void matrixchain(){
	int i,r,j,k;

	for(r=2;r<=n;r++){
		for(i=1;i<=n-r+1;i++){
			j=i+r-1;
			m[i][j]=m[i+1][j]+p[i-1]*p[i]*p[j];
			for(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;
				}
			}
		}
	}
}

 
int main(){
	cin>>n;

	for(int i=0;i<=n;i++)
		cin>>p[i];

	matrixchain(); 

	cout<<m[1][n]<<endl;
	
	return 0;
}

当size设为101时，段错误，需要把值设大，200也不行，不知道为什么。