Packing Rectangles

  这题应该算是chapter 里最难的一题了，个人是看了之后表示一点想法都没有，告诉我是搜索的我都不知道怎么去搜，后面终于下决心要A掉它，网上找了好多的解题报告，看了它们的思路，我个人对它们的正确性是没有疑问，就是为什么这样做就包含了所有的情况呢？后来还是自己独自想通的。大家注意我的分析过程都是分类穷举，因此没有遗漏，不过可能也还有不正确的地方，，希望大家能够指出。

 这里和大家分享一下,希望大家能尽快的形成自己的思路。不要像我一样无限纠结。

网上对这题大部分映像不是很好，个人感觉这题对我们的逻辑思维还是很有挑战性的。

 

大致的思路：

 

一、首先确定矩形的位置

 

1.四个并排 

 

 

 

2.下面一个，上面三个(下面三个，上面一个）

 

下面一个和上面三个的组合情况有三种

• 1-3（图1）
• 1-2（图2）
• 1-1（图3，图4）因此3和4本质是一样的。
• 

3.上下各两个

 

 

之所以这里要讨论是因为最小矩形（答案）的宽是受高的限制的。

 

不会做图，就只些分析过程了：

0,1,2,3分别表示从左至右，从上至下的编号。

有四种情况

• h[3]<h[2]

 1.h[3]+h[1]<=h[2]

 2.h[3]+h[1]>h[2] &&h[1]+h[3]<=h[0]+h[2]

• h[3]>h[2]

3.h[1]+h[3]<=h[0]+h[2]

4.h[3]>h[0]+h[2]

 

 

 

 

二、每一种位置方式都会有4!种矩形的的摆放

三、每一种摆放再枚举宽高的组合有2^4种

 

实现过程时可以先执行二、三步骤（得到矩形的摆放，以及宽高的组合），再执行步骤一（枚举出现的位置） 

 

/* 让我卡的第十七组数据！！！就一个小小的错误让我纠结了差不多一天。纪念一下。 4 4 5 5 3 4 6 3 ans: 72 8 9 */ #include<cstdio> #include<iostream> #include<algorithm> #include<string.h> #include<climits> using namespace std; #define max(x,y) (x>y?x:y) #define M 24*16*6//4!*2^4*(6case) #define N 4 struct rect{ int w,h; }ans[M];int cnt; bool cmp(const rect &a,const rect &b){ return a.w<b.w; } int org[N][2],comR[N],mini; int comH[4];//各个矩形高的序号的组合 int w[5],h[5];//w0,w1,w2,w3,w,h0,h1,h2,h3,h; void getw_h(){ int i; for(i=0;i<N;i++){ w[i]=org[ comR[i] ][ comH[i]^1 ]; h[i]=org[ comR[i] ][ comH[i] ]; } } void getAns(){ if(mini>w[N]*h[N])mini=w[N]*h[N],cnt=0; if(mini>=w[N]*h[N]){ if(w[N]>h[N])swap(w[N],h[N]); ans[cnt].h=h[N];ans[cnt].w=w[N]; cnt++; } } void case0(){ w[N]=w[0]+w[1]+w[2]+w[3]; h[N]=max( max(h[0],h[1]) , max(h[2],h[3]) ); } void case1(){ w[N]=max( w[0]+w[1]+w[2],w[3] ); h[N]=max( max(h[0],h[1]) ,h[2] )+h[3]; } void case2(){ w[N]=max( w[0]+w[1],w[3] )+w[2]; h[N]=max( max(h[0],h[1])+h[3], h[2] ); } void case3(){ w[N]=max(w[0],w[3])+w[1]+w[2]; h[N]=max ( max( h[0]+h[3],h[1] ) , h[2] ); } void case4(){ h[N]=max( h[0]+h[2] , h[1]+h[3] ); if(h[2]>=h[1]+h[3]){ w[N]=max (max (w[1],w[3]) + w[2] , w[0] ); } else if(h[3]<h[0]+h[2]){ if(h[3]<h[2]){ w[N]=max( w[0],w[2] )+w[1];//这种情况比较特殊!!! w[N]=max( w[N],w[2]+w[3] ); } else{ w[N]=max( w[0],w[2] )+w[3]; w[N]=max( w[N],w[0]+w[1] ); } } else{ w[N]=max( max(w[0],w[2])+w[3],w[1] ); } } void comHeight(int k){//高度组合 if(k==N){ getw_h(); case0();getAns(); case1();getAns(); case2();getAns(); case3();getAns(); case4();getAns(); } else{ for(int i=0;i<2;i++){ comH[k]=i; comHeight(k+1); } } } bool used[N]; void comRect(int k){//矩形排列 if(k==N){ comHeight(0); } else{ for(int i=0;i<N;i++){ if(!used[i]){ comR[k]=i; used[i]=true; comRect(k+1); used[i]^=1; } } } } int main(){ freopen("packrec.in","r",stdin); freopen("packrec.out","w",stdout); int i,j; for(i=0;i<N;i++){ for(j=0;j<2;j++){ cin>>org[i][j]; } } mini=INT_MAX; cnt=0; memset(used,0,sizeof(used)); comRect(0); sort(ans,ans+cnt,cmp); cout<<mini<<endl; cout<<ans[0].w<<' '<<ans[0].h<<endl; for(i=1;i<cnt;i++)if(ans[i].w!=ans[i-1].w)cout<<ans[i].w<<' '<<ans[i].h<<endl; return 0; }