[bzoj1433][ZJOI2009]假期的宿舍 二分图匹配 匈牙利算法/网络流

本文探讨了一种宿舍分配问题,即如何合理安排在校学生及其来访的朋友的住宿,确保每个人都能够有床位。通过构建特定的数据结构和使用匈牙利算法或网络流算法,解决了这一问题并提供了有效的解决方案。

Description

学校放假了······有些同学回家了,而有些同学则有以前的好朋友来探访,那么住宿就是一个问题。比如A

和B都是学校的学生,A要回家,而C来看B,C与A不认识。我们假设每个人只能睡和自己直接认识的人的床。那么一

个解决方案就是B睡A的床而C睡B的床。而实际情况可能非常复杂,有的人可能认识好多在校学生,在校学生之间也

不一定都互相认识。我们已知一共有n个人,并且知道其中每个人是不是本校学生,也知道每个本校学生是否回家

。问是否存在一个方案使得所有不回家的本校学生和来看他们的其他人都有地方住。

Input

第一行一个数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,但是显然自己还是可以睡自己的床),

认识的关系是相互的。

1 ≤ n ≤ 50,1 ≤ T ≤ 20

Output

对于每组数据,如果存在一个方案则输出“^_^”(不含引号)否则输出“T_T”(不含引号)。

(注意输出的都是半角字符,即三个符号的ASCII码分别为94,84,95)

Sample Input

1
3
1 1 0
0 1 0
0 1 1
1 0 0
1 0 0

Sample Output

ˆ ˆ

HINT

 

Source

 

思维难度:NOIP+

代码难度:NOIP

算法:匈牙利算法/网络流

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int Maxn=105;
int n,h[Maxn],cnt,in[Maxn],must[Maxn],vis[Maxn],now[Maxn],fir[Maxn],mat[Maxn];
struct node{
    int v,next;
}e[Maxn*Maxn*2];
void add(int u,int v){
    cnt++;
    e[cnt].v=v;
    e[cnt].next=h[u];
    h[u]=cnt;
}
int next(int x){
    if(fir[x]==0){fir[x]=1;return h[x];}
    return e[now[x]].next;
}
bool dfs(int x){
    for(int i=next(x);i;i=e[i].next){
        now[x]=i;
        int v=e[i].v;
        if(!vis[v]){
            vis[v]=x;
            return 1;
        }
        else{
            if(dfs(vis[v])){
                vis[v]=x;
                return 1;
            }
        }
    }
    return 0;
}
void work(){
    int u;
    memset(h,0,sizeof(h));cnt=0;
    memset(vis,0,sizeof(vis));memset(fir,0,sizeof(fir));memset(now,0,sizeof(now));
    memset(in,0,sizeof(in));memset(must,0,sizeof(must));memset(e,0,sizeof(e));
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d",&in[i]);
        if(in[i])add(i,i+n);
    }
    for(int i=1;i<=n;i++){
        scanf("%d",&u);
        if(in[i]){
            if(u==0){
                must[i]=1;
            }
        }
        else{
            must[i]=1;
        }
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            scanf("%d",&u);
            if(u==1&&in[j])add(i,j+n);
            if(u==1&&in[i])add(j,i+n);
        }
    }
    for(int i=n;i>=1;i--){
        if(must[i]){
            if(!dfs(i)){
                printf("T_T\n");
                return;
            }
        }
    }
    printf("^_^\n");
}
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        work();
    } 
    return 0;
}

 

人工智能关于最大二分图的程序代码: #include #include main() { bool map[100][300]; int i,i1,i2,num,num1,que[300],cou,stu,match1[100],match2[300],pque,p1,now,prev[300],n; scanf("%d",&n); for(i=0;i<n;i++) { scanf("%d%d",&cou;,&stu;); memset(map,0,sizeof(map)); for(i1=0;i1<cou;i1++) { scanf("%d",&num;); for(i2=0;i2<num;i2++) { scanf("%d",&num1;); map[i1][num1-1]=true; } } num=0; memset(match1,int(-1),sizeof(match1)); memset(match2,int(-1),sizeof(match2)); for(i1=0;i1<cou;i1++) { p1=0; pque=0; for(i2=0;i2<stu;i2++) { if(map[i1][i2]) { prev[i2]=-1; que[pque++]=i2; } else prev[i2]=-2; } while(p1<pque) { now=que[p1]; if(match2[now]==-1) break; p1++; for(i2=0;i2=0) { match1[match2[prev[now]]]=now; match2[now]=match2[prev[now]]; now=prev[now]; } match2[now]=i1; match1[i1]=now; num++; } if(num==cou) printf("YES\n"); else printf("NO\n"); } } dfs实现过程: #include #include #define MAX 100 bool map[MAX][MAX],searched[MAX]; int prev[MAX],m,n; bool dfs(int data) { int i,temp; for(i=0;i<m;i++) { if(map[data][i]&&!searched[i]) { searched[i]=true; temp=prev[i]; prev[i]=data; if(temp==-1||dfs(temp)) return true; prev[i]=temp; } } return false; } main() { int num,i,k,temp1,temp2,job; while(scanf("%d",&n)!=EOF&&n!=0) { scanf("%d%d",&m,&k); memset(map,0,sizeof(map)); memset(prev,int(-1),sizeof(prev)); memset(searched,0,sizeof(searched)); for(i=0;i<k;i++) { scanf("%d%d%d",&job;,&temp1;,&temp2;); if(temp1!=0&&temp2;!=0) map[temp1][temp2]=true; } num=0; for(i=0;i<n;i++) { memset(searched,0,sizeof(searched)); dfs(i); } for(i=0;i<m;i++) { if(prev[i]!=-1) num++; } printf("%d\n",num); } }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值