HDOJ 3605 - Escape 状态压缩缩点+最大流

                题意:

                         有N(<=100000)个人要逃离到M(<=10)个星球..每个人有喜欢的星球与不喜欢的星球..而每个星球的容量是有限的..问能否让所有人都在自己喜爱的星球上..

                题解:

                         上次使用二分图多重匹配..类匈牙利算法做的.这次用最大流来做...构图很easy了...但是...有100000个人..直接搞会超时...要注意.船至多只有10艘..那么最多对应了2^10=1024种状态...好了..那么100000个人就可以压缩成最多1024个点了..跑最大流~AC

Program:

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<queue>
#define MAXN 2005
#define MAXM 20000005
#define oo 1000000007
#define ll long long
using namespace std;
int sum[1200];    
struct Dinic  
{  
      struct node   
      {  
             int x,y,c,next;  
      }line[MAXM];     
      int Lnum,_next[MAXN],dis[MAXN],Q[MAXN],Qnum;  
      void initial(int n)   
      {  
             for (int i=0;i<=n;i++) _next[i]=-1;  
             Lnum=-1;  
      }   
      void addline(int x,int y,int c)  
      {  
             line[++Lnum].next=_next[x],_next[x]=Lnum;  
             line[Lnum].x=x,line[Lnum].y=y,line[Lnum].c=c;  
             line[++Lnum].next=_next[y],_next[y]=Lnum;  
             line[Lnum].x=y,line[Lnum].y=x,line[Lnum].c=0;               
      }  
      bool BFS(int s,int e)  
      {   
             queue<int> Q;  
             while (!Q.empty()) Q.pop();  
             memset(dis,0,sizeof(dis));  
             dis[s]=1,Q.push(s);  
             while (!Q.empty())  
             {  
                   int h,k;  
                   h=Q.front(),Q.pop();  
                   if (h==e) return dis[e];  
                   for (k=_next[h];k!=-1;k=line[k].next)  
                      if (line[k].c && !dis[line[k].y])  
                         dis[line[k].y]=dis[h]+1,Q.push(line[k].y);                   
             }   
             return false;  
      }  
      int dfs(int x,int flow,int e)    
      {        
             if (x==e) return flow;     
             int temp,cost=0;    
             for (int k=_next[x];k!=-1;k=line[k].next)    
             if (line[k].c && dis[line[k].y]==dis[x]+1)    
             {    
                    temp=dfs(line[k].y,min(flow-cost,line[k].c),e);     
                    if (temp)    
                    {    
                           line[k].c-=temp,line[k^1].c+=temp;    
                           cost+=temp;    
                           if (flow==cost) return cost;    
                    }else dis[line[k].y]=-1;    
             }    
             return cost;    
      }    
      int MaxFlow(int s,int e)  
      {  
             int MaxFlow=0;   
             while (BFS(s,e)) MaxFlow+=dfs(s,oo,e);   
             return MaxFlow;  
      }  
}T;   
int main()
{    
      int n,m,p,i,j,x,y,s,e;  
      while (~scanf("%d%d",&n,&m))
      {
               memset(sum,0,sizeof(sum));
               for (i=1;i<=n;i++)
               {
                      y=0;
                      for (j=0;j<m;j++)
                      {
                             scanf("%d",&x);
                             if (x) y+=1<<j;
                      }
                      sum[y]++;
               } 
               s=1500,e=s+1,T.initial(e);
               p=(1<<m)+1;
               for (i=0;i<(1<<m);i++)
               {
                      T.addline(s,i,sum[i]);
                      for (x=0;x<m;x++) 
                          if (i&(1<<x)) T.addline(i,x+p,oo); 
               }
               for (i=0;i<m;i++) 
                   scanf("%d",&x),T.addline(i+p,e,x); 
               if (T.MaxFlow(s,e)==n) printf("YES\n");
                                else  printf("NO\n");
      }
      return 0;
}