HDOJ 3605 Escape 解题报告 二分图匹配 dinic 网络流

【Problem Description】
2012 If this is the end of the world how to do? I do not know how. But now scientists have found that some stars, who can live, but some people do not fit to live some of the planet. Now scientists want your help, is to determine what all of people can live in these planets.
题目大意：
这道题就是给你n个人，m个星球，这n个人对m个星期有不同的意向。然而，每个星球都有自己的承受能力。我们想知道能否全部满足这n个人的意向。
【解题报告】
我们不难想出这个网络流怎么建：
源点S连接每个人，边权为1，每个人同星球相连，边权为1，每个星球通汇点T相连，边权即承受能力。最后再跑一边dinic。
然而人们发现，这样一来铁定会超时（点太多了）。所以，我们采取了一个缩点的策略，将去同一个星球的人缩成了一个点，这样一来，每个人同星球的边权自然就不是1，而是人数。
下面我们看看代码：

#include<cstdio>  
#include<cstring>  
#include<string>  
#include<algorithm>  
#include<vector>  
#include<queue>  
#include<cmath>  
using namespace std;  
#define Del(a,b) memset(a,b,sizeof(a))  
const int N=1210;  
const int inf=0x3f3f3f3f;  
const double esp=1e-9;  
int n,m;  
struct Node  
{  
    int from,to,cap,flow;  
};  
vector<int> v[N];  
int sum[N];  
vector<Node> e;  
int vis[N];  //构建层次图  
int cur[N];  
int MIN(int x,int y)  
{  
    return x>y?y:x;  
}  
void add_Node(int from,int to,int cap)  
{  
    Node tmp1,tmp2;  
    tmp1.from=from,tmp1.to=to,tmp1.cap=cap,tmp1.flow=0;  
    e.push_back(tmp1);  
    tmp2.from=to,tmp2.to=from,tmp2.cap=0,tmp2.flow=0;  
    e.push_back(tmp2);  
    int tmp=e.size();  
    v[from].push_back(tmp-2);  
    v[to].push_back(tmp-1);  
}  
bool bfs(int s,int t)  
{  
    Del(vis,-1);  
    queue<int> q;  
    q.push(s);  
    vis[s] = 0;  
    while(!q.empty())  
    {  
        int x=q.front();  
        q.pop();  
        int i;  
        for(i=0;i<v[x].size();i++)  
        {  
            Node tmp = e[v[x][i]];  
            if(vis[tmp.to]<0 && tmp.cap>tmp.flow)  //第二个条件保证  
            {  
                vis[tmp.to]=vis[x]+1;  
                q.push(tmp.to);  
            }  
        }  
    }  
    if(vis[t]>0)  
        return true;  
    return false;  
}  
int dfs(int o,int f,int t)  
{  
    if(o==t || f==0)  //优化  
        return f;  
    int a = 0,ans=0;  
    int &i = cur[o];  
    for(;i<v[o].size();i++) //注意前面 ’&‘,很重要的优化  
    {  
        Node &tmp = e[v[o][i]];  
        if(vis[tmp.to]==(vis[o]+1) && (a = dfs(tmp.to,MIN(f,tmp.cap-tmp.flow),t))>0)  
        {  
            tmp.flow+=a;  
            e[v[o][i]^1].flow-=a; //存图方式  
            ans+=a;  
            f-=a;  
            if(f==0)  //注意优化  
                break;  
        }  
    }  
    return ans;  //优化  
}  
int dinci(int s,int t)  
{  
    int ans=0;  
    while(bfs(s,t))  
    {  
        Del(cur,0);  
        int tm=dfs(s,inf,t);  
        ans+=tm;  
    }  
    return ans;  
}  
int main()  
{  
    int n,m;  
    while(~scanf("%d%d",&n,&m))  
    {  
        int i,j;  
        memset(sum,0,sizeof(sum));  
        int s=0,x,t=m+(1<<m)+3,nn=(1<<m)+2;  
        for(i=1;i<=n;i++)  
        {  
            int tmp=0;  
            for(j=0;j<m;j++)  
            {  
                scanf("%d",&x);  
                if(x)  
                    tmp+=(1<<j);  
            }  
            sum[tmp]++;  
        }  
        for(i=0;i<(1<<m);i++)  
        {  
            if(sum[i])  
            {  
                add_Node(s,i+1,sum[i]);  
                for(j=0;j<m;j++)  
                {  
                    if(i&(1<<j))  
                        add_Node(i+1,nn+j,sum[i]);  
                }  
            }  
        }  
        for(i=0;i<m;i++)  
        {  
            scanf("%d",&x);  
            add_Node(nn+i,t,x);  
        }  
        int ans=dinci(s,t);  
        //printf("%d\n",ans);  
        if(ans==n)  
            printf("YES\n");  
        else  
            printf("NO\n");  
        for(i=0;i<=t;i++)  
            v[i].clear();  
        e.clear();  
    }  
    return 0;  
} 

以上