算法上机（三） 动态规划解决矩阵链问题

问题描述

给定矩阵序列A1, A2, …, An ，用普通的矩阵乘法，只加括号改变乘法顺序的情况下要求它们的积需要的最少乘法次数。
如A1：10100；
A2:1005;
A3:550
A1A2A3有两种情况：(A1A2)A3：101005次乘法
A1(A2A3）：1005*50次乘法。

动态规划解决

用p[0]p[1]代表第1个矩阵的阶数,p[1]p[2]代表第2个矩阵的阶数
m[i][j]表示AiAi+1…Aj的最少乘法次数.则有：
m[i][j]=min{m[i][k]+m[k+1][j]+Pi-1Pk*Pj
再自底向下的求出m[1][n]
代码如下：

#include <stdio.h>
#include <stdlib.h>

void  matrixChain(int n,int *p,int **m,int **s){
	for(int i=0;i<n;i++){
		m[i][i] = 0;
	}
	for(int h=1;h<n;h++){	//h表示i与j之差
		 for(int i=0;i<n-h;i++){	//除去n-1，i最大为n-2 
		 	int j=i+h;
		 	m[i][j]=999999999;	//初始值设为最大 
			for(int k=i;k<j;k++){
				int temp = m[i][k]+m[k+1][j]+p[i]*p[k+1]*p[j+1];
				if(temp<m[i][j]){
					m[i][j] = temp;
					s[i][j] = k;
				}
			}	 	
		} 
	}
}


int main (void){
	int n = 6;
	//p[0]*p[1]代表第1个矩阵的阶数,p[1]*p[2]代表第2个矩阵的阶数 
	int p[n+1] = {30,35,15,5,10,20,25};
	//m[i][j]表示第i+1个矩阵乘到第j+1个矩阵所需要的最少运算次数 
	int **m = (int **)malloc(n*sizeof(int *));
	for(int i=0;i<n;i++){
		m[i]=(int *)malloc(n*sizeof(int));
	}
	//s[i][j]表示第i+1个矩阵乘到第j+1个矩阵应在k=s[i][j]处分开 
	int **s = (int **)malloc(n*sizeof(int *));
	for(int i=0;i<n;i++){
		s[i]=(int *)malloc(n*sizeof(int));
	}
	matrixChain(n,p,m,s);
	printf("%d %d",m[1][3],s[1][3]);
}