[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。

多组数据！！！

---

sb当前弧优化回我青春。

三个小时三个人的艰苦奋斗找到了一句话写错...

---

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
inline int read(){
    int res=0;char ch=getchar();
    while(!isdigit(ch))ch=getchar();
    while(isdigit(ch)){res=(res<<3)+(res<<1)+(ch^48);ch=getchar();}
    return res; 
}
int Ti, n;
int S, T;
struct edge{
    int nxt, to, w;
}ed[2505];
int head[2505], cnt = 1;
inline void add(int x, int y, int z)
{
    ed[++cnt] = (edge){head[x], y, z};
    head[x] = cnt;
}
inline void ADD(int x, int y, int z)
{
    add(x, y, z), add(y, x, 0);
}
bool back[55], is[55];
bool know[55][55];
int d[105], cur[105];

inline bool bfs()
{
    queue <int> q;
    while(!q.empty()) q.pop();
    memset(d, 0, sizeof d);
    d[S] = 1;
    q.push(S);
    while(!q.empty())
    {
        int x = q.front();
        q.pop();
        for (int i = head[x] ; i ; i = ed[i].nxt)
        {
            int to = ed[i].to;
            if (!d[to] and ed[i].w > 0)
            {
                d[to] = d[x] + 1;
                q.push(to);
                if (to == T) return 1;
            }            
        }
    }
    return 0;
}
inline void init()
{
    S = n + n + 1,T = n + n + 2;
    memset(head, 0, sizeof head);
    cnt = 1;
    memset(back, 0, sizeof back);
    memset(is, 0, sizeof is);
}

int dfs(int x, int flow)
{
    if (x == T) return flow;
    for (int i = head[x] ; i ; i = ed[i].nxt)
    {
        int to = ed[i].to;
        if (d[to] == d[x] + 1 and ed[i].w)
        {
            int di = dfs(to, min(flow, ed[i].w));
            if (di){
                ed[i].w -= di;
                ed[i^1].w += di;
                return di;
            }
        }
    }
    return 0;
}

int main()
{
    Ti = read();
    while(Ti--)
    {
        n = read();
        init();
        for (int i = 1 ; i <= n ; i ++) is[i] = read();
        for (int i = 1 ; i <= n ; i ++) {back[i] = read();if(!is[i]) back[i] = 0;}
        for (int i = 1 ; i <= n ; i ++)
            for (int j = 1 ; j <= n ; j ++)
                know[i][j] = read();
        for (int i = 1 ; i <= n ; i ++)
            for (int j = 1 ; j <= n ; j ++)
                if (i == j or know[i][j] == 1) ADD(i, j+n, 1);
        for (int i = 1 ; i <= n ; i ++) if(!back[i]) ADD(S, i, 1);
        for (int i = 1 ; i <= n ; i ++) if (is[i] == 1) ADD(i + n, T, 1);
        int ans = 0;
        while(bfs())
        {
            for (int i = 1 ; i <= n + n ; i ++) cur[i] = head[i];
            while(int di = dfs(S, 1e9)) ans += di; 
        }
        int tot = 0;
        for (int i = 1 ; i <= n ; i ++) if(!back[i]) tot++;
        if (tot != ans) puts("T_T");
        else puts("^_^");
    }
    return 0;
}