uva437 动态规划

题目链接点击打开链接

题意：有n（n<=30）种立方体，每种都有无穷多个。要求选一些立方体落成一个尽量高的柱子（可自行选择哪一条边作为高） ，使得每个立方体的底面长宽分别严格小于它下方立方体的底面长度。

分析：用二元组（a , b）来表示“顶面尺寸为a*b的状态” ， 每增加一个长方体后顶面的长宽都会严格减小，所以这个图是DAG，套用DAG最长路算法。因为a,b值过大，所以不能直接用d（a,b）表示状态值。  用（idx , k）来“间接”表示这个状态， idx 为顶面立方体的序号，k是高的序号。

状态转移方程   d(id1 , k1) = max{ d(id2 , k2) + a[id1][k2] , ( (id1,k1)  , (id2 , k2) ) 属于E} ， d[i][j] 初始值为 a[i][j] ;

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
using namespace std;
const int maxn = 35 ;
int d[maxn][3] , a[maxn][3] , G[maxn][3][maxn][3] , visit[maxn][3];
int n ;

int dp(int id1 , int k1)  //记忆化搜索，用visit数组记录有没有访问过比较方便
{
    int &ans = d[id1][k1] ;
    if(visit[id1][k1]) return ans ;
    else
    {
        for(int id2 = 1 ; id2 <= n ; id2++ )
        {
            for(int k2 = 0 ; k2 <3 ; k2 ++)
                if(G[id1][k1][id2][k2]) //有边的话，
                  ans = max(ans , dp(id2 , k2) + a[id1][k1]) ; //状态转移
        }
        visit[id1][k1] = 1 ;
    }
    return ans ;
}

void Go(int i , int j)  //判断是否有边 ，三种情况下有边
{
   for(int x = 0 ; x < 3 ; x ++)
   {
       for(int y = 0 ; y < 3 ; y++)
       {
           if(a[i][(x + 1)%3] < a[j][(y + 1) % 3]  && a[i][(x+2)%3] < a[j][(y+2)%3]) 
              G[i][x][j][y] = 1 ;
           if(a[i][(x + 2)%3] < a[j][(y + 1) % 3]  && a[i][(x+1)%3] < a[j][(y+2)%3])
              G[i][x][j][y] = 1 ;
           if(a[i][(x + 1)%3] < a[j][(y + 2) % 3]  && a[i][(x+2)%3] < a[j][(y+1)%3])
              G[i][x][j][y] = 1 ;
       }
   }
}

void Print()
{
    for(int i = 1 ; i <= n ; i ++)
    {
        for(int x = 0 ; x <3 ; x ++)
        {
            for(int j = 1 ; j <= n ; j++)
            {
                for(int y = 0 ; y < 3 ; y ++)
                    if(G[i][x][j][y])
                    printf("i=%d x=%d j=%d y=%d %d\n" , i , x , j , y ,G[i][x][j][y] ) ;
            }
        }
    }
}

void Build()
{
    memset(G , 0 , sizeof(G)) ;
    for(int i = 1 ; i <= n ; i ++)
    {
        for(int j = 1 ; j <= n ; j ++)
        {
           Go(i , j) ;
        }
    }
}

int main()
{
    //freopen("a.txt" , "r" , stdin) ;
    int kase = 0  ;
    while(scanf("%d" , &n) != EOF && n)
    {
        for(int i = 1 ; i <= n ; i ++)
            scanf("%d%d%d" , &a[i][0] , &a[i][1] , &a[i][2]) ;
        Build() ;
        //Print() ;
        memset(visit , 0 , sizeof(visit)) ;
        for(int i = 1 ; i <= n ; i++)  //初始化，注意值为a[i][j] ，不为0
          for(int j = 0 ; j <3 ; j ++)
            d[i][j] = a[i][j] ;
        for(int i = 1 ; i <= n ; i++)
        {
            for(int j = 0 ; j <3 ; j ++)
                dp(i , j) ;
        }
        int Max = 0 ;
        for(int i = 1 ; i <= n ; i ++) //求最大值
        {
            for(int j = 0 ; j < 3 ; j ++)
                Max = max(Max , d[i][j]) ;
        }
        printf("Case %d: maximum height = %d\n" , ++kase , Max) ;
    }
    return 0;
}