1. 程式人生 > >矩陣鏈乘(Matrix Chain Multiplication)

矩陣鏈乘(Matrix Chain Multiplication)

矩陣鏈 alpha names namespace ror cati 次數 [0 expr

輸入n個矩陣的維度和一些矩陣鏈乘表達式,輸出乘法的次數。如果乘法無法進行,則輸出error。假定A是m*n矩陣,B是n*p矩陣,那麽A*B是m*p矩陣,乘法次數為m*n*p。如果A的列數不等於B的行數,則乘法無法進行。

例如,A是50*10的,B是10*20的,C是20*5的,則(A(BC))的乘法次數為10*20*5(BC的乘法次數) +50*10*5(A(BC)的乘法次數) = 3500。

#include<cstdio>
#include<stack>
#include<string>
#include<iostream>
using namespace std;


struct Matrix {
	int a, b;
	Matrix(int a = 0, int b = 0) :a(a), b(b) {}

}m[26];

stack<Matrix> s;
int main() {
	int n;
	cin >> n;
	for (int i = 0; i < n; i++) {
		string name;
		cin >> name;
		int k = name[0] - ‘A‘;
		cin >> m[k].a >> m[k].b;
	}
	string expr;
	while (cin >> expr) {
		int len = expr.length();
		bool error = false;
		int ans = 0;
		for (int i = 0; i < len; i++) {
			if (isalpha(expr[i])) s.push(m[expr[i] - ‘A‘]);
			else if (expr[i] == ‘)‘) {
				Matrix m2 = s.top(); s.pop();
				Matrix m1 = s.top(); s.pop();
				if (m1.b != m2.a) { error = true; break; }
				ans += m1.a *m1.b *m2.b;
				s.push(Matrix(m1.a, m2.b));
			}
		}
		if (error) printf("error\n"); else printf("%d\n", ans);
	}
	return 0;
}

  

矩陣鏈乘(Matrix Chain Multiplication)