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的时候，记得是用每个点的总父亲找。带花树的细节真是恶心人。

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<cmath>
 6 #include<queue>
 7 using namespace std;
 8 #define LL long long
 9 const int maxn=730010;
10 const int maxm=610;
11 int n,m,siz;
12 struct nod{
13     int y,next;
14 }e[maxn];int head[maxm]={},tot=0;
15 int q[maxn]={},tl,tr;
16 int bel[maxm]={},tly,d[maxm]={},pre[maxm]={},tp[maxm]={},fa[maxm]={};
17 inline void init(int x,int y){e[++tot].y=y;e[tot].next=head[x];head[x]=tot;}
18 inline int getfa(int x){
19     return fa[x]==x?x:fa[x]=getfa(fa[x]);
20 }
21 inline int lca(int x,int y){
22     ++tly;
23     for(;;){
24         if(x){
25             x=getfa(x);
26             if(bel[x]==tly)return x;
27             else{bel[x]=tly;x=pre[d[x]];}
28         }swap(x,y);        
29     }
30 }
31 inline void mlink(int x,int y,int pa){
32     while(getfa(x)!=pa){
33         pre[x]=y;y=d[x];
34         if(tp[y]==2){tp[y]=1;q[++tr]=y;}
35         if(getfa(x)==x)fa[x]=pa;
36         if(getfa(y)==y)fa[y]=pa;
37         x=pre[y];
38     }
39 }
40 int doit(int s){
41     for(int i=1;i<=siz;++i)fa[i]=i;
42     memset(tp,0,sizeof(tp));memset(pre,0,sizeof(pre));
43     tl=tr=1;q[1]=s;tp[s]=1;
44     while(tl<=tr){
45         int x=q[tl];++tl;
46         for(int i=head[x];i;i=e[i].next){
47             int y=e[i].y;
48             if(getfa(y)==getfa(x)||tp[y]==2)continue;//cout<<i<<endl;
49             if(!tp[y]){
50                 tp[y]=2;pre[y]=x;//cout<<y<<endl;
51                 if(!d[y]){
52                     for(int now=y,j,las;now;now=las){
53                         j=pre[now];las=d[pre[now]];
54                         d[j]=now;d[now]=j;
55                     }//cout<<i<<1<<endl;
56                     return 1;
57                 }
58                 tp[d[y]]=1;q[++tr]=d[y];
59             }
60             else{
61                 int pa=lca(x,y);
62                 mlink(x,y,pa);mlink(y,x,pa);
63             }
64         }
65     }
66     return 0;
67 }
68 int main(){
69     int T;scanf("%d",&T);
70     while(T-->0){
71         memset(bel,0,sizeof(bel));memset(d,0,sizeof(d));
72         memset(head,0,sizeof(head));
73         tot=tly=0;
74         int p;
75         scanf("%d%d%d",&n,&m,&p);siz=n+m*3;
76         int x,y;
77         for(int i=1;i<=m;++i){
78             init(i,i+m);init(i+m,i);
79             init(i,i+m*2);init(i+m*2,i);
80             init(i+m,i+m*2);init(i+m*2,i+m);
81         }
82         for(int i=1;i<=p;++i){
83             scanf("%d%d",&x,&y);
84             init(m*3+x,y);init(y,m*3+x);
85             init(m*3+x,y+m);init(y+m,m*3+x);
86             init(m*3+x,y+m*2);init(y+m*2,m*3+x);
87         }
88         int ans=0;
89         for(int i=1;i<=siz;++i)if(!d[i])ans+=doit(i);
90         printf("%d\n",ans-n);
91     }
92     return 0;
93 }

View Code