uva 348

本文介绍了一个解决矩阵链乘法问题的动态规划算法,并详细展示了如何找到最优乘法顺序及其对应的路径。通过递归分解的方式,该算法不仅能够计算出最小的运算次数,还能明确展示出矩阵相乘的具体顺序。

dp还是比较容易  关键是输出路径 ....

#include <cstdlib>
#include <cstdio>
#include <cstring>
#define inf 0x7fffffff
using namespace std;
int n;
struct cc
{
    int a;
    int b;
    void f(int x, int y)
    {
        a = x;
        b = y;
    }
};
cc m[105];
int dp[105][105];
int path[105][105];

int DP(int l, int r)
{
    if(dp[l][r] != 0)
        return dp[l][r];
    if(l+1 == r)
    {
        dp[l][r] = m[l].a*m[l].b*m[r].b;
        return dp[l][r];
    }
    if(l == r)
        return 0;
    int mins = inf;
    for(int i = l; i < r; i++)
    {
        int t = DP(l, i) + DP(i+1, r) + m[l].a*m[i].b*m[r].b;
        if(t < mins)
        {
            mins = t;
            path[l][r] = i;
        }
    }
    dp[l][r] = mins;
    return mins;
}
void print(int l, int r)
{
    if(l == r)
    {
        printf("A%d",r);
        return;
    }
    if(l+1 == r)
    {
        printf("(A%d x A%d)",l,r);
        return;
    }
    printf("(");
    print(l, path[l][r]);
    printf(" x ");
    print(path[l][r]+1, r);
    printf(")");
}
int main()
{
    int ca = 1;
    while(scanf("%d",&n) == 1 && n)
    {
        for(int i = 1; i <= n; i++)
            {
                int x, y;
                scanf("%d %d",&x,&y);
                m[i].f(x, y);
            }
        memset(dp, 0, sizeof(dp));
        memset(path, 0, sizeof(path));
        DP(1,n);
        printf("Case %d: ", ca++);
        print(1,n);
        puts("");
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值