HDU-3605 Escape

题目大意：

在地球上有n个人，但是地球快要game over了，所以这n个人要去别的星球生存，有m个适宜居住的星球，然后每个人对这m个星球都有一个适宜度，1代表适宜，0代表不适宜，最后给出m个星球每个能容纳多少人，问你这m个星球能否让这n个人搬去住

解题思路：

这是一道网络流的题，据说可以用二分图匹配来做，（萌新还没搞二分图

第一个思路：每个人与每个星球建图。这样TLE，MLE拿好不谢

第二个思路：先计算每个星球有多少人可以去，再与每个星球能容纳多少人比较。这样标准姿势的wa

第三个思路：有最多十个星球，这样把所有能去相同星球的人归为一类，比如a是0 0 0 0 1 0 0 0 1 1，b也是0 0 0 0 1 0 0 0 1 1，这样a和b是同一类，这样可以把10w个人分为最多1024类，然后这1024类分别与各个能够去的星球建图，设置超级源点和超级汇点，然后裸网络流，搞，ac

代码：

#include <map>
#include <cstdio>
#include <vector>
#include <cstring>
using namespace std;
#define Min(a, b) ((a) < (b) ? (a) : (b))

typedef struct node{
    int to, cap, rev;
    node(int a = 0, int b = 0, int c = 0){
        to = a; cap = b; rev = c;
    }
}Edge;

const int maxn = 15;
const int INF = 0x3f3f3f3f;

int a[maxn];
bool vis[2000];
vector<Edge> vec[2000];

void add(int from, int to, int cap){
    vec[from].push_back(Edge(to, cap, vec[to].size()));
    vec[to].push_back(Edge(from, 0, vec[from].size() - 1));
}
int dfs(int u, int v, int f){
    if(u == v) return f;
    vis[u] = 1;
    int len = vec[u].size();
    for(int i = 0; i < len; ++i){
        Edge &e = vec[u][i];
        if(!vis[e.to] && e.cap > 0){
            int d = dfs(e.to, v, Min(f, e.cap));
            if(d > 0){
                e.cap -= d;
                vec[e.to][e.rev].cap += d;
                return d;
            }
        }
    }
    return 0;
}
bool solve(int s, int t, int n){
    int flow = 0;
    while(true){
        memset(vis, 0, sizeof(vis));
        int d = dfs(s, t, INF);
        if(d == 0){
            if(flow == n) return 1;
            else return 0;
        }
        flow += d;
    }
}
int main(){
    int n, m, st, tt, val;
    while(~scanf("%d%d", &n, &m)){
        memset(a, 0, sizeof(a));
        map<int, int> mp; mp.clear();
        for(int i = 0; i < 2000; ++i) vec[i].clear();
        for(int i = 0; i < n; ++i){
            st = 0;
            for(int j = 0; j < m; ++j){
                st <<= 1;
                scanf("%d", &val);
                st |= val;
            }
            mp[st] += 1;
        }
        for(int i = 0; i < m; ++i){
            scanf("%d", &a[i]);
        }

        int status = (1 << m);
        for(int i = 1; i < (1 << m); ++i){
            if(mp[i] == 0) continue;
            st = i; val = mp[i]; tt = 0;
            add(0, i, val);
            while(st){
                if(st & 1) add(i, status + tt, val);
                st >>= 1;
                ++tt;
            }
        }

        for(int i = 0; i < m; ++i){
            add(status + i, status + m, a[i]);
        }

        if(solve(0, status + m, n)) puts("YES");
        else puts("NO");
    }
    return 0;
}