网络流 HDU 3605

建图   源点    ->   1024类人   ->   星球   ->     汇点

权             每类人数目       星球容量     星球容量

列举 0~1024  一位是1 那么和对应的星球建边

 #include<stdio.h>
 #include<algorithm>
 #include<string.h>
 #include<queue>
 #include<math.h>

 using namespace std;
 #define inf 100000000
 int z[];
 int vis[];
 int head[];
 int y[];
 int cnt,S,T;
 struct edg
 {
     int w,next,to;

 }x[];
 void add(int u,int v,int w)
 {
     x[cnt].next=head[u];
     x[cnt].to=v;
     x[cnt].w=w;
     head[u]=cnt++;
 }
 int bfs()
 {
     memset(vis,-,sizeof(vis));
     queue<int>q1;
     q1.push(S);
     vis[S]=;
     int j;

     while(!q1.empty())
     {
         int now=q1.front();
         q1.pop();
         for(j=head[now];j!=-;j=x[j].next)
         {
             if(vis[x[j].to]<&&x[j].w)
             {

                 vis[x[j].to]=vis[now]+;
                 q1.push(x[j].to);
             }
         }
     }
     return vis[T]!=-;
 }
 int dfs(int u,int w)
 {
     int ans=;

     if(u==T)
         return w;
     int j;
     for(j=head[u];j!=-;j=x[j].next)
     {
         if(vis[x[j].to]==vis[u]+&&x[j].w)
         {

             int b=dfs(x[j].to,min(w-ans,x[j].w));
             ans=ans+b;
             x[j].w-=b;
             x[j^].w+=b;
         }
     }
     return ans;
 }
 int main()
 {
     int n,m;

     while(scanf("%d%d",&n,&m)!=EOF)
     {
         int i,j;
         memset(z,,sizeof(z));

         for(i=;i<=n;i++)
         {
             int b=;

             for(j=;j<=m;j++)
             {
                 int a;
                 scanf("%d",&a);
                 if(a==)b=b+(<<(j-));
             }
             z[b]++;
         }
         for(i=;i<=m;i++)
             scanf("%d",&y[i]);
         S=;
         T=;  //这2个点只要没用过就行
         cnt=;
         memset(head,-,sizeof(head));
         for(i=;i<=;i++)
         {
             if(z[i])
             {
                 add(S,i,z[i]);
                 add(i,S,);
             }
         }

         //0 1024  人的下标  1025 1035 星球下标
         for(i=;i<=;i++)
         {
             if(z[i])
             {
                 for(j=;j<=m;j++)
                 {
                     if(<<(j-)&i)
                     {
                         add(i,+j,y[j]);
                         add(+j,i,);
                     }
                 }
             }
         }

         for(j=;j<=m;j++)
         {
             add(+j,T,y[j]);
             add(T,+j,);
         }

         int ans=;
         while(bfs())
             ans+=dfs(S,inf);
         if(ans>=n)
             printf("YES\n");
         else
             printf("NO\n");
     }

     return ;
 }