矩阵链乘问题——求最小乘数

题目:

6个矩阵相乘,求需要乘法的最少数目。

A30X35  A35X15  A15X5  A5X10  A10X20 A20X25 (分别对应A1  A2  A3  A4  A5  A6, 后面是下标)

参考代码:

#include<stdio.h>
const int P[7]= {30,35,15,5,10,20,25};
const int length=6;

int m[7][7]={0}; //  m[i][j] 代码第i个数乘到第j个数(第一个数为i=1)
int s[7][7]={0};


void Pri(int s[][7],int i,int j)
{
 if(i==j)
  printf("A%d",i);
 else
 {
  printf("(");
  Pri(s,i,s[i][j]);
  printf(") (");
  Pri(s,s[i][j]+1,j);
  printf(") ");
 }
}

int main()
{
 int i,j,k;
 int min;
 i=j=k=0;
 i=1;
 for(i=1;i<=6;i++)  //遍历共有多少个数
 {
  for(j=i;j>=1;j--)
  {
   if(i==j)
   {
    m[j][i]=0;
   }
   else
   {
    int temp;
    min = 10000000;
    for(k=j;k<i;k++)
    {
     temp= m[j][k]+m[k+1][i]+ P[j-1]*P[k]*P[i];
     if(min>temp)
     {
      min=temp;
      s[j][i] = k;
     }
    }
    m[j][i]=min;
   }
  }

 }
 printf("%d\n",m[1][6]);
 Pri(s,1,6);
 return 0;
}

 

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
矩阵链乘问题是一个经典的动态规划问题,其目标是找到一种最优的方式来计算给定的一组矩阵的连乘积。这个问题可以通过动态规划算法来解决。 动态规划算法的基本思想是将问题分解成更小的子问题,并使用已知的信息来计算更大的问题。在矩阵链乘问题中,我们可以将问题分解成计算两个矩阵的乘积的子问题,并使用已知的信息来计算更大的问题。 具体来说,我们可以定义一个二维数组m,其中m[i][j]表示从第i个矩阵到第j个矩阵的最小计算代价。我们还可以定义一个二维数组s,其中s[i][j]表示从第i个矩阵到第j个矩阵的最优计算次序。 接下来,我们可以使用以下递归公式来计算m和s: m[i][j] = 0 (i = j) m[i][j] = min{m[i][k] + m[k+1][j] + ri*ck*cm} (i <= k < j) 其中,ri和ci分别表示第i个矩阵的行数和列数,cm表示两个矩阵相乘的计算代价。 使用上述递归公式,我们可以计算出所有的m[i][j]和s[i][j]。最终,我们可以通过s数组来构造出最优的计算次序,并使用m数组来计算最小的计算代价。 下面是一个Python实现的例子: ```python def matrix_chain_order(p): n = len(p) - 1 m = [[0] * n for i in range(n)] s = [[0] * n for i in range(n)] for l in range(2, n+1): for i in range(n-l+1): j = i + l - 1 m[i][j] = float('inf') for k in range(i, j): q = m[i][k] + m[k+1][j] + p[i]*p[k+1]*p[j+1] if q < m[i][j]: m[i][j] = q s[i][j] = k return m, s def print_optimal_parens(s, i, j): if i == j: print("A{}".format(i+1), end='') else: print("(", end='') print_optimal_parens(s, i, s[i][j]) print_optimal_parens(s, s[i][j]+1, j) print(")", end='') p = [30, 35, 15, 5, 10, 20, 25] m, s = matrix_chain_order(p) print_optimal_parens(s, 0, len(p)-2) print("\nMinimum cost:", m[0][len(p)-2]) ``` 输出结果为: ``` ((A1(A2A3))((A4A5)A6)) Minimum cost: 15125 ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值