DP算法之矩阵链乘

最新推荐文章于 2020-10-22 16:26:20 发布

原创最新推荐文章于 2020-10-22 16:26:20 发布 · 1.3k 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#dp算法 #矩阵链称

开心的mdd

时间限制：1000 ms | 内存限制：65535 KB

难度：3

描述

himdd有一天闲着无聊，随手拿了一本书，随手翻到一页，上面描述了一个神奇的问题，貌似是一个和矩阵有关的东西。

给出三个矩阵和其行列A1（10*100），A2（100*5），A3（5*50）。现在himdd要算出计算矩阵所要的乘法次数，他发现不同的计算次序，所要的乘法次数也不一样，

如：

(A1*A2)*A3 : 10*100*5+5*10*50=7500;

A1*(A2*A3) : 5*100*50+10*100*50 =75000;

他想知道计算矩阵所要的最少乘法次数是多少，很快一个解法就诞生了，有点小happy~~现在他想问问你是否也能找出一个解法呢？

注意：矩阵不可改变顺序。

输入

有多组测试数据（<=100)，每组表述如下：
第一行，有一个整数n矩阵的个数（1<=n<=100)
接下来有n行
第i行有两整数，r,c表示第i个矩阵的行列;(1<=r,c<=100)

输出

输出计算矩阵所要的最少乘法次数。

样例输入

样例输出

AC代码

#include <stdio.h>
#include <string.h>
int min(int a,int b)
{
    return a==0 ? b : (a < b ? a : b);
}
int main()
{
    int n,p[102],dp[102][102];
    while(~scanf("%d",&n))
    {
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++)                       //第i的矩阵的行号和列号分别是p[i-1],p[i]
            scanf("%d%d",&p[i-1],&p[i]);
        for(int d=1;d<n;d++)                        //控制计算链长  实际链长为 d+1
            for(int i=1;i<=n-d;i++)                 //链的起始值    同是也是二维数字的行号
            { 
                int j = i+d;                        //链的结尾值    同时也是二维数组的列号
                for(int k=i;k<j;k++)                //根据状态转移式求出dp[i][j]
                    dp[i][j] = min(dp[i][j],dp[i][k]+dp[k+1][j]+p[i-1]*p[k]*p[j]);
            }
        printf("%d\n",dp[1][n]);                    //dp[1][n]表示从1到n的最少乘法次数
    }
	return 0;
}

dp值并为全满

如dp[0][0]并没有使用，从下标1开始只为整体方便。