动态规划之矩阵链乘法

问题

给定n个矩阵的链<A 1 ,A 2 ,…,A n >，其中i=1,2,…,n，矩阵A i的维数为p i-1 ×p i 。求一个完全“括号化方案”，使得计算乘积A1 A2 …A n 所需的标量乘法次数最小。

解析

　　对于矩阵链乘法问题，我们将所有对于1≤i≤j≤n确定Ai Ai+1 … j 的最小代价括号方案作为子问题。令m[i,j]表示计算矩阵A i,j 所需要的标量乘法的次数最小值，则最优解就是计算A i…n所需的最低代价就是m[1,n]。
① 对于i=j的情况下，显然有m=0，不需要做任何标量乘法运算。所以，对于所有的i=1、2…n，m[i,i] =0。
②当i < j的情况，就按照最优括号化方案的结构特征进行计算m[i,j]。假设最优括号化方案的分割点在矩阵Ak和Ak+1之间，那么m的值就是Ai…k和Ak+1…j的代价加上两者量程的代价的最小值。即。该公式的假设是最优分割点是已知的，但是实际上不知道。然而，k只有j-i中情况取值。由于最优分割点k必定在i~j内取得，只需要检查所有可能的情况，找到最优解即可。可以得出一个递归公式
　　m只是给出了子问题最优解的代价，但是并未给出构造最优解的足够信息(即分割点的位置信息)。所以，在此基础之上，我们使用一个二维数组s[i,j]来保存 Ai Ai+1 …Aj 的分割点位置k。

设计

分析

T(n)=O(n³)

源码

#include<bits/stdc++.h>
using namespace std;
const int maxn=1010;
const int inf=0x3f3f3f3f;

int dp[maxn][maxn],p[maxn],s[maxn][maxn];

void print_PARENS(int i,int j){
	if(i==j) printf("A%d",i);
	else{
		printf("(");
		print_PARENS(i,s[i][j]);
		print_PARENS(s[i][j]+1,j);
		printf(")");
	}
}

int main(){
	int n;
	scanf("%d",&n);
	for(int i=0;i<n;i++) scanf("%d",&p[i]);
	for(int i=1;i<=n;i++){
		dp[i][i]=0;
	}
	for(int len=2;len<=n;len++){
		for(int i=1;i+len-1<=n;i++){
			int j=i+len-1;
			dp[i][j]=inf;
			for(int k=i;k<j;k++){
				int tmp=dp[i][k]+dp[k+1][j]+p[i-1]*p[k]*p[j];
				if(tmp<dp[i][j]){
					dp[i][j]=tmp;
					s[i][j]=k;
				}
			}
		}
	}
	printf("最优解：%d\n",dp[1][n-1]);
	print_PARENS(1,n-1);
} 
/*
7
30 35 15 5 10 20 25
*/

Github:
https://github.com/Geedhayb/Geed/blob/master/MATRIX_CHAIN_ORDER.cpp