BZOJ 1933 Shoi2007 Bookcase 书柜的尺寸

1933: [Shoi2007]Bookcase 书柜的尺寸

Time Limit: 5 Sec Memory Limit: 64 MB
Submit: 638 Solved: 251
[Submit][Status][Discuss]

Description

Tom不喜欢那种一字长龙式的大书架，他只想要一个小书柜来存放他的系列工具书。Tom打算把书柜放在桌子的后面，这样需要查书的时候就可以不用起身离开了。显然，这种书柜不能太大，Tom希望它的体积越小越好。另外，出于他的审美要求，他只想要一个三层的书柜。为了物尽其用，Tom规定每层必须至少放一本书。现在的问题是，Tom怎么分配他的工具书，才能让木匠造出最小的书柜来呢？ Tom很快意识到这是一个数学问题。每本书都有自己的高度hi和厚度ti。我们需要求的是一个分配方案，也就是要求把所有的书分配在S1、S2和S3三个非空集合里面的一个，不重复也不遗漏，那么，很明显，书柜正面表面积（S）的计算公式就是：

由于书柜的深度是固定的（显然，它应该等于那本最宽的书的长度），所以要求书柜的体积最小就是要求S最小。Tom离答案只有一步之遥了。不过很遗憾，Tom并不擅长于编程，于是他邀请你来帮助他解决这个问题。

Input

文件的第一行只有一个整数n（3≤n≤70），代表书本的本数。接下来有n行，每行有两个整数hi和ti，代表每本书的高度和厚度，我们保证150≤hi≤300，5≤ti≤30。

Output

只有一行，即输出最小的S。

Sample Input

4
220 29
195 20
200 9
180 30

Sample Output

18000

HINT

Source

Day2

抄自cydiater大神

很容易看出来这是一道DP题，那么怎么设置状态就成了这道题的关键。本题有点特殊的地方是有两个维度的状态，而每个维度又有三个部分的参数，如果全部设置出来的话肯定会MLE。首先对书的厚度状态简化。

书的厚度是求和的，这个显然不能作为状态的值，作为状态的参数是比较好的， 30*70=2100 2100^3是内存无法接受的，简化状态求出前i本书的前缀和sum[i]，如果第一层的厚度是i,第二层的厚度是j，那么第三层的状态显然是sum[i]-j-k，bingo，内存的问题解决了。显然两个维度一个维度表状态另一个维度显然表示权值。但是问题是，书是三层的高度，不可能把每一层都暴力表示出来，所以巧妙的地方是把书的高度递减排序。然后这个问题就显然很好解决了。

#include <bits/stdc++.h>
#define ll long long
#define inf 10000010
#define eps 1e-7
using namespace std;
inline int read(){int x=0;int f=1;char ch=getchar();while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}return x*f;
}
const int MAXN=2105;
namespace zhangenming{struct node{int h,v;}a[80];inline bool mycmp(node xx,node yy){return xx.h>yy.h;}   int n,lim=0,f[2][MAXN][MAXN];int sum[80];int sta=0;void init(){n=read();memset(sum,0,sizeof(sum));for(int i=1;i<=n;i++){a[i].h=read();a[i].v=read();lim+=a[i].v;}sort(a+1,a+n+1,mycmp);for(int i=1;i<=n;i++){sum[i]=sum[i-1]+a[i].v;}}void dp(){memset(f,10,sizeof(f));f[sta][0][0]=0;int t,h;for(int i=1;i<=n;i++){sta^=1;memset(f[sta],10,sizeof(f[sta]));t=a[i].v;h=a[i].h;for(int j=0;j<=sum[i-1];j++){for(int k=0;k<=sum[i-1];k++){if(f[sta^1][j][k]>1000000) continue;if(j+k>sum[i-1]) continue;if(k==0) f[sta][j][t]=min(f[sta][j][t],f[sta^1][j][k]+h);else f[sta][j][k+t]=min(f[sta][j][k+t],f[sta^1][j][k]);if(j==0) f[sta][t][k]=min(f[sta][t][k],f[sta^1][j][k]+h);else f[sta][j+t][k]=min(f[sta][j+t][k],f[sta^1][j][k]);if(sum[i-1]-j-k==0) f[sta][j][k]=min(f[sta][j][k],f[sta^1][j][k]+h);else f[sta][j][k]=min(f[sta][j][k],f[sta^1][j][k]);}}}}void output(){int ans=10000000;for(int i=1;i<=lim;i++){for(int j=1;j<=lim;j++){if(i+j<sum[n]&&f[sta][i][j]<=10000){ans=min(ans,max(max(i,j),sum[n]-i-j)*f[sta][i][j]);}}}cout<<ans<<endl;}
}
int main(){using namespace zhangenming;init();dp();output();return 0;
}

转载于:https://www.cnblogs.com/something-for-nothing/p/7888831.html