hdu3605（最大流+状态压缩）

传送门：Escape

题意：给出每个人适合住的星球信息和该星球能住多少人 ，第一行给出n m 代表有 n 个人 m 个星球，然后接下来n行每行m个数字 1代表适合第 i 个星球 0 代表不适合第 i 个星球，最后一行m个数表示第 i 个星球最多可以住多少个人，问是不是所有人都可以住到星球上。

分析：很裸的最大流，但刚开始直接去建图，100W级别的边MLE了，就算没MLE也会TLE的，因此需要状态压缩，m<=10，最多1024种状态，边数降成1W级别，最大流妥妥的A了。

#pragma comment(linker,"/STACK:1024000000,1024000000")
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <limits.h>
#include <iostream>
#include <algorithm>
#include <queue>
#include <cstdlib>
#include <stack>
#include <vector>
#include <set>
#include <map>
#define LL long long
#define mod 100000000
#define inf 0x3f3f3f3f
#define eps 1e-6
#define N 10010
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define PII pair<int,int>
using namespace std;
inline int read()
{
    char ch=getchar();int x=,f=;
    while(ch>''||ch<''){if(ch=='-')f=-;ch=getchar();}
    while(ch<=''&&ch>=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
int n,m,sum,vs,vt,tot,NV;
int pre[N],cur[N],h[N],q[N];
struct edge
{
    int u,v,w,next;
    edge() {}
    edge(int u,int v,int w,int next):u(u),v(v),w(w),next(next) {}
} e[N*];
void addedge(int u,int v,int w)
{
    e[tot]=edge(u,v,w,pre[u]);
    pre[u]=tot++;
    e[tot]=edge(v,u,,pre[v]);
    pre[v]=tot++;
}
void init()
{
    memset(pre,-,sizeof(pre));
    tot=;
}

/*******************dinic************************/
int bfs()
{
    int head=,tail=;
    memset(h,-,sizeof(h));
    q[]=vs;h[vs]=;
    while(head!=tail)
    {
        int u=q[head++];
        for(int i=pre[u];~i;i=e[i].next)
        {
            int v=e[i].v,w=e[i].w;
            if(w&&h[v]==-)
            {
                h[v]=h[u]+;
                q[tail++]=v;
            }
        }
    }
    return h[vt]!=-;
}

int dfs(int u,int flow)
{
    if(u==vt)return flow;
    int used=;
    for(int i=cur[u];~i;i=e[i].next)
    {
        int v=e[i].v,w=e[i].w;
        if(h[v]==h[u]+)
        {
            w=dfs(v,min(flow-used,w));
            e[i].w-=w;e[i^].w+=w;
            if(e[i].w)cur[u]=i;
            used+=w;
            if(used==flow)return flow;
        }
    }
    if(!used)h[u]=-;
    return used;
}
int dinic()
{
    int res=;
    while(bfs())
    {
        for(int i=vs;i<=vt;i++)cur[i]=pre[i];
        res+=dfs(vs,inf);
    }
    return res;
}
/********************dinic***********************/
int x,state,a[N];
void build()
{
    vs=;vt=(<<m)+m+;
    memset(a,,sizeof(a));
    for(int i=;i<=n;i++)
    {
        state=;
        for(int j=;j<m;j++)
        {
            x=read();
            if(x)state|=(<<j);
        }
        a[state]++;
    }
    for(int i=;i<(<<m);i++)
    {
        if(a[i])
        {
            addedge(vs,i+,a[i]);
            for(int j=;j<m;j++)
            {
                if(i&(<<j))
                    addedge(i+,j++(<<m),a[i]);
            }
        }
    }
    for(int i=;i<m;i++)
    {
        x=read();
        addedge((<<m)+i+,vt,x);
    }
    if(dinic()==n)puts("YES");
    else puts("NO");
}
int main()
{
    while(scanf("%d%d",&n,&m)>)
    {
        init();
        build();
    }
}