dp_1_M

最近开始刷奇奇神的dp专题,呃,我是弱菜,啥都不会,现在才开始刷dp,
在nocLyt神的熏陶下,觉着区间dp有些感觉了
不过还是调半天才调出来
今天做的这个
M - Optimal Array Multiplication Sequence
Time Limit:3000MS     Memory Limit:0KB     64bit IO Format:%lld & %llu

Description


 Optimal Array Multiplication Sequence 

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 tex2html_wrap_inline67 array, and B is a tex2html_wrap_inline71 array, it will take tex2html_wrap_inline73 , 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 tex2html_wrap_inline103 array, Y is a tex2html_wrap_inline67 array, and Z is a tex2html_wrap_inline111 array. Let's look at the number of multiplications required to compute the product using the two different sequences:

(X Y) Z

  • tex2html_wrap_inline119 multiplications to determine the product (X Y), a tex2html_wrap_inline123 array.
  • Then tex2html_wrap_inline125 multiplications to determine the final result.
  • Total multiplications: 4500.

X (Y Z)

  • tex2html_wrap_inline133 multiplications to determine the product (Y Z), a tex2html_wrap_inline139 array.
  • Then tex2html_wrap_inline141 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 multiplications 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 tex2html_wrap_inline157 . 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))

虽说这个是水题吧,但好歹是自己用心写的一个dp,就放在这了,

#include<stdio.h>
#include
<string.h>
#include
<math.h>
#define inf 0x7ffffff
#define maxn 15
int f[maxn][maxn],path[maxn][maxn];
int l[maxn],r[maxn];
int n;
int min(int a,int b)
{
    
return a<b?a:b;
}
void print(int h,int t)
{
    printf(
"(");
    
if (t-h==1)
    {
        printf(
"A%d x A%d",h,t);
    }
    
else
    {
        
int tmp=path[h][t];
        
if(tmp-h==0)
        {
            printf(
"A%d x ",h);
             
if(t-tmp-1==0) printf("A%d",t);else print(tmp+1,t);
        }
        
else if(t-tmp==0)
        {
            
            
if(h-tmp-1==0) printf("A%d",h);else print(h,tmp-1);
            printf(
" x A%d",t);
        }
        
else
        {
            
if(h-tmp==0) printf("A%d",h);else print(h,tmp);
            printf(
" x ");
            
if(t-tmp-1==0) printf("A%d",t);else print(tmp+1,t);
        }
    }
    printf(
")");
}
int main()
{
    
int i,j,k,times;
    times
=0;
    
while(scanf("%d",&n)!=EOF&&n!=0)
    {
        times
++;
        
for(i=1; i<=n; i++) scanf("%d%d",&l[i],&r[i]);
        memset(f,
0,sizeof(f));
        
for(i=1; i<=n-1; i++)
            f[i][i
+1]=l[i]*r[i]*r[i+1];//,printf("%d %d %d\n",i,i+1,f[i][i+1]);
        for(i=n-2; i>=1; i--)
        {
            
for(j=i+2; j<=n; j++)
            {
                f[i][j]
=inf;
                
for(k=i; k<=j; k++)
                    
if(f[i][k]+f[k+1][j]+l[i]*r[k]*r[j]<f[i][j])
                    {
                        
//printf("%d %d %d %d %d %d\n",i,k,j,f[i][k],f[k+1][j],l[i]*r[k]*r[j]);
                        f[i][j]=f[i][k]+f[k+1][j]+l[i]*r[k]*r[j];
                        path[i][j]
=k;
                    }
                
//printf("%d %d %d %d\n",i,j,f[i][j],path[i][j]);
            }
        }
        
//printf("%d\n",f[1][n]);
        printf("Case %d: ",times);
        print(
1,n);
        printf(
"\n");
    }
    
return 0;
}

呃,我是water,那个这个还可以用记忆化dp来实现,不过我没写过很好的记忆化搜索,干脆每次都写的递推实现
不过效率还好……

posted on 2012-06-02 12:25 jh818012 阅读(124) 评论(0)  编辑 收藏 引用


只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   博问   Chat2DB   管理


<2024年12月>
24252627282930
1234567
891011121314
15161718192021
22232425262728
2930311234

导航

统计

常用链接

留言簿

文章档案(85)

搜索

最新评论