USACO4.1 Fence rails DFSID

所有题在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，它们的顺序是不重要的，所以可以做如下优化：
优化后代码

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;
}
 优化后提交，果断AC。

经实验发现，这两种优化缺少任何一种都只能过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;
}