bzoj4405 [wc2016]挑战NPC(一般图最大匹配)

我们把每一个筐子拆成三个,互相连边,得到一个三元环。x向y拆成的三个点连边。我们不难发现,每一个匹配都对应了一种分配方法。如果一个筐子至多放了1个球,则它的三元环中自己还有且仅有一对匹配。所以半空的筐子数=最大匹配-n
然后是个一般图,要用带花树来求最大匹配,模板见uoj79

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define N 610
#define ll long long
#define inf 0x3f3f3f3f
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,E,h[N],num=0,fa[N],col[N],ans=0,mate[N],tot;
int q[N],qh,qt,pre[N],tim=0,vis[N];
struct edge{
    int to,next;
}data[181000];
inline void add(int x,int y){
    data[++num].to=y;data[num].next=h[x];h[x]=num;
    data[++num].to=x;data[num].next=h[y];h[y]=num;
}
inline int find(int x){return x==fa[x]?x:fa[x]=find(fa[x]);}
inline int lca(int x,int y){
    x=find(x),y=find(y);++tim;
    for(;;swap(x,y))
        if(x){
            if(vis[x]==tim) return x;
            vis[x]=tim;
            if(mate[x]) x=find(pre[mate[x]]);
            else x=0;
        }
}
inline void blossom(int x,int y,int t){
    while(find(x)!=t){
        pre[x]=y;
        y=mate[x];
        if(col[y]==1) col[y]=0,q[++qt]=y;
        fa[x]=fa[y]=t;x=pre[y];
    }
}
inline bool match(int s){
    memset(col,-1,sizeof(col));memset(pre,0,sizeof(pre));
    for(int i=1;i<=tot;++i) fa[i]=i;qh=1;qt=0;
    col[s]=0;q[++qt]=s;
    while(qh<=qt){
        int x=q[qh++];
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;
            if(col[y]==-1){
                col[y]=1;pre[y]=x;
                if(!mate[y]){
                    for(int lst;x;y=lst,x=pre[y])
                        lst=mate[x],mate[x]=y,mate[y]=x;return 1;
                }else col[q[++qt]=mate[y]]=0;
            }else if(!col[y]&&find(x)!=find(y)){
                int t=lca(x,y);
                blossom(x,y,t);blossom(y,x,t);
            }
        }
    }return 0;
}
int main(){
//  freopen("npc.in","r",stdin);
    int tst=read();
    while(tst--){
        num=0;memset(h,0,sizeof(h));memset(mate,0,sizeof(mate));
        n=read();m=read();E=read();ans=0;tim=0;memset(vis,0,sizeof(vis));
        while(E--){
            int x=read(),y=read();
            add(x,y+n);add(x,y+n+m);add(x,y+n+2*m);
        }for(int i=1;i<=m;++i) add(i+n,i+n+m);tot=n+3*m;
        for(int i=1;i<=tot;++i)
            if(!mate[i]&&match(i)) ans++;
        printf("%d\n",ans-n);
        for(int i=1;i<=n;++i){
            int y=mate[i]-n;while(y>m) y-=m;
            printf("%d%c",y,i==n?'\n':' ');
        }
    }return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值