hdu 5352MZL's City(15多校第五场1010)

本文介绍了一种解决二分图中一对多匹配问题的方法,通过使用匈牙利算法来寻找最大匹配。针对特定的问题场景,如点的选择可能重叠且每次选择独立,文章详细介绍了如何构造二分图并实现匹配过程。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

//15多校第五场
//hdu 5352
//对于每次选点,可能会有重叠
//每次选点是独立的,待选择点也是独立的
//所以可以把所有的选择情况看成一张二分图
//对每一个X点,最多匹配K个Y点
//匈牙利算法即可。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include <fstream>
#include<set>
using namespace std;
const int INF=0x7fffffff;
typedef long long LL;
const int maxm=40010,maxn=600,maxnn=100010;//最多500条双向边
int mp[maxn][maxn];

int parent[maxn];

int root(int x){
    if(parent[x]==-1)return x;
    return parent[x]=root(parent[x]);
}
void merge(int a,int b){
    int m=root(a),n=root(b);
    if(m!=n)
        parent[m]=n;
}
struct edge{
    int u,v,next;
}edge[maxm];
int cnt,head[maxn];
void addedge(int u,int v){
    edge[cnt].u=u;edge[cnt].v=v;edge[cnt].next=head[u];head[u]=cnt++;
}
void buildedge(int sta,int x,int n){
    for(int i=1;i<=n;i++){
        if(root(i)==root(x)){
            addedge(sta,i);
        }
    }
}
void del(int n){
    memset(parent,-1,sizeof(parent));
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            if(i==j)continue;
            if(mp[i][j]==1)merge(i,j);
        }
    }
}
int line[maxn];
int times[maxn];//记录X点匹配的数量
//普通匈牙利算法中的used数组是处理当前y点是否尝试匹配过的情况
//如果只是匹配一个,那么匹配哪个都一样,可以跳过。
//但是,如果要匹配多个,就可能存在这个点改变匹配情况,出现一条新的增广路的情况
//所以这里不能用used数组标记
int find(int u,int ny,int k){//一对多
    int ret=0;
    for(int i=head[u];i!=-1;i=edge[i].next,head[u]=i){//每条边只用一次,每次循环最多加一
        int v=edge[i].v;
        if(times[u]>=k)break;//匹配数大于等于N跳出
        if(line[v]==-1){//
            line[v]=u;
            times[u]++;
            ret+=1;
        }
        else {
            times[line[v]]--;//匹配数先减一,防止跳出
            if(find(line[v],ny,k)){
                    line[v]=u;
                    times[u]++;
                    ret+=1;
                  }
                 else times[line[v]]++;//无法增广时,匹配数还原。不用考虑可以增广的情况
                }
    }
    return ret;
}
int num[maxn];
int Hungry(int nx,int ny,int k){
    int all=0;
    memset(line, -1, sizeof(line));
    memset(times, 0, sizeof(times));
    for(int i=nx;i>=1;i--){//保证优先选后面的
        all+=find(i,ny,k);
    }
    memset(num, 0, sizeof(num));
    for(int i=1;i<=ny;i++){
        if(line[i]==-1)continue;
        num[line[i]]++;
    }
    return all;
}
int main(){
    int T;
    scanf("%d",&T);
    for(int i=0;i<T;i++){
        int n,m,k;
        int nx,ny;
       
        memset(parent,-1,sizeof(parent));
        memset(mp, 0, sizeof(mp));
        memset(head, -1, sizeof(head));
        cnt=0;
        scanf("%d%d%d",&n,&m,&k);
        ny=n;
        nx=0;
        for(int i=0;i<m;i++){
            int op,x,y;
            scanf("%d",&op);
            if(op==1){
                scanf("%d",&x);
                buildedge(++nx,x,n);
            }
            else if(op==2){
                scanf("%d%d",&x,&y);
                mp[x][y]=mp[y][x]=1;
                merge(x,y);
            }
            else {
                scanf("%d",&x);
                int u,v;
                for(int i=0;i<x;i++){
                    scanf("%d%d",&u,&v);
                    mp[u][v]=mp[v][u]=0;
                }
                del(n);
            }
        }
        printf("%d\n",Hungry(nx,ny,k));
        for(int i=1;i<=nx;i++){
            if(i!=1)printf(" ");
            printf("%d",num[i]);
        }
        printf("\n");
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值