HDU 3605

http://acm.hdu.edu.cn/showproblem.php?pid=3605

用最大流做的，G++超时，C++可以过，看别人写的叫二分图多重匹配，还不会这玩意一会学学

显然的最大流模型，n是100w，直接做最大流一定超时，但是注意到m只有10，所以可以对所有点进行状态压缩，状态相同的点可以进行合并，这样n被压缩到1024个点，水水过

#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
#include <map>
using namespace std ;

const int INF=0xfffffff ;
struct node
{
    int s,t,cap,nxt ;
}e[] ;
int m,n,cnt,head[],level[],q[] ;
void add(int s,int t,int cap)
{
    e[cnt].s=s ;e[cnt].t=t ;e[cnt].cap=cap ;e[cnt].nxt=head[s] ;head[s]=cnt++ ;
    e[cnt].s=t ;e[cnt].t=s ;e[cnt].cap= ;e[cnt].nxt=head[t] ;head[t]=cnt++ ;
}
bool build(int s,int t)
{
    int front=,rear= ;
    memset(level,-,sizeof(level)) ;
    q[rear++]=s ;
    level[s]= ;
    while(front<rear)
    {
        int u=q[front++] ;
        for(int i=head[u] ;i!=- ;i=e[i].nxt)
        {
            int tt=e[i].t ;
            if(level[tt]==- && e[i].cap>)
            {
                level[tt]=level[u]+ ;
                if(tt==t)return true ;
                q[rear++]=tt ;
            }
        }
    }
    return false ;
}
int find(int s,int t,int flow)
{
    if(s==t)return flow ;
    int ret=,a ;
    for(int i=head[s] ;i!=- ;i=e[i].nxt)
    {
        int tt=e[i].t ;
        if(level[tt]==level[s]+ && e[i].cap>)
        {
            a=find(tt,t,min(e[i].cap,flow-ret)) ;
            e[i].cap-=a ;
            e[i^].cap+=a ;
            ret+=a ;
            if(ret==flow)
                return ret ;
        }
    }
    if(!ret)level[s]=- ;
    return ret ;
}
int dinic(int s,int t)
{
    int flow,ret= ;
    while(build(s,t))
        while(flow=find(s,t,INF))
            ret+=flow ;
    return ret ;
} 

int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        map <int,int> mp ;
        memset(head,-,sizeof(head)) ;
        cnt= ;
        for(int i= ;i<n ;i++)
        {
            int res= ;
            for(int j= ;j<m ;j++)
            {
                int x ;
                scanf("%d",&x) ;
                if(x)res+=(<<(m-j-)) ;
            }
            mp[res]++ ;
        }
        int s=(<<m) ;
        int S,T ;
        S= ;T=s+m+ ;
        for(int i= ;i<s ;i++)
        {
            if(mp[i])add(S,i+,mp[i]) ;
            int temp=i ;
            int st=s+m ;
            while(temp)
            {
                if(temp&)add(i+,st,mp[i]) ;
                temp>>= ;
                st-- ;
            }
        }
        for(int i= ;i<m ;i++)
        {
            int x ;
            scanf("%d",&x) ;
            add(s+i+,T,x) ;
        }
        int ans=dinic(S,T) ;
        if(ans==n)puts("YES") ;
        else puts("NO") ;
    }
    return  ;
}