矩陣鏈乘(Matrix Chain Multiplication)
阿新 • • 發佈:2017-06-24
矩陣鏈 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)