所有题在A不掉之前都是难题,而当你把它A掉后,它便成了水题,而在你A掉它之前,你永远不知道它有多水。
这道题,我最先想到的是DP,觉得就是一个典型背包问题,不过是背包多了一点而已,然而再仔细一想N=50便意味着50个背包,128^50的运算量,不用细算也知道这样的程序跑出来太阳都熄火了,所以这样的DP是走不通的。想到这个section的主题是讲搜索优化的,那就试试搜索吧,直接DFS加剪枝?嗯,想了很久没想出来怎么做DFS,就在这个时候ZZY看完了刚更新的日剧,坐在旁边摇头晃脑满脸的意犹未尽,时不时往我这边漂上几眼,看到我一副苦逼的表情,点了几下头然后便笑而不语,过了许久才阴森的冒出一句话“这个要用DFSID”,我靠,什么是DFSID啊,完全不懂啊,听起来好高端啊...赶紧google之.....
搜索DFSID,wiki上没有,××没有,看了别人的博客,大致的感觉就是先确定一个界限再做DFS,而在这道题里面感觉倒有点像是枚举加DFS。在这个题中,用DFS找到能切割出的最大rail数不好怎么下手,但是,用DFS确定是否能切割出K块rail是相对容易的,所以我们就可以通过枚举出各个K值然后DFS得出能切割出的最大rail数,这里可以用二分的思想来做优化。
现在的主要问题就是如何DFS了。 既然是要确定一堆board能否切割出k块rail,我们不妨用一个全局的数组来储存各个board 的长度,用bool DFS(int k)来表示k块rail能否被切割出来,如果这堆board能切割出k块rail那么肯定能切割出k块最短的rail,所以,我们应该先给rail排个序。
朴素DFS代码:
bool DFS(int k)
{
if(k<0) return true;
for(int i=0;i<boards;i++)
{
if(remain[i]>=rarr[k])//remain[i]记录了第i块board的剩余长度
{
remain[i]-=rarr[k];
if(DFS(k-1)) return true;//先对rail排了序,rail从大到小搜
remain[i]+=rarr[k];
}
}
return false;
}
这样没有做任何优化的DFS 只能过3个点,所以动手优化吧。
在网上看了大牛们的博客,都说对于这样的切割问题,用剩余材料做优化可以减少大量的冗余搜索,其大致思路是这样的,如果boad的总长度是board_sum,需要切割出来的rail总长度是rail_sum,那么最大的可浪费材料是max_waste,如果某一种切割方式在切割完之前已产生了waste>max_waste的浪费,那么显然这种方式是不可行的。
优化后代码
bool DFS(int k) { if(k<0) return true;
//用浪费材料做剪枝
int waste=0;
for(int i=0;i<boards;i++) if(remain[i]<rarr[0]) waste+=remain[i]; if(waste>max_waste) return false; for(int i=start;i<boards;i++) { if(remain[i]>=rarr[k]) { if(k-1>=0 && rarr[k]==rarr[k-1]) ns=i; remain[i]-=rarr[k]; if(DFS(k-1)) return true; remain[i]+=rarr[k]; } } return false; }
优化后,竟然还是只能过3个点,太受打击了。由于rail最多有1023块,而rail的长度是<=128的正整数,之意味着很多rail的长度是一样的,而对于长度一样的rail,它们的顺序是不重要的,所以可以做如下优化:
优化后代码
优化后提交,果断AC。bool DFS(int k,int start) { if(k<0) return true; int waste=0; for(int i=0;i<boards;i++) if(remain[i]<rarr[0]) waste+=remain[i]; if(waste>max_waste) return false; int ns=0; for(int i=start;i<boards;i++) { if(remain[i]>=rarr[k]) { if(k-1>=0 && rarr[k]==rarr[k-1]) ns=i; else ns=0; remain[i]-=rarr[k]; if(DFS(k-1,ns)) return true; remain[i]+=rarr[k]; } } return false; }
经实验发现,这两种优化缺少任何一种都只能过3个点,只有两种优化同时采用才能AC。
下面是完整的源代码:
#include<cstdio> #include<algorithm> using namespace std; FILE *in,*out; int barr[60],remain[60],rarr[1050],rarr_sum[1050],max_waste,board_sum,boards; int search(int s,int e); bool check(int k); bool DFS(int k,int start); int main() { in=fopen("fence8.in","r"); out=fopen("fence8.out","w"); int rails,i; fscanf(in,"%d",&boards); for(i=0;i<boards;i++) fscanf(in,"%d",&barr[i]); for(i=0;i<boards;i++) barr[i]=-barr[i]; sort(barr,barr+boards); for(i=0;i<boards;i++) barr[i]=-barr[i]; fscanf(in,"%d",&rails); for(i=0;i<rails;i++) fscanf(in,"%d",&rarr[i]); sort(rarr,rarr+rails); for(i=0;i<boards;i++) board_sum+=barr[i]; for(i=0;i<rails;i++) rarr_sum[i]=rarr[i]; for(i=1;i<rails;i++) rarr_sum[i]+=rarr_sum[i-1]; max_result=0; for(i=1;i<rails;i*=2) if(!check(i)) break; if(i>rails) i=rails; fprintf(out,"%d\n",search(i>>1,i)); fclose(in); fclose(out); return 0; } int search(int s,int e) { int mid=(s+e)>>1; if(e-s==1) mid=e; bool flag=check(mid); if(e-s==1 && flag) return e; if(e-s==1 && !flag) return s; if(e==s && flag) return e; if(flag) return search(mid,e); else return search(s,mid); } bool check(int k) { if(k<=0) return true; for(int i=0;i<boards;i++) remain[i]=barr[i]; max_waste=board_sum-rarr_sum[k-1]; if(max_waste<0) return false; return DFS(k-1,0); } bool DFS(int k,int start) { if(k<0) return true; int waste=0; for(int i=0;i<boards;i++) if(remain[i]<rarr[0]) waste+=remain[i]; if(waste>max_waste) return false; int ns=0; for(int i=start;i<boards;i++) { if(remain[i]>=rarr[k]) { if(k-1>=0 && rarr[k]==rarr[k-1]) ns=i; else ns=0; remain[i]-=rarr[k]; if(DFS(k-1,ns)) return true; remain[i]+=rarr[k]; } } return false; }