动态规划——矩阵连乘

#include <iostream>using namespace std;void MatrixMul(int *p,int n,int**m ,int **s)///p记录矩阵维数，m最少乘法次数，s分法{    for(int i=0;i<n;i++)        m[i][i]=0;    for(int r=1;r<n;r++)        for(int i=0;i<n-r;i++)        {          int j=i+r;          m[i][j]=m[i+1][j]+p[i]*p[i+1]*p[j+1];          s[i][j]=i;          for(int k=i+1;k<j;k++)          {              int t=m[i][k]+m[k+1][j]+p[i]*p[k]*p[j+1];              if(t<m[i][j])              {                  m[i][j]=t;                  s[i][j]=k;              }          }        }}int main(){    int *p,**m,**s,n;    cin>>n;    m=new int*[n];    s=new int*[n];    for(int i=0;i<n;i++)    {        m[i]=new int[n];        s[i]=new int[n];    }    p=new int[n+1];    for(int i=0;i<n+1;i++)        cin>>p[i];    MatrixMul(p,n,m ,s);    cout<<m[0][n-1]<<endl;    return 0;}