poj 1651 矩阵连乘

今天终于学会了矩阵连乘~~~矩阵连乘第一次AC~纪念一下，并且感谢这位博主的讲解http://blog.sina.com.cn/s/blog_64018c250100s123.html

如果不会矩阵连乘的同学可以去看看~~。DP经典之一

思路： 把这n个数写入 p[0]......p[n-1]这个数组中，然后  就相当于有n-1个矩阵连乘，

p[0]是第一个矩阵的 行 ， p[1]是第一个矩阵的  列  ......这样推下去，就可以出结果了

附上代码：

#include <iostream>
#include <stdio.h>
#define N 150
using namespace std;
int m[N][N],p[N],n;
int dp_multiply()
{
 n-=1;                                  //构建n-1个矩阵连乘
 for(int i=1;i<=n;i++)          //m[i][i]表示自己与自己相乘，那么最小的乘法次数就为0
   m[i][i]=0;
 for(int r=2;r<=n;r++)      //r表示矩阵的个数，  之所以从下面往上面推，是因为 有个状态转移方程 ;   m[i][j]=(m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j])i<=k<j
  {
    for(int i=1;i<=n-r+1;i++)//表示 r个矩阵个数已定的情况下，i的其实
     { 
       int j=i+r-1;                   
       m[i][j]=m[i+1][j]+p[i-1]*p[i]*p[j];      
       for(int k=i+1;k<j;k++)                              //dp
       {
         int tem=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
         m[i][j]=tem<m[i][j]?tem:m[i][j];
       }
     }   
  
  } 


 return m[1][n];
}
int main()
{
 while(cin>>n)
 {
  memset(m,0,sizeof(m));
  memset(p,0,sizeof(p));
  for(int i=0;i<n;i++)
    cin>>p[i];
  cout<<dp_multiply()<<endl;  
 }
 return 0;   
}