dp

如果我们直接定义状态：

\(dp[i][t1][t2][t3][h1][h2][h3]\)表示前i个，第一层宽度为t1,，第二层宽度为t2，第三层宽度为t3，第一层高度为h1，第二层高度为h2，第三层高度为h3的最小面积。

如果直接这样定义，你会发现，你不仅内存炸飞，时间也会T的飞起。

考虑优化状态。

1.首先，你会发现，面积可以直接用t1,t2,t3,h1,h2,h3算出来，所以我们不妨砍掉一维

2.列一波状态转移方程，你会发现，i只会从i-1转移过来，于是又可以把第一维滚动

3.不难发现，\(t1+t2+t3=\sum _{j=1}^{j \le i} t_j\)于是只用知道t1，t2，t3中的任意两个，就可以推出第三个

那么，状态就优化成了:\(dp[0/1][t1][t2][h1][h2]\)表示前i个，第一层宽度为t1,，第二层宽度为t2，第一层高度为h1，第二层高度为h2，的第三层最小高度。

然后，你又会发现，内存和时间依旧承受不住。。。

仔细琢磨一下，不难观察到，每一层的高度是这一层中所放书本的最大值，那么如果按照一定顺序插入书本，高度不就可以省略掉了吗？

因此，我们先把书本按高度从大到小排个序，这样每层的高度就是第一次插入到这层书的高度，于是状态又优化成了:\(dp[0/1][t1][t2]\)表示前i个，第一层宽度为t1，第二层宽度为t2，的最小三层总高度

转移的话就不用多说了吧。。。

代码：

注意一点，每层至少得有1本书，因此最后算面积时要排除某一层没有书的情况。（第二组样例已经良心地说明了这点）

#include<bits/stdc++.h>
using namespace std;
int n,dp[2][2110][2110],sum[101],ans=(1<<30);
struct node{int t,h;
}A[1010];
bool cmp(node A,node B){return A.h>B.h;
}
int main(){scanf("%d",&n);for(int i=1;i<=n;i++)scanf("%d %d",&A[i].h,&A[i].t);memset(dp,63,sizeof(dp));const int oo=dp[0][0][0];sort(A+1,A+1+n,cmp);for(int i=1;i<=n;i++)sum[i]=sum[i-1]+A[i].t;dp[0][0][0]=0;for(int i=1;i<=n;i++){int now=i&1,la=!now;memset(dp[now],63,sizeof(dp[now]));for(int j=0;j<=sum[i-1];j++){for(int k=0;k<=sum[i-1];k++){int o=sum[i-1]-j-k;if(o<0)continue;if(dp[la][j][k]==oo)continue;if(j==0)dp[now][j+A[i].t][k]=min(dp[now][j+A[i].t][k],dp[la][j][k]+A[i].h);else dp[now][j+A[i].t][k]=min(dp[now][j+A[i].t][k],dp[la][j][k]);if(k==0)dp[now][j][k+A[i].t]=min(dp[now][j][k+A[i].t],dp[la][j][k]+A[i].h);else dp[now][j][k+A[i].t]=min(dp[now][j][k+A[i].t],dp[la][j][k]);if(o==0)dp[now][j][k]=min(dp[now][j][k],dp[la][j][k]+A[i].h);else dp[now][j][k]=min(dp[now][j][k],dp[la][j][k]);}}}for(int i=1;i<=sum[n];i++){for(int j=1;j<=sum[n];j++){int o=sum[n]-i-j;if(o<=0)continue;if(dp[n&1][i][j]==oo)continue;ans=min(ans,max(max(i,j),o)*dp[n&1][i][j]);}}cout<<ans;return 0;
}