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=;
 const int maxm=;
 int n,m,siz;
 struct nod{
     int y,next;
 }e[maxn];int head[maxm]={},tot=;
 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]==){tp[y]=;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=;i<=siz;++i)fa[i]=i;
     memset(tp,,sizeof(tp));memset(pre,,sizeof(pre));
     tl=tr=;q[]=s;tp[s]=;
     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]==)continue;//cout<<i<<endl;
             if(!tp[y]){
                 tp[y]=;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 ;
                 }
                 tp[d[y]]=;q[++tr]=d[y];
             }
             else{
                 int pa=lca(x,y);
                 mlink(x,y,pa);mlink(y,x,pa);
             }
         }
     }
     return ;
 }
 int main(){
     int T;scanf("%d",&T);
     while(T-->){
         memset(bel,,sizeof(bel));memset(d,,sizeof(d));
         memset(head,,sizeof(head));
         tot=tly=;
         int p;
         scanf("%d%d%d",&n,&m,&p);siz=n+m*;
         int x,y;
         for(int i=;i<=m;++i){
             init(i,i+m);init(i+m,i);
             init(i,i+m*);init(i+m*,i);
             init(i+m,i+m*);init(i+m*,i+m);
         }
         for(int i=;i<=p;++i){
             scanf("%d%d",&x,&y);
             init(m*+x,y);init(y,m*+x);
             init(m*+x,y+m);init(y+m,m*+x);
             init(m*+x,y+m*);init(y+m*,m*+x);
         }
         int ans=;
         for(int i=;i<=siz;++i)if(!d[i])ans+=doit(i);
         printf("%d\n",ans-n);
     }
     return ;
 }