USACO Section 5.3 Milk Measuring - DFSID＋DP...

这道题要是不要写具体方案数就很普通的可重复背包了..对一个bool数组进行DP...从1做到Pnum..每次当k空间的k-P[i]为可得到时..k空间则可得到...每次扫描空间从0到Q.最后当Q为true..这些物品各种可重复的组合说能得到Q...

　但要求具体的方案数就eggache了..开始我想用一个当在做重复背包的时候跟着判断并传递...每个空间不为bool...而是一个struct..包括到达这个空间的最小物品数..以及按字典序最小的那几个..也就是

　　　　struct node

{

int num,a[105];

}

应该是行得通的...空间来说Q最大20000...每个空间挂100的int数组..那么相当于开了个2000000的int 数组.在USACO允许的空间内..

　既然这节的TEXT就是写的搜索...我还是练了下DFSID...2分枚举当且查找方案的所用物品个数...然后搜能否得到解...而当判断有无解时就是一样的可重复背包..能背出Q就是可行解...不难发现若先搜较小的数更容易提前搜出答案...所以对P先排序...

　这里有个很强劲的剪枝...就是边往深层DFS一种方案时..当枚举的当且某物品单个大小前面就能背出来..那就没必要再以这个物品往下搜了..因为这个物品的单个价值前面能得到..那么它与所有的数怎么怎么样相加都不能对背包空间做出任何更新....

　还有个比较重要的剪枝..当搜(1,2)..显然得到的结果与搜(2,1)一样..所以在DFS搜一种方案时就要保证每一层所选择的物品必须在上一层的后面...

Program:

/*  
ID: zzyzzy12   
LANG: C++   
TASK: milk4
*/      
#include<iostream>      
#include<istream>  
#include<stdio.h>     
#include<string.h>      
#include<math.h>      
#include<stack>
#include<map>
#include<algorithm>      
#include<queue>   
#define oo 2000000005  
#define ll long long  
#define pi (atan(2)+atan(0.5))*2 
using namespace std;
int Q,N,P[105],num,ans[105];
int now[105],nownum;
bool can[30001];
bool DFS(int i,int k)
{
      int j;
      bool temp[30001];
      if (can[Q])
      {
            if (num>nownum)
            {
                  num=nownum;
                  for (j=1;j<=num;j++) ans[j]=now[j];              
            }     
            return true;          
      }
      if (!k) return false;
      for (;i<=N;i++)
      if (!can[P[i]])
      {
            for (j=0;j<=Q;j++) temp[j]=can[j];
            for (j=0;j<=Q-P[i];j++)
                if (can[j]) can[j+P[i]]=true;
            now[++nownum]=P[i];
            if (DFS(i+1,k-1)) return true;
            nownum--;
            for (j=0;j<=Q;j++) can[j]=temp[j];
      }
      return false;
}
bool ok(int k)
{
      if (k>N) return false;
      int i; 
      memset(can,false,sizeof(can));
      can[0]=true;
      nownum=0;
      if (DFS(1,k)) return true; 
      return false;
}
void getanswer()
{
      int l,r,mid;
      num=oo;
      r=1;
      while (!ok(r)) r*=2;
      l=r/2;
      while (r-l>1)
      {
            mid=(r+l)/2;
            if (ok(mid)) l=mid;
              else r=mid;
      }
      return;
}
int main()  
{  
      freopen("milk4.in","r",stdin);   
      freopen("milk4.out","w",stdout);
      scanf("%d%d",&Q,&N);
      int i;
      for (i=1;i<=N;i++) scanf("%d",&P[i]);
      sort(P+1,P+1+N);
      getanswer();
      printf("%d",num);
      for (i=1;i<=num;i++) printf(" %d",ans[i]);
      printf("\n");
      return 0;     
}