「BZOJ1433」[ZJOI2009] 假期的宿舍(二分图,网络流)

题目描述

学校放假了 · · · · · · 有些同学回家了,而有些同学则有以前的好朋友来探访,那么住宿就是一个问题。比如 A 和 B 都是学校的学生,A 要回家,而 C 来看B,C 与 A 不认识。我们假设每个人只能睡和自己直接认识的人的床。那么一个解决方案就是 B 睡 A 的床而 C 睡 B 的床。而实际情况可能非常复杂,有的人可能认识好多在校学生,在校学生之间也不一定都互相认识。我们已知一共有 n 个人,并且知道其中每个人是不是本校学生,也知道每个本校学生是否回家。问是否存在一个方案使得所有不回家的本校学生和来看他们的其他人都有地方住。

输入输出格式

输入格式:

 

第一行一个数 T 表示数据组数。接下来 T 组数据,每组数据第一行一个数n 表示涉及到的总人数。接下来一行 n 个数,第 i 个数表示第 i 个人是否是在校学生 (0 表示不是,1 表示是)。再接下来一行 n 个数,第 i 个数表示第 i 个人是否回家 (0 表示不回家,1 表示回家,注意如果第 i 个人不是在校学生,那么这个位置上的数是一个随机的数,你应该在读入以后忽略它)。接下来 n 行每行 n 个数,第 i 行第 j 个数表示 i 和 j 是否认识 (1 表示认识,0 表示不认识,第 i 行 i 个的值为 0,但是显然自己还是可以睡自己的床),认识的关系是相互的。

 

输出格式:

 

对于每组数据,如果存在一个方案则输出 “^_^”(不含引号) 否则输出“T_T”(不含引号)。(注意输出的都是半角字符,即三个符号的 ASCII 码分别为94,84,95)

 

输入输出样例

输入样例#1: 复制
1
3
1 1 0
0 1 0
0 1 1
1 0 0
1 0 0
输出样例#1: 复制
^_^

说明

对于 30% 的数据满足 1 ≤ n ≤ 12。

对于 100% 的数据满足 1 ≤ n ≤ 50,1 ≤ T ≤ 20。

多组数据!!!

题解

  比较裸的二分图匹配

  左边为人,右边为床,然后看最大匹配是否等于左边的人数

  如果一个人留校,那么将他与自己的床连边

  如果a与b认识且b是学生,那么把a与b的床连边

  然后求一个最大匹配

  我是用网络流dinic,也可以直接匈牙利(然而我不会匈牙利)

 1 //minamoto
 2 #include<iostream>
 3 #include<cstdio>
 4 #include<cstring>
 5 #include<queue>
 6 #define inf 0x3f3f3f3f
 7 using namespace std;
 8 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
 9 char buf[1<<21],*p1=buf,*p2=buf;
10 inline int read(){
11     #define num ch-'0'
12     char ch;bool flag=0;int res;
13     while(!isdigit(ch=getc()))
14     (ch=='-')&&(flag=true);
15     for(res=num;isdigit(ch=getc());res=res*10+num);
16     (flag)&&(res=-res);
17     #undef num
18     return res;
19 }
20 const int N=105,M=10005;
21 int ver[M],Next[M],head[N],edge[M],tot;
22 int dep[N],cur[N],n,m,s,t;
23 int is[N],cnt;
24 queue<int> q;
25 inline void add(int u,int v,int e){
26     ver[++tot]=v,Next[tot]=head[u],head[u]=tot,edge[tot]=e;
27     ver[++tot]=u,Next[tot]=head[v],head[v]=tot,edge[tot]=0;
28 }
29 bool bfs(){
30     while(!q.empty()) q.pop();
31     memset(dep,-1,sizeof(dep));
32     for(int i=s;i<=t;++i) cur[i]=head[i];
33     q.push(s),dep[s]=0;
34     while(!q.empty()){
35         int u=q.front();q.pop();
36         for(int i=head[u];i;i=Next[i]){
37             int v=ver[i];
38             if(dep[v]<0&&edge[i]){
39                 dep[v]=dep[u]+1,q.push(v);
40                 if(v==t) return true;
41             }
42         }
43     }
44     return false;
45 }
46 int dfs(int u,int limit){
47     if(!limit||u==t) return limit;
48     int flow=0,f;
49     for(int i=cur[u];i;i=Next[i]){
50         int v=ver[i];cur[u]=i;
51         if(dep[v]==dep[u]+1&&(f=dfs(v,min(limit,edge[i])))){
52             flow+=f,limit-=f;
53             edge[i]-=f,edge[i^1]+=f;
54             if(!limit) break;
55         }
56     }
57     if(!flow) dep[u]=-1;
58     return flow;
59 }
60 int dinic(){
61     int flow=0;
62     while(bfs()) flow+=dfs(s,inf);
63     return flow;
64 }
65 int main(){
66     int T=read();
67     while(T--){
68         n=read(),cnt=0,tot=1;
69         memset(head,0,sizeof(head));
70         s=0,t=2*n+1;
71         for(int i=1;i<=n;++i) is[i]=read();
72         for(int i=1;i<=n;++i){
73             int k=read();
74             if(!is[i]) add(s,i,1),++cnt;
75             else if(!k) add(s,i,1),add(i,i+n,1),++cnt;
76         }
77         for(int i=1;i<=n;++i)
78         for(int j=1;j<=n;++j){
79             int x=read();
80             if(x&&is[j]) add(i,n+j,1);
81         }
82         for(int i=1;i<=n;++i) add(i+n,t,1);
83         puts((dinic()==cnt)?"^_^":"T_T");
84     }
85     return 0;
86 }

 

posted @ 2018-08-21 16:35  bztMinamoto  阅读(93)  评论(0编辑  收藏
Live2D