最优矩阵链乘

附链接：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;
}