codevs2777 栅栏的木料

题目描述 Description

农民John准备建一个栅栏来围住他的牧场。他已经确定了栅栏的形状，但是他在木料方面有些问题。当地的杂货储存商扔给John一些木板，而John必须从这些木板中找出尽可能多所需的木料。

当然，John可以切木板。因此，一个9英尺的木板可以切成一个5英尺和一个4英尺的木料 (当然也能切成3个3英尺的，等等)。John有一把（完美的）梦之锯，因此他在切木料时，不会有木料的损失。

所需要的栅栏长度可能会有重复（比如，一个3英尺和另一个3英尺长的栅栏可能同时都需要）。所需要的木料规格都已经给定。你不必切出更多木料，那没有用。

输入描述 Input Description

第1行: N (1 <= N <= 50), 表示提供的木板的数目

第2行到第N+1行: N行，每行包括一个整数，表示各个木板的长度。

第N+2行: R (1 <= R <= 1023), 所需木料的数目

第N+3行到第N+R+1行: R行，每行包括一个整数(1 <= ri <= 128)表示所需木料的长度。

输出描述 Output Description

只有一行，一个数字，表示能切出的最多的所需木料的数目。当然，并不是任何时候都能切出所有所需木料。

样例输入 Sample Input

4
30
40
50
25
10
15
16
17
18
19
20
21
25
24
30

样例输出 Sample Output

7

数据范围及提示 Data Size & Hint

USACO译题 4.1.2

 

解题思路：

一、一个发现：

我们很容易就可以发现若能满足a块木料，则必能满足最小的a块木料；所以当我们将木料排序之后，被切出的木料将在一个前缀区间里，所以我们发现这个答案是具有单调性的，所以我的第一直觉是二分答案，判断前a块木料是否能被切出。

那么如何判断呢？我的第一直觉就是贪心。

二、尝试贪心：

1、策略：我的贪心策略是把最大的木料给能切出它的最小的木板切，直到木料用尽或不存在这样的木板为止。

2、反例：可是当我兴高采烈把程序交上去时，才发现它只能得75分。（竟然得了75分！）反例来自当最大的木料与最小的木板之差小于最小的木料时，它会造成冗余的浪费。如：用6,7切出2,3,3,5，这样贪心会得到3，虽然实际上它显然是4。

3、继续尝试，我们却发现我们想不到更好的贪心了。。

三、DFS：

于是我又开始尝试暴力搜索，但是在DFS的时候却忘了一得到的结论，实际上它能很好地缩小搜索树的形态的。

我当时的做法是：暴力枚举每个木料，令其被N个木板切或不切，直到搜出最优解。

想到的一个最优化剪枝是若当前剩余的可用木板（不小于尚未被切的最小木料）的总和小于更新更优解所需的木料，则剪枝；这样的话，需要我们按照木料从小到大的数据搜索。

这样的话。。只能得67分，还不如贪心呢。

当然。。这其实也很有可能是因为我犯了一个错误，就是实际上，我每次判断当前剩余可用木板实际上都应用O(n)的时间复杂度扫一遍才对，但我在这里。。犯了一个比较奇葩的错误，导致实际上削弱了这个剪枝的强度。

四、对DFS的改进：

1、来自一的改进，将N个木板升序排序，从第一个木板开始搜索，不存在不被切的状态了，这样搜索树的深度就是答案了。

2、来自二的改进，先贪出一个较优解，再用这个较优解做最优化剪枝。

3、来自三中的剪枝可以得以继续延用。

4、来自题解的剪枝：我们（题解）发现，这个题木料一共有1024块，但却落在[1,128]∩N的范围内，这说明有很多块木料会是一样的，DFS的时候会算出它们的排列，这浪费了我们大量时间却没有任何意义。所以我们干脆让它有序好了，即令若干块长度相同的木料被切的木板序号不下降或不上升即可。

5、但是有了上述四个剪枝以后，我们还是会TLE，所以题解又提供了一个蛋疼剪枝，就是将每一个状态的木板也排序，对于若干个相同的木板，只试切第一个即可，即若令剪枝3中序列不下降，就切编号最小的那个；若令剪枝3中序列不上升，就切编号最大的那个。当然实际上这个剪枝之所以够强，我想主要是由于数据的原因。

那么有了以上五个剪枝，我们是可以A掉这道题了的，但尤其是第四个剪枝，难免有卡数据的嫌疑。实际上，我们还有真正优秀的算法，并不需要这种来自数据的剪枝。

五、迭代加深搜索：

1、让我们继续沿用一中的思路，如果是二分+DFS判定呢？对！这样一来，搜索最优解就转变成了搜索可行解，这就为我们提供了新的优化方向。

2、四中改进1,2,3,4显然可以继续沿用，当然如果已贪出较优解就不需要再二分了，但这还是不够的，那么更牛的优化来自哪里呢？是来自搜索顺序的，搜索顺序对搜索最优解并不能造成太大的影响，但对于搜索可行解可就大不一样了，改变搜索顺序往往能起到非常惊人的效果。比如本题中，先切大木料显然要比先切小木料更优！这样的话我们发现剪枝3可以做某种程度的改动了，即其并不再需要每次都O(n)的扫描了，而是再切的时候判断剩余木料和最小木料的大小即可，因为最小木料一定是最后被切出来的。这样，究竟是使它增强了还是减弱了呢。。。还真不好说，毕竟搜索顺序都改变了。

3、但是这里同样需要注意一个问题，就是再搜索可行解的时候，需要中间强行return，这时千万不要忘了回溯，我在这个地方蛋疼了好久。

综，至此我们已经非常完美的解决本题了，即使是最大数据，运行时间也在0.004s.

#include<iostream>
using namespace std;
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
int ans,N,M,s[],a[],b[],sum;
inline bool dfs(int x,int pred){
    if(!x)return ;
    if(sum<s[x])return ;
    sum-=b[x];
    for(int i=x<ans&&b[x]==b[x+]?pred:;i<N;++i){
        if(i<N&&a[i]==a[i+])continue;
        if(a[i]>=b[x]){
            a[i]-=b[x];
            if(a[i]<b[])sum-=a[i];
            if(dfs(x-,i)){
                if(a[i]<b[])sum+=a[i];
                sum+=b[x];
                a[i]+=b[x];
                return ;
            }
            if(a[i]<b[])sum+=a[i];
            a[i]+=b[x];
        }
    }
    sum+=b[x];
    return ;
}
inline bool check(short m){
    short x,minx,j,i,tmp[];
    for(i=;i<N;++i)tmp[i]=a[i];
    for(i=m;i;--i){
        x=-,minx=0x7fff;
        for(j=;j<N;++j)
            if(tmp[j]>=b[i]&&tmp[j]<minx){
                minx=tmp[j];
                x=j;
            }
        if(x>-)tmp[x]-=b[i];
        else return ;
    }
    return ;
}
char * ptr=(char *)malloc();
inline void in(int &x){
    while(*ptr<''||*ptr>'')++ptr;
    x=;
    while(*ptr>&&*ptr<)x=x*+*ptr++-'';
}
int main(){
    fread(ptr,,,stdin);
    in(N);
    int i;
    for(i=;i<N;++i)in(a[i]),sum+=a[i];
    in(M);
    for(i=;i<=M;++i)in(b[i]);
    sort(b+,b+M+);
    sort(a,a+N);
    b[M+]=-;
    for(i=;i<=M;++i)s[i]=s[i-]+b[i];
    short l=,r=M;
    if(check(r)){
        printf("%hd",r);
        return ;
    }
    while(r-l>){
        short m=(l+r)>>;
        if(check(m))l=m;
        else r=m;
    }
    ans=r;
    while(ans<=M&&dfs(ans,N))++ans;
    printf("%d",ans-);
}