Dinic二分图匹配 || Luogu P3386

Posted 2021-12-14 alenanuna

题面：【模板】二分图匹配

思路：Dinic实现二分图匹配，要建一个超级源点（S）和超级汇点（T），分别定为N+M+1和N+M+2

然后S去和N中的数建正边和反边，正边权值为1，反边权值为0；M中的数去与T建正边和反边，正边权值为1。

N、M之间的数建图一样。

然后就去跑最大流。

注意：在Dinic函数中每次更新Cur的值时，要把S和T的Cur也更新了。

代码：

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<iostream>
 4 #define min(a,b) ((a)<(b)?(a):(b))
 5 using namespace std;
 6 const int maxn=1050,maxm=maxn,maxe=maxm*maxn;
 7 int N,M,u,v,E,num_edge=-1,edge_head[maxn+maxm],Q[maxn+maxm],f1,f2,Dep[maxn+maxm];
 8 int S,T,Cur[maxn+maxm];
 9 struct Edgeint to,nx,dis;edge[maxe];
10 inline void Add_edge(int from,int to,int dis)
11     edge[++num_edge].nx=edge_head[from];
12     edge[num_edge].to=to;
13     edge[num_edge].dis=dis;
14     edge_head[from]=num_edge;
15     return;
16 
17 inline bool Bfs()
18     memset(Dep,0,sizeof(Dep));
19     f1=f2=1;
20     Dep[S]=1;
21     Q[f2++]=S;
22     while(f1<f2)
23         int x=Q[f1++];
24         for(int i=edge_head[x];i!=-1;i=edge[i].nx)
25             int y=edge[i].to;
26             if(edge[i].dis&&Dep[y]==0)
27                 Dep[y]=Dep[x]+1;
28                 Q[f2++]=y;
29             
30         
31     
32     if(Dep[T])return 1;
33     return 0;
34 
35 inline int Dfs(int x,int fw)
36     if(x==T)return fw;
37     for(int &i=Cur[x];i!=-1;i=edge[i].nx)
38         int y=edge[i].to;
39         if(Dep[y]==Dep[x]+1&&edge[i].dis)
40             int p=Dfs(y,min(edge[i].dis,fw));
41             if(p>0)
42                 edge[i].dis-=p;
43                 edge[i^1].dis+=p;
44                 return p;
45             
46         
47     
48     return 0;
49 
50 inline int Dinic()
51     int ans=0;
52     while(Bfs())
53         int toi=N+M+2;
54         for(int i=1;i<=toi;i++)Cur[i]=edge_head[i];
55         while(int k=Dfs(S,1<<30))ans+=k;
56     
57     return ans;
58 
59 int main()
60     memset(edge_head,-1,sizeof(edge_head));
61     scanf("%d%d%d",&N,&M,&E);
62     S=N+M+1;T=N+M+2;
63     for(int i=1;i<=N;i++)
64         Add_edge(S,i,1);
65         Add_edge(i,S,0);
66     
67     int toi=N+M;
68     for(int i=N+1;i<=toi;i++)
69         Add_edge(i,T,1);
70         Add_edge(T,i,0);
71     
72     for(int i=1;i<=E;i++)
73         scanf("%d%d",&u,&v);
74         if(v>M||u>N)continue;
75         v+=N;
76         Add_edge(u,v,1);
77         Add_edge(v,u,0);
78     
79     printf("%d\n",Dinic());
80     return 0;
81 

---

By:AlenaNuna