BZOJ 4405 [wc2016]挑战NPC 带花树 一般图最大匹配

https://www.lydsy.com/JudgeOnline/problem.php?id=4405

这道题大概就是考场上想不出来，想出来也调不出来的题。

把每个桶拆成三个互相有边的点，每个球向它连接的桶的三个点分别连边。

0球1桶，匹配数为1；1球1桶，匹配数为2；2球一桶，匹配数为2；3球一桶，匹配数为3；

发现每种半桶的情况下匹配数都比该桶中放的球数大1，那么ans=最大匹配数-球数。

带花树找lca的时候，记得是用每个点的总父亲找。带花树的细节真是恶心人。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
#define LL long long
const int maxn=730010;
const int maxm=610;
int n,m,siz;
struct nod{
    int y,next;
}e[maxn];int head[maxm]={},tot=0;
int q[maxn]={},tl,tr;
int bel[maxm]={},tly,d[maxm]={},pre[maxm]={},tp[maxm]={},fa[maxm]={};
inline void init(int x,int y){e[++tot].y=y;e[tot].next=head[x];head[x]=tot;}
inline int getfa(int x){
    return fa[x]==x?x:fa[x]=getfa(fa[x]);
}
inline int lca(int x,int y){
    ++tly;
    for(;;){
        if(x){
            x=getfa(x);
            if(bel[x]==tly)return x;
            else{bel[x]=tly;x=pre[d[x]];}
        }swap(x,y);        
    }
}
inline void mlink(int x,int y,int pa){
    while(getfa(x)!=pa){
        pre[x]=y;y=d[x];
        if(tp[y]==2){tp[y]=1;q[++tr]=y;}
        if(getfa(x)==x)fa[x]=pa;
        if(getfa(y)==y)fa[y]=pa;
        x=pre[y];
    }
}
int doit(int s){
    for(int i=1;i<=siz;++i)fa[i]=i;
    memset(tp,0,sizeof(tp));memset(pre,0,sizeof(pre));
    tl=tr=1;q[1]=s;tp[s]=1;
    while(tl<=tr){
        int x=q[tl];++tl;
        for(int i=head[x];i;i=e[i].next){
            int y=e[i].y;
            if(getfa(y)==getfa(x)||tp[y]==2)continue;//cout<<i<<endl;
            if(!tp[y]){
                tp[y]=2;pre[y]=x;//cout<<y<<endl;
                if(!d[y]){
                    for(int now=y,j,las;now;now=las){
                        j=pre[now];las=d[pre[now]];
                        d[j]=now;d[now]=j;
                    }//cout<<i<<1<<endl;
                    return 1;
                }
                tp[d[y]]=1;q[++tr]=d[y];
            }
            else{
                int pa=lca(x,y);
                mlink(x,y,pa);mlink(y,x,pa);
            }
        }
    }
    return 0;
}
int main(){
    int T;scanf("%d",&T);
    while(T-->0){
        memset(bel,0,sizeof(bel));memset(d,0,sizeof(d));
        memset(head,0,sizeof(head));
        tot=tly=0;
        int p;
        scanf("%d%d%d",&n,&m,&p);siz=n+m*3;
        int x,y;
        for(int i=1;i<=m;++i){
            init(i,i+m);init(i+m,i);
            init(i,i+m*2);init(i+m*2,i);
            init(i+m,i+m*2);init(i+m*2,i+m);
        }
        for(int i=1;i<=p;++i){
            scanf("%d%d",&x,&y);
            init(m*3+x,y);init(y,m*3+x);
            init(m*3+x,y+m);init(y+m,m*3+x);
            init(m*3+x,y+m*2);init(y+m*2,m*3+x);
        }
        int ans=0;
        for(int i=1;i<=siz;++i)if(!d[i])ans+=doit(i);
        printf("%d
",ans-n);
    }
    return 0;
}