zoj 1276 Optimal Array Multiplication Sequence
2009-05-15 21:15
351 查看
ZOJ Problem Set - 1276
Optimal Array Multiplication Sequence
Time Limit: 1 Second Memory Limit: 32768 KB Special Judge
Given two arrays A and B, we can determine the array C = A B using the standard definition of matrix multiplication:
The number of columns in the A array must be the same as the number of rows in the B array. Notationally, let's say that rows(A) and columns(A) are the number of rows and columns, respectively, in the A array. The number of individual multiplications required to compute the entire C array (which will have the same number of rows as A and the same number of columns as B) is then rows(A) columns(B) columns(A). For example, if A is a 10 x 20 array, and B is a 20 x 15 array, it will take 10 x 15 x 20, or 3000 multiplications to compute the C array.
To perform multiplication of more than two arrays we have a choice of how to proceed. For example, if X, Y, and Z are arrays, then to compute X Y Z we could either compute (X Y) Z or X (Y Z). Suppose X is a 5 x 10 array, Y is a 10 x 20 array, and Z is a 20 x 35 array. Let's look at the number of multiplications required to compute the product using the two different sequences:
(X Y) Z
5 x 20 x 10 = 1000 multiplications to determine the product (X Y), a 5 x 20 array.
Then 5 x 35 x 20 = 3500 multiplications to determine the final result.
Total multiplications: 4500.
X (Y Z)
10 x 35 x 20 = 7000 multiplications to determine the product (Y Z), a 10 x 35 array.
Then 5 x 35 x 10 = 1750 multiplications to determine the final result.
Total multiplications: 8750.
Clearly we'll be able to compute (X Y) Z using fewer individual multiplications.
Given the size of each array in a sequence of arrays to be multiplied, you are to determine an optimal computational sequence. Optimality, for this problem, is relative to the number of individual multiplcations required.
Input
For each array in the multiple sequences of arrays to be multiplied you will be given only the dimensions of the array. Each sequence will consist of an integer N which indicates the number of arrays to be multiplied, and then N pairs of integers, each pair giving the number of rows and columns in an array; the order in which the dimensions are given is the same as the order in which the arrays are to be multiplied. A value of zero for N indicates the end of the input. N will be no larger than 10.
Output
Assume the arrays are named A1, A2, ..., AN. Your output for each input case is to be a line containing a parenthesized expression clearly indicating the order in which the arrays are to be multiplied. Prefix the output for each case with the case number (they are sequentially numbered, starting with 1). Your output should strongly resemble that shown in the samples shown below. If, by chance, there are multiple correct sequences, any of these will be accepted as a valid answer.
Sample Input
3
1 5
5 20
20 1
3
5 10
10 20
20 35
6
30 35
35 15
15 5
5 10
10 20
20 25
0
Sample Output
Case 1: (A1 x (A2 x A3))
Case 2: ((A1 x A2) x A3)
Case 3: ((A1 x (A2 x A3)) x ((A4 x A5) x A6))
Source: North Central North America 1996
Submit Status
//1858191 2009-05-07 11:07:29 Accepted 1276 C++ 0 184 Wpl
#include <iostream>
#define MAX 11
using namespace std;
int dp[MAX][MAX],path[MAX][MAX],data[MAX];
int n;
void Init()
{
int i=0,j=1,m=n;
while(m--)
{
scanf("%d%d",&data[i],&data[j]);
i++;
j++;
}
}
void MatrixChain()
{
int i,j,k,r;
for(i=1;i<=n;i++)
dp[i][i]=0;
for(r=2;r<=n;r++)
for(i=1;i<=n-r+1;i++)
{
j=i+r-1;
dp[i][j]=dp[i][i]+dp[i+1][j]+data[i-1]*data[i]*data[j];
path[i][j]=i; //默认在i处断开
for(k=i+1;k<j;k++)
{
int temp=dp[i][k]+dp[k+1][j]+data[i-1]*data[k]*data[j];
if(temp<dp[i][j])
{
dp[i][j]=temp;
path[i][j]=k;
}
}
}
}
void TraceBack(int i,int j)
{
if(i==j)
{
cout<<"A"<<i;
}
else if(i==j-1)
{
cout<<"A"<<i<<" x A"<<j;
}
else
{
if(i!=path[i][j])
cout << "(";
TraceBack(i,path[i][j]);
if (i!=path[i][j])
cout<<") x ";
else
cout<<" x ";
if (path[i][j]+1!=j)
cout<<"(";
TraceBack(path[i][j]+1,j);
if (path[i][j]+1!=j)
cout << ")";
}
}
int main()
{
int zz=1;
while(scanf("%d",&n)!=EOF&&n!=0)
{
Init();
MatrixChain();
printf("Case %d: (",zz++);
TraceBack(1,n);
printf(")\n");
}
return 0;
}
Optimal Array Multiplication Sequence
Time Limit: 1 Second Memory Limit: 32768 KB Special Judge
Given two arrays A and B, we can determine the array C = A B using the standard definition of matrix multiplication:
The number of columns in the A array must be the same as the number of rows in the B array. Notationally, let's say that rows(A) and columns(A) are the number of rows and columns, respectively, in the A array. The number of individual multiplications required to compute the entire C array (which will have the same number of rows as A and the same number of columns as B) is then rows(A) columns(B) columns(A). For example, if A is a 10 x 20 array, and B is a 20 x 15 array, it will take 10 x 15 x 20, or 3000 multiplications to compute the C array.
To perform multiplication of more than two arrays we have a choice of how to proceed. For example, if X, Y, and Z are arrays, then to compute X Y Z we could either compute (X Y) Z or X (Y Z). Suppose X is a 5 x 10 array, Y is a 10 x 20 array, and Z is a 20 x 35 array. Let's look at the number of multiplications required to compute the product using the two different sequences:
(X Y) Z
5 x 20 x 10 = 1000 multiplications to determine the product (X Y), a 5 x 20 array.
Then 5 x 35 x 20 = 3500 multiplications to determine the final result.
Total multiplications: 4500.
X (Y Z)
10 x 35 x 20 = 7000 multiplications to determine the product (Y Z), a 10 x 35 array.
Then 5 x 35 x 10 = 1750 multiplications to determine the final result.
Total multiplications: 8750.
Clearly we'll be able to compute (X Y) Z using fewer individual multiplications.
Given the size of each array in a sequence of arrays to be multiplied, you are to determine an optimal computational sequence. Optimality, for this problem, is relative to the number of individual multiplcations required.
Input
For each array in the multiple sequences of arrays to be multiplied you will be given only the dimensions of the array. Each sequence will consist of an integer N which indicates the number of arrays to be multiplied, and then N pairs of integers, each pair giving the number of rows and columns in an array; the order in which the dimensions are given is the same as the order in which the arrays are to be multiplied. A value of zero for N indicates the end of the input. N will be no larger than 10.
Output
Assume the arrays are named A1, A2, ..., AN. Your output for each input case is to be a line containing a parenthesized expression clearly indicating the order in which the arrays are to be multiplied. Prefix the output for each case with the case number (they are sequentially numbered, starting with 1). Your output should strongly resemble that shown in the samples shown below. If, by chance, there are multiple correct sequences, any of these will be accepted as a valid answer.
Sample Input
3
1 5
5 20
20 1
3
5 10
10 20
20 35
6
30 35
35 15
15 5
5 10
10 20
20 25
0
Sample Output
Case 1: (A1 x (A2 x A3))
Case 2: ((A1 x A2) x A3)
Case 3: ((A1 x (A2 x A3)) x ((A4 x A5) x A6))
Source: North Central North America 1996
Submit Status
//1858191 2009-05-07 11:07:29 Accepted 1276 C++ 0 184 Wpl
#include <iostream>
#define MAX 11
using namespace std;
int dp[MAX][MAX],path[MAX][MAX],data[MAX];
int n;
void Init()
{
int i=0,j=1,m=n;
while(m--)
{
scanf("%d%d",&data[i],&data[j]);
i++;
j++;
}
}
void MatrixChain()
{
int i,j,k,r;
for(i=1;i<=n;i++)
dp[i][i]=0;
for(r=2;r<=n;r++)
for(i=1;i<=n-r+1;i++)
{
j=i+r-1;
dp[i][j]=dp[i][i]+dp[i+1][j]+data[i-1]*data[i]*data[j];
path[i][j]=i; //默认在i处断开
for(k=i+1;k<j;k++)
{
int temp=dp[i][k]+dp[k+1][j]+data[i-1]*data[k]*data[j];
if(temp<dp[i][j])
{
dp[i][j]=temp;
path[i][j]=k;
}
}
}
}
void TraceBack(int i,int j)
{
if(i==j)
{
cout<<"A"<<i;
}
else if(i==j-1)
{
cout<<"A"<<i<<" x A"<<j;
}
else
{
if(i!=path[i][j])
cout << "(";
TraceBack(i,path[i][j]);
if (i!=path[i][j])
cout<<") x ";
else
cout<<" x ";
if (path[i][j]+1!=j)
cout<<"(";
TraceBack(path[i][j]+1,j);
if (path[i][j]+1!=j)
cout << ")";
}
}
int main()
{
int zz=1;
while(scanf("%d",&n)!=EOF&&n!=0)
{
Init();
MatrixChain();
printf("Case %d: (",zz++);
TraceBack(1,n);
printf(")\n");
}
return 0;
}
相关文章推荐
- zju/zoj 1276 Optimal Array Multiplication Sequence(dp)
- ZOJ 1276 Optimal Array Multiplication Sequence
- 矩阵连乘积 ZOJ 1276 Optimal Array Multiplication Sequence
- ZOJ 1276 Optimal Array Multiplication Sequence
- ZOJ 1276 Optimal Array Multiplication Sequence
- UVA 348 & ZOJ 1276 Optimal Array Multiplication Sequence(dp , 矩阵链相乘问题)
- zoj 1276 Optimal Array Multiplication Sequence
- ZOJ1276 Optimal Array Multiplication Sequence dp
- ZOJ 1276 Optimal Array Multiplication Sequence(矩阵连乘)
- Optimal Array Multiplication Sequence+uva+简单区间dp
- UVA 348 Optimal Array Multiplication Sequence 最优矩阵链乘 dp经典
- uva 348 - Optimal Array Multiplication Sequence
- uva 348 Optimal Array Multiplication Sequence
- uva 348 Optimal Array Multiplication Sequence
- UVA - 348 Optimal Array Multiplication Sequence
- uva 348 - Optimal Array Multiplication Sequence
- UVa 348 Optimal Array Multiplication Sequence (DP 最优矩阵链乘)
- UVa 348 - Optimal Array Multiplication Sequence
- uva348 - Optimal Array Multiplication Sequence(简单动规)
- uva348 Optimal Array Multiplication Sequence