[网络流]假期的宿舍

题目描述

学校放假了
有些同学回家了，而有些同学则有以前的好朋友来探访，那么住宿就是一个问题。
比如 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。

分析

题意贼模糊
把第i个人的床编号为i+n
那么明显在校学生就连床，需要床的（不回家的）连源点，床全部连汇点
i认识j，就连j的床
然后跑网络流

#include <iostream>
#include <cstdio>
#include <queue>
#include <memory.h>
#define rep(i,a,b) for (i=a;i<=b;i++)
using namespace std;
int T,n;
int s,t;
int cnt,scnt;
int list[105],next[10005],u[10005],v[10005],w[10005];
int cur[105],d[105];
int ans;
void add(int u1,int v1,int w1)
{
    u[++cnt]=u1;v[cnt]=v1;next[cnt]=list[u1];list[u1]=cnt;w[cnt]=w1;
    u[++cnt]=v1;v[cnt]=u1;next[cnt]=list[v1];list[v1]=cnt;w[cnt]=0;
}
void init()
{
    int i,j,p;
    cnt=1;scnt=0;
    bool b[105],c[105];
    memset(list,0,sizeof(list));
    memset(u,0,sizeof(u));
    memset(v,0,sizeof(v));
    memset(w,0,sizeof(w));
    memset(next,0,sizeof(next));
    memset(cur,0,sizeof(cur));
    scanf("%d",&n);
    rep(i,1,n)
    {
        scanf("%d",&b[i]);
        if (b[i])
        add(i+n,2*n+1,1);
    }
    rep(i,1,n)
    {
        scanf("%d",&c[i]);
        if (!b[i]) c[i]=0;
        if (!c[i])
        {
            add(0,i,1);
            scnt++;
        }
    }
    rep(i,1,n)
    if (!c[i])
    rep(j,1,n)
    {
        scanf("%d",&p);
        if (p||i==j)
        add(i,j+n,1);
    }
    else
    rep(j,1,n)
    scanf("%d",&p);
    s=0;t=2*n+1;
}
bool bfs()
{
    int i,x;
    queue<int> q;
    while (!q.empty()) q.pop();
    memset(d,0,sizeof(d));
    d[s]=1;q.push(s);
    while (!q.empty())
    {
        x=q.front();q.pop();
        for (i=list[x];i;i=next[i])
        if (w[i]&&!d[v[i]])
        {
            d[v[i]]=d[x]+1;
            if (v[i]==t) return 1;
            q.push(v[i]);
        }
    }
    return 0;
}
int dfs(int p,int mf)
{
    int ret=0,f;
    if (p==t||!mf) return mf;
    for (int &i=cur[p];i;i=next[i])
    if (w[i]&&d[p]+1==d[v[i]])
    {
        f=dfs(v[i],min(w[i],mf-ret));
        w[i]-=f;
        w[i^1]+=f;
        ret+=f;
        if (mf==ret) break;
    }
    return ret;
}
void dinic()
{
    int i;
    ans=0;
    while (bfs())
    {
        rep(i,s,t) cur[i]=list[i];
        ans+=dfs(s,2147483647);
    }
}
int main()
{
    scanf("%d",&T);
    for (;T;T--)
    {
        init();
        dinic();
        if (ans>=scnt) printf("^_^\n");
        else printf("T_T\n");
    }
}