【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
3
1 1 0
0 1 0
0 1 1
1 0 0
1 0 0

样例输出

^_^

---

题解

建模：

把床位与源点连边，把需要床位的人与汇点连边，把需要床位的人与可以住的床连边，跑最大流。

其中人与床编号相同，只需要把其中一个变成 i + n 就行了。

#include<bits/stdc++.h>
using namespace std;
#define ll long long

const int maxn=50*2+50;
const int maxm=10000*2+50;

int T,n,s=0,t=101,tot,d[maxn],x,maxflow;
int fir[maxn],nex[maxm],to[maxm],wi[maxm],ecnt;
bool in[maxn],go[maxn],know[maxn][maxn];

void add(int u,int v,int w){
    nex[ecnt]=fir[u];fir[u]=ecnt;to[ecnt]=v;wi[ecnt++]=w;
}

queue<int> q;

int bfs(){
    memset(d,0,sizeof(d));
    d[s]=1;
    q.push(s);
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int e=fir[u];~e;e=nex[e]){
            int v=to[e];
            if(wi[e]&&!d[v]){
                d[v]=d[u]+1;
                q.push(v);
            }
        }
    }
    return d[t];
}

int dfs(int u,int flow){
    if(u==t) return flow;
    if(d[u]>=d[t]) return 0;
    for(int e=fir[u];~e;e=nex[e]){
        int v=to[e];
        if(wi[e]&&d[v]==d[u]+1&&(x=dfs(v,min(flow,wi[e])))){
            wi[e]-=x;
            wi[e^1]+=x;
            return x;
        }
    }
    return 0;
}

template<typename T>inline void read(T &aa){
    ll ff=1;char cc=getchar();aa=0;
    while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
    if(cc=='-') ff=-1,cc=getchar();
    while(cc<='9'&&cc>='0') aa=aa*10+cc-48,cc=getchar();
    aa*=ff;
}

int main(){
    read(T);
    while(T--){
        read(n);
        memset(fir,-1,sizeof(fir));
        memset(to,0,sizeof(to));
        memset(nex,0,sizeof(nex));
        memset(wi,0,sizeof(wi));ecnt=0;
        tot=0;maxflow=0;x=0;
        for(int i=1;i<=n;i++){
            read(in[i]);
            if(in[i]) add(s,i,1),add(i,s,0);
        }
        for(int i=1;i<=n;i++){
            read(go[i]);
            if((in[i]&&!go[i])||!in[i]) add(i+n,t,1),add(t,i+n,0),tot++;
        }
        for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++){
            read(know[i][j]);
            if(i==j||know[i][j]){
                add(i,j+n,1);add(j+n,i,0);
            }
        }
        while(bfs()) while(dfs(s,2e9)) maxflow+=x;
        if(maxflow>=tot) cout<<"^_^"<<endl;
        else cout<<"T_T"<<endl;
    }
    return 0;
}