题目链接点击打开链接
题意:有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;
}