您的位置:首页 > 其它

FZU 1061 矩阵连乘

2016-05-04 16:27 344 查看

Problem 1061 矩阵连乘

Accept: 445 Submit: 1699

Time Limit: 1000 mSec Memory Limit : 32768 KB



Problem Description

给定n个矩阵{A1,A2,...,An},考察这n个矩阵的连乘积A1A2...An。由于矩阵乘法满足结合律,故计算矩阵的连乘积可以有许多不同的计算次序,这种计算次序可以用加括号的方式来确定。

矩阵连乘积的计算次序与其计算量有密切关系。例如,考察计算3个矩阵{A1,A2,A3}连乘积的例子。设这3个矩阵的维数分别为10*100,100*5,和5*50。若按(A1A2)A3计算,3个矩阵连乘积需要的数乘次数为10*100*5+10*5*50 = 7500。若按A1(A2A3)计算,则总共需要100*5*50+10*100*50 = 75000次数乘。

现在你的任务是对于一个确定的矩阵连乘方案,计算其需要的数乘次数。



Input

输入数据由多组数据组成。每组数据格式如下:

第一行是一个整数n (1≤n≤26),表示矩阵的个数。

接下来n行,每行有一个大写字母,表示矩阵的名字,后面有两个整数a,b,分别表示该矩阵的行数和列数,其中1<a,b<100。

第n+1行是一个矩阵连乘的表达式,由括号与大写字母组成,没有乘号与多余的空格。如果表达式中没有括号则按照从左到右的顺序计算,输入的括号保证能够配对。



Output

对于每组数据,输出仅一行包含一个整数,即将该矩阵连乘方案需要的数乘次数。如果运算过程中出现不满足矩阵乘法法则的情况(即左矩阵列数与右矩阵的行数不同),则输出“error”。



Sample Input

3A 10 100B 100 5C 5 50A(BC)



Sample Output

75000

栈应用

#include <iostream>
#include <string.h>
#include <stdlib.h>
#include <algorithm>
#include <math.h>
#include <stdio.h>
#include <map>

using namespace std;
struct Node
{
int a;
int b;
}e[105];
map<char,int> m;

char c;
char d[2005];
Node s[2005];
Node s2[2005];
int top2;
int top;
long long int sum;
bool ans;
int main()
{
int n;
while(scanf("%d",&n)!=EOF)
{
sum=0;
ans=true;
getchar();
for(int i=1;i<=n;i++)
{
scanf("%c ",&c);
m[c]=i;
scanf("%d%d",&e[i].a,&e[i].b);
getchar();
}
scanf("%s",d);
int len=strlen(d);
top=-1;top2=-1;
for(int i=0;i<len;i++)
{
if(d[i]=='(')
{
Node term;
term.a=0;term.b=0;
s[++top]=term;
}
else if(d[i]==')')
{
top2=-1;
while(s[top].a!=0&&top!=-1)
{
s2[++top2]=s[top];
top--;
}

if(top2==0)
s[top]=s2[top2];
else
{
if(s2[top2].b!=s2[top2-1].a)
ans=false;
sum+=s2[top2].a*s2[top2].b*s2[top2-1].b;
Node temp;temp.a=s2[top2].a;temp.b=s2[top2-1].b;
top2-=2;
while(top2!=-1)
{
if(temp.b!=s2[top2].a)
ans=false;
sum+=temp.a*temp.b*s2[top2].b;
temp.a=temp.a;temp.b=s2[top2].b;
top2--;
}
top--;
s[++top]=temp;
}
}
else
s[++top]=e[m[d[i]]];
}
top2=-1;
while(top!=-1)
{
s2[++top2]=s[top];
top--;
}
if(top2==0)
{
if(!ans)
printf("error\n");
else
printf("%lld\n",sum);
}
else
{
sum+=s2[top2].a*s2[top2].b*s2[top2-1].b;
if(s2[top2].b!=s2[top2-1].a)
ans=false;
Node temp;temp.a=s2[top2].a;temp.b=s2[top2-1].b;
top2-=2;
while(top2!=-1)
{
if(temp.b!=s2[top2].a)
ans=false;
sum+=temp.a*temp.b*s2[top2].b;
temp.a=temp.a;temp.b=s2[top2].b;
top2--;
}
if(!ans)
printf("error\n");
else
printf("%lld\n",sum);
}
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: