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;
}