【NOJ1048】计算矩阵连乘积【区间dp】

计算矩阵连乘积

时限:1000ms 内存限制:10000K  总时限:3000ms

描述

在科学计算中经常要计算矩阵的乘积。矩阵A和B可乘的条件是矩阵A的列数等于矩阵B的行数。若A是一个p×q的矩阵,B是一个q×r的矩阵,则其乘积C=AB是一个p×r的矩阵。计算C=AB总共需要p×q×r次乘法。
现在的问题是,给定n个矩阵{A1,A2,…,An}。其中Ai与Ai+1是可乘的,i=1,2,…,n-1。
要求计算出这n个矩阵的连乘积A1A2…An最少需要多少次乘法。

输入

输入数据的第一行是一个整树n(0 < n <= 10),表示矩阵的个数。
接下来的n行每行两个整数p,q( 0 < p,q < 100),分别表示一个矩阵的行数和列数。

输出

输出一个整数:计算连乘积最少需要乘法的次数。

输入样例

10
1 2
2 3
3 4
4 5
5 6
6 7
7 8
8 9
9 10
10 11

输出样例

438


递推关系:
      设计算A[i:j],1≤i≤j≤n,所需要的最少数乘次数m[i,j],则原问题的最优值为m[1,n]。

  •       当i=j时,A[i:j]=Ai,因此,m[i][i]=0,i=1,2,…,n
  •       当i<j时,若A[i:j]的最优次序在Ak和Ak+1之间断开,i<=k<j,则:m[i][j]=m[i][k]+m[k+1][j]+pi-1pkpj。由于在计算是并不知道断开点k的位置,所以k还未定。不过k的位置只有j-i个可能。因此,k是这j-i个位置使计算量达到最小的那个位置。

      综上,有递推关系如下:
  

#include <iostream>
#include <memory.h>
#include <algorithm>

using namespace std;

int main()
{
    int n;
    int p[11];
    int i,j,k;
    int temp;

    cin>>n;
    for(i=0;i<n;i++)
    {
        cin>>p[i]>>temp;
    }
    p[n]=temp;

    int dp[11][11];
    memset(dp,0,sizeof(dp));
    for(int d=1;d<=n-1;d++)  //区间长度
    {
        for(i=0;i<n-d;i++)  //起始位置
        {
            j=i+d;
            dp[i][j]=dp[i+1][j]+p[i]*p[i+1]*p[j+1];
            for(k=i+1;k<j;k++)//在k处断开
            {
                dp[i][j]=min(dp[i][j],dp[i][k]+dp[k+1][j]+p[i]*p[k+1]*p[j+1]);
            }
        }
    }

    cout<<dp[0][n-1]<<endl;
    return 0;
}

 

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值