Description
给定n个矩阵{ A1,A2,…,An },保证Ai与Ai+1是可乘的,i = 1,2,…,n-1。考察这n个矩阵的连乘积A1A2…An。由于矩阵乘法满足结合律,故计算矩阵的连乘积可以有许多不同的计算次序。这种计算次序可以用加括号的方式来确定。现要求设计一个高效的算法,对给定的n个矩阵确定一个计算次序使得总的乘法次数最少,并输出该最优值。
Input
输入的第一行是单独一个整数T,表示案例的数目。每个案例的第一行是单独一个n ( 1 ≤ n ≤ 600 ),表示矩阵的个数。接下来第n行,依序分别对应第i个矩阵,每行包括两个整数xi,yi (1 ≤ i ≤ n,1 ≤ xi , yi ≤ 100 ),表示该矩阵的行数和列数。保证n个矩阵依序是可乘的。
Output
每个案例输出一个整数,表示最少需要的乘法次数。 运算过程及结果不会超出int范围。
Sample Input
1
4
50 10
10 40
40 30
30 5
Sample Output
10500
代码如下:
#include<stdio.h>
void PM(int p[],int N);
int main(){
int T,N,p[610];
scanf("%d",&T);
while(T--){
scanf("%d",&N);
for(int i=1;i<=N;i++){ //输入行列值
int a,b;
scanf("%d%d",&a,&b);
if(i==1){
p[i-1]=a;
p[i]=b;
}else{
p[i]=b;
}
}
PM(p,N);
}
return 0;
}
void PM(int p[],int N){
int m[600][600];
for(int i=1;i<=N;i++){//单个矩阵和为0
m[i][i]=0;
}
for(int r=2;r<=N;r++){//矩阵的规模(1-N)
for(int i=1;i<=N-r+1;i++){//规模r首矩阵的位置(i~N-r+1)
int j=i+r-1;//规模r末矩阵的位置
m[i][j]=m[i][i]+m[i+1][j]+p[i-1]*p[i]*p[j];//从第一个矩阵开始划分m(1)+m(2,5),设为最小,再和m(1,2)+m(3,5)....等比较,的出规模r的最小值
for(int k=i+1;k<j;k++){//判断有没有更小的
int mm=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
if(mm<m[i][j]){
m[i][j]=mm;
}
}
}
}
printf("%d\n",m[1][N]);
}