标签:
Time Limit: 1000MS | Memory Limit: 10000K | |
Total Submissions: 18103 | Accepted: 6439 |
Description
Input
Output
Sample Input
100
2
Sample Output
68
Hint
Source
又见花式大剪枝。
1、从半径=n+1开始搜
2、当前已有表面积,加上之后层的预估最小表面积,若大于最优解,减掉。
3、当前已有体积,加上之后层的预估最小体积,若大于最优解,减掉。
4、搜索途中,若体积超出限制,减掉(大概不会生效)
只用上面这些,只能拿到70分(在洛谷测的)
然后遍寻题解,找到了神奇的剪枝5:(目标体积-已有体积)/r*2+已有表面积>=ans,加上以后0ms AC,简直可怕。
1 /*by SilverN*/ 2 #include<iostream> 3 #include<algorithm> 4 #include<cstring> 5 #include<cstdio> 6 #include<cmath> 7 using namespace std; 8 const int mxn=30; 9 int ms[mxn],mv[mxn]; 10 int n,m; 11 int ans; 12 void DFS(int smm,int v,int dep,int h,int r){//表面积 体积 层数 高 半径 13 if(v>n)return; 14 if(!dep){ 15 if(v==n && smm<ans)ans=smm; 16 return; 17 } 18 if(v+mv[dep-1]>n || smm+ms[dep-1]>ans || (n-v)/r*2+smm>=ans)return;//剪枝 19 for(int i=r-1;i>=dep;--i){//半径 20 if(dep==m)smm=i*i;//上表面的总面积 21 int mxh=min((n-v-mv[dep-1])/(i*i),h-1); 22 for(int j=mxh;j>=dep;--j){//高度 23 DFS(smm+2*i*j ,v+i*i*j ,dep-1,j,i); 24 } 25 } 26 return; 27 } 28 void init(){ 29 for(int i=1;i<=20;i++){ 30 ms[i]=ms[i-1]+i*i*2;//最小侧面积 31 mv[i]=mv[i-1]+i*i*i;//最小体积 32 } 33 return; 34 } 35 int main(){ 36 init (); 37 int i,j; 38 scanf("%d%d",&n,&m); 39 ans=1e9; 40 DFS(0,0,m,n+1,n+1); 41 if(ans==1e9)ans=0; 42 printf("%d\n",ans); 43 return 0; 44 }
标签:
原文地址:http://www.cnblogs.com/SilverNebula/p/5971785.html