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[400005] ; int m,n,cnt,head[100005],level[100005],q[100005] ; 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=0 ;e[cnt].nxt=head[t] ;head[t]=cnt++ ; } bool build(int s,int t) { int front=0,rear=0 ; memset(level,-1,sizeof(level)) ; q[rear++]=s ; level[s]=1 ; while(front<rear) { int u=q[front++] ; for(int i=head[u] ;i!=-1 ;i=e[i].nxt) { int tt=e[i].t ; if(level[tt]==-1 && e[i].cap>0) { level[tt]=level[u]+1 ; 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=0,a ; for(int i=head[s] ;i!=-1 ;i=e[i].nxt) { int tt=e[i].t ; if(level[tt]==level[s]+1 && e[i].cap>0) { a=find(tt,t,min(e[i].cap,flow-ret)) ; e[i].cap-=a ; e[i^1].cap+=a ; ret+=a ; if(ret==flow) return ret ; } } if(!ret)level[s]=-1 ; return ret ; } int dinic(int s,int t) { int flow,ret=0 ; 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,-1,sizeof(head)) ; cnt=0 ; for(int i=0 ;i<n ;i++) { int res=0 ; for(int j=0 ;j<m ;j++) { int x ; scanf("%d",&x) ; if(x)res+=(1<<(m-j-1)) ; } mp[res]++ ; } int s=(1<<m) ; int S,T ; S=0 ;T=s+m+1 ; for(int i=0 ;i<s ;i++) { if(mp[i])add(S,i+1,mp[i]) ; int temp=i ; int st=s+m ; while(temp) { if(temp&1)add(i+1,st,mp[i]) ; temp>>=1 ; st-- ; } } for(int i=0 ;i<m ;i++) { int x ; scanf("%d",&x) ; add(s+i+1,T,x) ; } int ans=dinic(S,T) ; if(ans==n)puts("YES") ; else puts("NO") ; } return 0 ; }