【匈牙利算法】

前提是你 提交于 2019-11-28 08:41:04

前置

二分图:二分图又称作二部图,是图论中的一种特殊模型。 设\(G=(V,E)\)是一个无向图,如果顶点V可分割为两个互不相交的子集\((A,B)\),并且图中的每条边\((i,j)\)所关联的两个顶点i和j分别属于这两个不同的顶点集\((i\;in A,j\;in B)\),则称图G为一个二分图。
简而言之,就是顶点集V可分割为两个互不相交的子集,并且图中每条边依附的两个顶点都分属于这两个互不相交的子集,两个子集内的顶点不相邻。
无向图G为二分图的充分必要条件是,G至少有两个顶点,且其所有回路的长度均为偶数。

匹配:给定一个二分图G,在G的一个子图M中,M的边集{E}中的任意两条边都不依附于同一个顶点,则称M是一个匹配
极大匹配(Maximal Matching)是指在当前已完成的匹配下,无法再通过增加未完成匹配的边的方式来增加匹配的边数。最大匹配(maximum matching)是所有极大匹配当中边数最大的一个匹配。选择这样的边数最大的子集称为图的最大匹配问题。


如果一个匹配中,图中的每个顶点都和图中某条边相关联,则称此匹配为完全匹配,也称作完备匹配。
求二分图最大匹配可以用最大流(Maximal Flow)或者匈牙利算法(Hungarian Algorithm)

实战

下面好几道(好像总共也没几道题)都是网络流的题 因为数据太水可以用匈牙利水过 emmm大概看一下吧 酒店之王有说明

[hdu2063]过山车

匈牙利算法模板

不断地进行匹配,如果当前匹配不可行,则尝试修改原来的匹配,使得最后的匹配数最大。

#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<algorithm> using namespace std; const int N=1000+5,M=500+5,inf=0x3f3f3f3f,P=19650827; int k,n,m,ans,match[N]; double link[N][N],vis[N]; template <class t>void rd(t &x){     x=0;int w=0;char ch=0;     while(!isdigit(ch)) w|=ch=='-',ch=getchar();     while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();     x=w?-x:x; }  bool dfs(int x){     for(int i=1;i<=n;++i){//扫描每个男生          if(link[x][i]&&!vis[i]){//如果能连接 并且i不在当前匈牙利树中              vis[i]=1;             if(!match[i]||dfs(match[i])){match[i]=x;return 1;}             //名花无主或者能腾出空位置来          }      }     return 0; }  int main(){     while(scanf("%d",&k)!=EOF&&k){         ans=0;         memset(link,0,sizeof(link));         memset(match,0,sizeof(match));         rd(m),rd(n);         for(int i=1,x,y;i<=k;++i) rd(x),rd(y),link[x][y]=1;         for(int i=1;i<=m;++i){             memset(vis,0,sizeof(vis));             if(dfs(i)) ++ans;         }         printf("%d\n",ans);     }      return 0; }
dfs版
>queue Q; int prev[__maxNodes]; int Hungarian() {     int ans = 0;     memset(matching, -1, sizeof(matching));     memset(check, -1, sizeof(check));     for (int i=0; i= 0) { // 此点为匹配点                             prev[matching[v]] = u;                         } else { // 找到未匹配点,即找到增广路                             flag = true;                             int d=u, e=v;                             while (d != -1) {//修改原来的匹配关系                                  int t = matching[d];                                 matching[d] = e;                                 matching[e] = d;                                 d = prev[d];                                 e = t;                             }                         }                     }                 }                 Q.pop();             }             if (matching[i] != -1) ++ans;         }     }     return ans; }

[SCOI2010]游戏

luogu1640 bzoj1854

把各个装备值看作左部点 装备编号看作右部点

(还可以用并查集做 具体看黄学长der

vis太大用memset会超时 从标程上学到 用一个数组gai记录变了的vis 然后初始化时只改变了的

#include<bits/stdc++.h> using namespace std #define Max(x,y) ((x)>(y)?(x):(y)) #define Min(x,y) ((x)<(y)?(x):(y)) const int N=1000000+5,M=10000+5,inf=0x3f3f3f3f,P=19650827; int n,mx,ans,match[N],gai[N]; double vis[N]; template <class t>void rd(t &x){     x=0;int w=0;char ch=0;     while(!isdigit(ch)) w|=ch=='-',ch=getchar();     while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();     x=w?-x:x; }  int head[M],tot=0; struct edge{int v,nxt;}e[N<<1]; void add(int u,int v){     e[++tot]=(edge){v,head[u]},head[u]=tot; }  bool dfs(int x){     for(int i=head[x],v;i;i=e[i].nxt)     if(!vis[v=e[i].v]){         vis[gai[++gai[0]]=v]=1;         if(!match[v]||dfs(match[v])){match[v]=x;return 1;}     }     return 0; }  int main(){     freopen("in2.txt","r",stdin);     //freopen("xor.out","w",stdout);     rd(n);mx=ans=0;     for(int i=1,x,y;i<=n;++i) rd(x),rd(y),add(x,i),add(y,i),mx=Max(mx,Max(x,y));     for(int i=1;i<=mx;++i){         for(int i=1;i<=gai[0];++i) vis[gai[i]]=0;gai[0]=0;         if(dfs(i)) ++ans;         else break;     }     printf("%d",ans);     return 0; }

luogu1402酒店之王

一个人有两个要求 若两个要求都满足 则这个人满意

从一个点同时向两边出发 若两个同时满足 则++ans 若不满足 就得把这个点增广出去的边给退回来

但是有bug...主要是这题数据水

(fromluogu讨论区

一组样例:

3 2 2 1 1 1 2 1 1 1 1 1 1 2 1
输出:
2
说明:

中间的节点表示奶牛,左右表示食物和饮料。

还一组样例:
3 2 2 1 1 1 1 1 1 1 2 1 1 2 1
输出:

2
说明:


两个样例输出都是2,但匈牙利算法的第二个程序输出结果为1.

结果的差异本质上是代码思路的问题,可以模拟一下程序运行第二组样例的过程

但改一下,可以把luogu2891水过,让连接边最多的牛先去赠广,这样可以避免发生一些不必要的冲突,例如A食物可以同时给牛a和牛b,但是牛a可以匹配更多食物,我们肯定会将食物给牛a,这样再后续的匹配中,答案更优(但还不是最优)

#include<bits/stdc++.h> using namespace std; #define Max(x,y) ((x)>(y)?(x):(y)) #define Min(x,y) ((x)<(y)?(x):(y)) const int N=100+5,M=10000+5,inf=0x3f3f3f3f,P=19650827; int n,p,q,ans,match1[N],match2[N],M1[N],M2[N]; bool vis[N],vis2[N],lka[N][N],lkb[N][N]; template <class t>void rd(t &x){     x=0;int w=0;char ch=0;     while(!isdigit(ch)) w|=ch=='-',ch=getchar();     while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();     x=w?-x:x; }  bool dfs1(int x){     for(int i=1;i<=p;++i)     if(lka[x][i]&&!vis[i]){         vis[i]=1;         if(!match1[i]||dfs1(match1[i])){match1[i]=x;return 1;}     }     return 0; } bool dfs2(int x){     for(int i=1;i<=q;++i)     if(lkb[x][i]&&!vis2[i]){         vis2[i]=1;         if(!match2[i]||dfs2(match2[i])){match2[i]=x;return 1;}     }     return 0; }  int main(){     //freopen("in2.txt","r",stdin);     //freopen("xor.out","w",stdout);     ans=0;     rd(n),rd(p),rd(q);     for(int i=1;i<=n;++i)     for(int j=1;j<=p;++j) rd(lka[i][j]);     for(int i=1;i<=n;++i)     for(int j=1;j<=q;++j) rd(lkb[i][j]);     for(int i=1;i<=n;++i){         memset(vis,0,sizeof(vis));         memset(vis2,0,sizeof(vis2));         for(int i=1;i<=p;++i) M1[i]=match1[i];         for(int i=1;i<=q;++i) M2[i]=match2[i];         if(dfs1(i)&&dfs2(i)) ++ans;         else{             for(int i=1;i<=p;++i) match1[i]=M1[i];             for(int i=1;i<=q;++i) match2[i]=M2[i];         }     }     printf("%d",ans);     return 0; }

[SCOI2001]小狗散步

luogu2526

每两个相遇点之间只能去一个兴趣点而且每个兴趣点只去一次

将每两个相邻的相遇点与其能去的兴趣点连起

然后从相遇点出发做最大二分图匹配 最后输出答案

#include<bits/stdc++.h> using namespace std; #define ll long long #define Max(x,y) ((x)>(y)?(x):(y)) #define Min(x,y) ((x)<(y)?(x):(y)) const int N=100+5,M=10000+5,inf=0x3f3f3f3f,P=19650827; int n,m,ans=0,match[N]; bool vis[N],lnk[N][N]; template <class t>void rd(t &x){     x=0;int w=0;char ch=0;     while(!isdigit(ch)) w|=ch=='-',ch=getchar();     while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();     x=w?-x:x; }  struct node {     int x,y; }a[N],b[N]; double qdis(node A,node B){return sqrt((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y));}  bool dfs(int x){     for(int i=1;i<=n;++i)     if(lnk[x][i]&&!vis[i]){         vis[i]=1;         if(!match[i]||dfs(match[i])){match[i]=x;return 1;}     }     return 0; }  int main(){     freopen("in2.txt","r",stdin);     //freopen("xor.out","w",stdout);     rd(n),rd(m);     for(int i=1;i<=n;++i) rd(a[i].x),rd(a[i].y);     for(int i=1;i<=m;++i) rd(b[i].x),rd(b[i].y);     for(int i=1;i<n;++i)     for(int j=1;j<=m;++j)     if(qdis(a[i],a[i+1])*2.0>=qdis(a[i],b[j])+qdis(b[j],a[i+1])) lnk[j][i]=1;     for(int i=1;i<=m;++i){         memset(vis,0,sizeof(vis));         if(dfs(i)) ++ans;     }     printf("%d\n",ans+n);     for(int i=1;i<n;++i){         printf("%d %d ",a[i].x,a[i].y);         if(match[i]) printf("%d %d ",b[match[i]].x,b[match[i]].y);     }     printf("%d %d",a[n].x,a[n].y);     return 0; }
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!