洛谷 P2055 [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
3
1 1 0
0 1 0
0 1 1
1 0 0
1 0 0
^_^
说明
对于 30% 的数据满足 1 ≤ n ≤ 12。
对于 100% 的数据满足 1 ≤ n ≤ 50,1 ≤ T ≤ 20。
多组数据!!!
Solution:
这道题一看就类似于二分图匹配的问题,既可以用匈牙利算法来做也可以用最大流做。
1、用匈牙利算法。思路就是把每个是学校里的学生和他们自己的床位相连且边权为1(实际实现时可以都连,只需要在增广时判断就行了),再统计需要在学校住的同学的人数(因为只有这些人需要床位),然后把相互认识的人和床位相连且边权为1(实际实现时这里也可以不必判断是否是学校的学生,只需要增广时特判就行了)。。。关键的增广时,记得判断:必须是留在学校的学生或者外来的学生,才能以他为起点进行增广。最后若增广累加的值等于一开始统计的需要在学校住的人数,则可以,否则不行。
代码:
1 #include<bits/stdc++.h> 2 #define il inline 3 using namespace std; 4 const int N=500,inf=2333333; 5 int T,n,l[N],k[N],cnt,pre[N],ans,tot,h[N]; 6 bool vis[N]; 7 struct edge{ 8 int to,net; 9 }e[N*2]; 10 il void add(int u,int v) 11 { 12 e[++cnt].to=v,e[cnt].net=h[u],h[u]=cnt; 13 } 14 il bool find(int u) 15 { 16 for(int i=h[u];i;i=e[i].net){ 17 int v=e[i].to; 18 if(!vis[v]){ 19 vis[v]=1; 20 if(!pre[v]||find(pre[v])){pre[v]=u;return 1;} 21 } 22 } 23 return 0; 24 } 25 int main() 26 { 27 scanf("%d",&T); 28 while(T--) 29 { 30 scanf("%d",&n); 31 ans=tot=0;cnt=0; 32 memset(h,0,sizeof(h)); 33 memset(pre,0,sizeof(pre)); 34 for(int i=1;i<=n;i++)scanf("%d",&l[i]); 35 for(int i=1;i<=n;i++){ 36 scanf("%d",&k[i]); 37 if(l[i]&&!k[i])add(i,i+n); 38 } 39 for(int i=1;i<=n;i++) 40 if(!l[i]||(l[i]&&!k[i]))tot++; 41 int x; 42 for(int i=1;i<=n;i++) 43 for(int j=1;j<=n;j++){ 44 scanf("%d",&x); 45 if(x&&l[j])add(i,j+n); 46 } 47 for(int i=1;i<=n;i++) 48 if(!l[i]||(l[i]&&!k[i])){ 49 memset(vis,0,sizeof(vis)); 50 if(find(i))ans++;} 51 // printf("%d %d\n",ans,tot); 52 if(ans==tot)printf("^_^\n"); 53 else printf("T_T\n"); 54 } 55 return 0; 56 }
2、用最大流做。思路与匈牙利算法的类似,但是附加源点S(连向需要住校的人并赋权值为1)和汇点T(由每个床位连向T并赋权值为1),然后与匈牙利算法相同,统计需要住校的人数,加边时注意必须是需要住校的且认识的才能与床位相连(该回家的肯定是回家好,这样不占床位),最后增广累加答案,若等于需要住校的人数则可以,否则不行。
代码:
1 #include<bits/stdc++.h> 2 #define il inline 3 using namespace std; 4 const int N=1005,inf=2333333; 5 int T,n,s,t=250,l[N],k[N],cnt=1,dis[N],h[N],ans,tot; 6 queue<int>q; 7 struct edge{ 8 int to,net,v; 9 }e[N*2]; 10 il void add(int u,int v,int w) 11 { 12 e[++cnt].to=v,e[cnt].v=w,e[cnt].net=h[u],h[u]=cnt; 13 e[++cnt].to=u,e[cnt].v=0,e[cnt].net=h[v],h[v]=cnt; 14 } 15 il bool bfs() 16 { 17 memset(dis,-1,sizeof(dis)); 18 q.push(s),dis[s]=0; 19 while(!q.empty()) 20 { 21 int u=q.front();q.pop(); 22 for(int i=h[u];i;i=e[i].net) 23 if(dis[e[i].to]==-1&&e[i].v>0)dis[e[i].to]=dis[u]+1,q.push(e[i].to); 24 } 25 return dis[t]!=-1; 26 } 27 il int dfs(int u,int op) 28 { 29 if(u==t)return op; 30 int flow=0,used=0; 31 for(int i=h[u];i;i=e[i].net) 32 { 33 int v=e[i].to; 34 if(dis[v]==dis[u]+1&&e[i].v>0) 35 { 36 used=dfs(v,min(op,e[i].v)); 37 if(!used)continue; 38 flow+=used,op-=used; 39 e[i].v-=used,e[i^1].v+=used; 40 if(!op)break; 41 } 42 } 43 if(!flow)dis[u]==-1; 44 return flow; 45 } 46 int main() 47 { 48 scanf("%d",&T); 49 while(T--) 50 { 51 scanf("%d",&n); 52 for(int i=0;i<=N;i++)h[i]=0,l[i]=0,k[i]=0; 53 //memset(e,0,sizeof(e)); 54 cnt=1;ans=tot=0; 55 int x; 56 for(int i=1;i<=n;i++)scanf("%d",&l[i]); 57 for(int i=1;i<=n;i++){ 58 scanf("%d",&x); 59 if(l[i]==1)k[i]=x; 60 } 61 for(int i=1;i<=n;i++){ 62 if(!k[i])add(s,i,1),tot++; 63 if(l[i])add(i+n,t,1); 64 } 65 for(int i=1;i<=n;i++) 66 for(int j=1;j<=n;j++){ 67 scanf("%d",&x); 68 if(x==1||i==j)add(i,j+n,1); 69 } 70 while(bfs())ans+=dfs(s,inf); 71 if(ans==tot)printf("^_^\n"); 72 else printf("T_T\n"); 73 } 74 return 0; 75 }