附链接:http://blog.csdn.net/zhaochengyuan/article/details/7863542
题目大意:
一个n*m矩阵由n行m列共n*m个数排列而成。两个矩阵A和B可以相乘当且仅当A的列数等于B的行数。一个N*M的矩阵乘以一个M*P的矩阵等于一个N*P的矩阵,运算量为nmp。
矩阵乘法满足结合律,A*B*C可以表示成(A*B)*C或者是A*(B*C),两者的运算量却不同。例如当A=2*3 B=3*4 C=4*5时,(A*B)*C=64而A*(B*C)=90。显然第一种顺序节省运算量。
现在给出N个矩阵,并输入N+1个数,第i个矩阵是a[i-1]*a[i]。
基本思路:
然后就是谁和谁先结合的问题,然后每一小的结合部分又存在谁和谁先结合的问题,这是很明显的最优子结构 ,然后很明显用记忆化搜索更为简单易懂,然后就是附上代码以及注释;
代码如下:
#include<iostream>
#include<iomanip>
#include<algorithm>
#include<sstream>
#include<string>
#include<vector>
#include<stack>
#include<map>
#include<list>
#include<set>
#include<queue>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
using namespace std;
const int maxn = 100+10;
int p[maxn],dp[maxn][maxn],s[maxn][maxn];
int chain(int i,int j)
{
if(dp[i][j]) return dp[i][j];//套路;
if(i==j) return 0;//套路;
int v=chain(i,i)+chain(i+1,j)+p[i-1]*p[i]*p[j];
s[i][j]=i;
for(int k=i+1;k<j;k++)//选k为断点;
{
int t=chain(i,k)+chain(k+1,j)+p[i-1]*p[k]*p[j];
if(t<v)//更新最优解;
{
v=t;
s[i][j]=k;
}
}
dp[i][j]=v;//套路
return v;
}
int main()
{
memset(dp,0,sizeof(dp));
int n;
cin>>n;
for(int i=0;i<=n;i++) cin>>p[i];
cout<<chain(1,n)<<endl;
return 0;
}