矩阵连乘最小计算次数 C语言

先上代码

#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <cmath>

const int N = 330;
struct pii
{
	int first, second;
};
pii matrix[N];
int f[N][N];

void init()
{
	memset(f, 0x3f, sizeof f);
	for (int i = 0; i < N; i++) f[i][i] = 0;
}

int min(int x, int y)
{
	if (x > y) return y;
	return x;
}

int main()
{
	init();
	int n;
	scanf("%d", &n);
	for (int i = 1; i <= n; i++)
	{
		scanf("%d%d", &matrix[i].first, &matrix[i].second);
	}
	for (int j = 1; j < n; j++)
	{
		for (int i = 1; i+j <= n; i++)
		{
			for (int k = i; k <= i + j; k++)
			{
				f[i][i+j] = min(f[i][i+j], f[i][k] + f[k + 1][i+j] + matrix[i].first * matrix[k].second * matrix[i+j].second);
			}
		}
	}
	printf("%d\n", f[1][n]);
	return 0;
}

思路

给出 i i i j j j两数 f [ i ] [ j ] f[i][j] f[i][j]表示区间 i i i i + j i+j i+j的连乘最小次数。用m的结构体(用pair也行),两相邻矩阵相乘的计算次数是 m [ i ] . f i r s t ∗ m [ i ] . s e c o n d ∗ m [ i + 1 ] . s e c o n d m[i].first * m[i].second * m[i+1].second m[i].firstm[i].secondm[i+1].second 或者 m [ i ] . f i r s t ∗ m [ i + 1 ] . f i r s t ∗ m [ i + 1 ] . s e c o n d m[i].first * m[i+1].first * m[i+1].second m[i].firstm[i+1].firstm[i+1].second 。推广开来就可以得出我们的动态规划f的层间关系。

层间关系

f [ i ] [ j ] = m i n ( f [ i ] [ j ] , f [ i ] [ k ] + f [ k + 1 ] [ j ] + m [ i ] . f i r s t ∗ m [ k ] . s e c o n d ∗ m [ j ] . s e c o n d ] f[i][j] = min(f[i][j],f[i][k]+f[k+1][j]+m[i].first * m[k].second * m[j].second] f[i][j]=min(f[i][j],f[i][k]+f[k+1][j]+m[i].firstm[k].secondm[j].second];
极限便是 f [ i ] [ i ] = 0 f[i][i] = 0 f[i][i]=0(自乘次数为 0 0 0)

for循环思路

因为这里增长的线性关系是 i i i i + j i+j i+j的搜索宽度,所以第一层for代表宽度,第二层代表起点,这两层就构造了所有的以i为起点的所有向后宽度的搜多范围,并且先看每一个起点,遍历完所有的起点后再增加宽度。

  • 2
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值