ZJU-1093

 

DP，先對輸入進行處理，每一個立方體處理爲3個面，如有尺寸雷同面，則合併。
保證每個面的x<y。

然後爲了減小無效搜索的次數，對於面進行了排序，以y爲第一關鍵字，x爲第二關鍵字，都是用quicksort，發現快速排序用起來真的是挺方便，不過不要排序直接DP應該也是沒有問題的。

（當然merge sort還能更快，就是代碼太多了。值得一提的是儘管merge sort和heap sort的最壞情況的複雜度都是比較排序的極限O(nlgn)但是其中隱含的參數還是有些區別的，merge sort要更快一些。）

排完了序以後，能放在第i個面上的面必定在0～i-1中，maxdone[]這個數組存儲已經計算出的每個面上能包括該面能壘的最高高度，所以符合DP的要求
1.最優子結構：如果第i個面已經是最高高度，則它上面的i-1這個面以上也是i-1這個面可以壘的最高高度。
2.重疊子問題：重疊的子問題的結果存儲在maxdone[]用來加速計算。

//2008-01-11 15:35:06 Accepted 1093 C++ 00:00.00 436K #include<stdio.h> struct piles{ int a,b; int h; }; struct piles pl[91]; int np; int maxdone[91]; int partb(int p,int r) { int x = pl[r].b; int i = p-1,j; struct piles t; for(j=p;j<=r-1;j++) { if(pl[j].b<=x) { i++; t=pl[j];pl[j]=pl[i];pl[i]=t; } } t=pl[r];pl[r]=pl[i+1];pl[i+1]=t; return i+1; } int parta(int p,int r) { int x = pl[r].a; int i = p-1,j; piles t; for(j=p;j<=r-1;j++) { if(pl[j].a<=x) { i++; t=pl[j];pl[j]=pl[i];pl[i]=t; } } t=pl[r];pl[r]=pl[i+1];pl[i+1]=t; return i+1; } void quicksort(int p,int r,int mode) { if(p<r) { int q; if(mode==2) q=partb(p,r); else if(mode==1) q=parta(p,r); quicksort(p,q-1,mode); quicksort(q+1,r,mode); } } int dpsolve(int key) { if(maxdone[key]!=0) return maxdone[key]; else { int i; int maxre=0,t,re=0; re += pl[key].h; for(i=key-1;i>=0;i--) { if(pl[i].a<pl[key].a && pl[i].b<pl[key].b) { t = dpsolve(i); if(t>maxre) maxre = t; } } maxdone[key] = re+maxre; return maxdone[key]; } } int main() { freopen("1093.txt","r",stdin); int n,x,y,z,d[3],i,j,k,t,min,count=0; int maxre,tre; while((scanf("%d",&n)!=EOF) && n!=0) { np=0; count++; for(k=0;k<n;k++) { scanf("%d %d %d",&d[0],&d[1],&d[2]); for(i=0;i<2;i++) { min=i; for(j=i+1;j<3;j++) { if(d[j]<d[min]) min = j; } if(min!=i) { t=d[i];d[i]=d[min];d[min]=t; } } x=d[0];y=d[1];z=d[2]; if(x==y && y==z) { pl[np].a = x;pl[np].b = y;pl[np].h = z;np++; } else if(x==y) { pl[np].a = x;pl[np].b = y;pl[np].h = z;np++; pl[np].a = x;pl[np].b = z;pl[np].h = y;np++; } else if(y==z) { pl[np].a = x;pl[np].b = y;pl[np].h = z;np++; pl[np].a = y;pl[np].b = z;pl[np].h = x;np++; } else { pl[np].a = x;pl[np].b = y;pl[np].h = z;np++; pl[np].a = x;pl[np].b = z;pl[np].h = y;np++; pl[np].a = y;pl[np].b = z;pl[np].h = x;np++; } } quicksort(0,np-1,2); k=0; for(i=1;i<np;i++) { if(pl[i].b==pl[k].b) continue; else { quicksort(k,i-1,1); k=i; } } quicksort(k,i-1,1); for(i=0;i<91;i++) { maxdone[i]=0; } maxre=0; for(i=0;i<np;i++) { tre = dpsolve(i); if(tre>maxre) maxre = tre; } printf("Case %d: maximum height = %ld/n",count,maxre); } fclose(stdin); return 0; }