UVa348 最优矩阵链乘 （记忆化搜搜&&dp&&输出路径）

 最优矩阵链乘：

一个n*m的矩阵由n行m列共n*m排列而成。两个矩阵A和B可以相乘当且仅当A的列数等于B的行数。一个n*m的矩阵乘m*p的矩阵，运算量为n*m*p。

矩阵乘法不满足分配律，但满足结合律。因此A*B*C既可以按顺序（A*B）C也可以按A（B*C）来进行。假设A、B、C分别是2*3、3*4、4*5的，则（A*B）C运算量是2*3*4+2*4*5=64，A（B*C）的运算量是3*4*5*2*3*5=90.显然第一种顺序节省运算量。给出n个矩阵组成的序列，设计一种方法把他们依次乘起来，使得总的运算量尽量小。假设第i个矩阵A[i]是P[i-1]*P[i]的。

 

记忆化搜索

#include<iostream>
#include<cstring>
using namespace std;
const int inf=0x3f3f3f3f;
int dp[20][20],pre[20][20],x[20],y[20];
int dfs(int a,int b){
    if(dp[a][b]!=-1)return dp[a][b];
    pre[a][b]=a;
    if(a==b)return dp[a][b]=0;
    dp[a][b]=inf;
    int t;
    for(int i=a;i<b;i++){
        t=dfs(a,i)+dfs(i+1,b)+x[a]*y[i]*y[b];
        if(t<dp[a][b]){
            dp[a][b]=t;
            pre[a][b]=i;
        }
    }
    return dp[a][b];
}
void print(int a,int b){
    if(a>b)return;
    if(a==b)
        cout<<"A"<<a+1;
    else{
        cout<<"(";
        print(a,pre[a][b]);
        cout<<" x ";
        print(pre[a][b]+1,b);
        cout<<")";
    }
}
int main(){
    int n,flag=0;
    while(cin>>n&&n){
        memset(dp,-1,sizeof(dp));
        for(int i=0;i<n;i++){
            cin>>x[i]>>y[i];
        }
        dfs(0,n-1);
        cout<<"Case "<<++flag<<": ";
        print(0,n-1);
        cout<<endl;
    }
    return 0;
}