动态规划之二

矩阵连乘 动态规划

　　题目描述：给定n个矩阵｛A1,A2,…,An｝，其中Ai与Ai+1是可乘的，i=1,2 ,…,n-1。如何确定计算矩阵连乘积的计算次序，使得依此次序计算矩阵连乘积需要的数乘次数最少。例如：

　　A1={30x35} ; A2={35x15} ;A3={15x5} ;A4={5x10} ;A5={10x20} ;A6={20x25} ;

最后的结果为：((A1(A2A3))((A4A5)A6))  最小的乘次为15125。

分析：先考虑一个问题，AiAi+1.....Aj，如果存在一个k使得AiAi+1...Ak和Ak+1.....Aj，这两部分的连乘使得AiAi+1....Aj的连乘次数最少，那么AiAi+1.....Ak这个子部分连乘次数最少的方案就是如此不加括号，Ak+1...Aj同样如此。现在就是考虑如何设置状态方程。用一个二维数组m[i][j]表示AiAi+1....Aj的最少次数，那么转移方程可得dp[i][j]=dp[i][k]+dp[k+1][j]+row[i]*col[k]*col[j]。row为行，col为列。当然了，对于i==j的情况，因为不存在连乘，故dp[i][j]==0(i==j)。

代码：

#include<iostream>
#include<algorithm>
#include<iomanip>
#define INF 0x7fffffff
using namespace std;
int row[105],col[105],dp[105][105];
int main()
{
int n,i=1,t;
cin>>n;
t=n;
while(t)
{
cin>>row[i]>>col[i];
i++;
t--;
}
for(int i=1;i<=n;i++)
{
dp[i][i]=0;
}
for(int l=2;l<=n;l++)//乘法链条长度 
{
for(int k=1;k<=n-l+1;k++)//这里的k之所以取n-l+1是因为s最大取值到倒数第二个矩阵，否则就会因为是单个矩阵而无法连乘。
{
int j=k+l-1;
dp[k][j]=INF;
for(int u=k;u<=j;u++)
{
dp[k][j]=min(dp[k][j],dp[k][u]+dp[u+1][j]+row[k]*col[u]*col[j]);
}
}
}
for(int z=1;z<=n;z++)
{
for(int v=1;v<=n;v++)
{
cout<<setw(5)<<left<<dp[z][v]<<" ";
}
cout<<endl;
}
cout<<dp[1][n];
}