矩阵连乘算法模板

什么是矩阵连乘问题？

*给定n个矩阵{A1，A2，A3，...,An}，求这n个矩阵的连乘积A1A2...An，乘法次数最小值和矩阵连乘最优计算次序。

*由于矩阵乘法满足结合律，故计算矩阵的连乘积可以有许多不同的计算次序。

*矩阵Ai的维度为 pi-1 X pi 。

/*
已知矩阵维数的一维数组，求矩阵连乘最优计算次序和最少次数的乘法。
*/
#include<cstdio>
#include<iostream>
using namespace std;
const int maxn =105;
int p[maxn];/*矩阵的维数*/
int m[maxn][maxn];/*最优值，最少乘法次数*/
int s[maxn][maxn];/*断开的位置*/
int n;/*p的长度*/

void matrixChain(){ /*计算最优值*/
    for(int i=1;i<=n;i++)
        m[i][i]=0;
    for(int r=2;r<=n;r++){
        for(int i=1;i<=n-r+1;i++){
            int j=r+i-1;
            m[i][j]=m[i+1][j]+p[i-1]*p[i]*p[j];
            s[i][j]=i;
            for(int k=i+1;k<j;k++){
                int t=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
                if(t<m[i][j]){
                    m[i][j]=t;
                    s[i][j]=k;
                }
            }
        }
    }
}

void trackBack(int i,int j){/*求最优计算次序,从1~n-1*/
    if(i==j) return;
    trackBack(i,s[i][j]);
    trackBack(s[i][j]+1,j);
    printf("A[%d:%d]*A[%d:%d]\n",i,s[i][j],s[i][j]+1,j);
}

int main()
{
  //  freopen("input.txt","r",stdin);
    while(~scanf("%d",&n)){
        for(int i=0;i<n;i++){
            scanf("%d",&p[i]);
        }
        matrixChain();
        trackBack(1,n-1);/*矩阵连乘最优计算次序*/
        printf("%d\n",m[1][n-1]);/*最少次数的乘法*/
    }
    return 0;
}

例题：Multiplication Puzzle