二分匹配的Hopcroft-Carp算法
生活随笔
收集整理的這篇文章主要介紹了
二分匹配的Hopcroft-Carp算法
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
題目:HDU2063 過山車
?
果然Hopcroft-Carp算法快,用匈牙利算法15ms,而Hopcorft-Carp卻0ms。因為匈牙利算法的時間復雜度為O(n*e),而Hopcroft-Carp算法O(sqrt(n)*e)
本算法適合處理大一點的數據。。。。。。
?
#include <stdio.h> #include <iostream> #include <string.h> #include <queue>const int N=1005; const int INF=1<<28;int g[N][N]; int Mx[N]; int My[N]; int dx[N]; int dy[N]; bool used[N];int Nx,Ny,dis;bool searchP() {dis=INF;int i,v,u;std::queue<int> Q;memset(dx,-1,sizeof(dx));memset(dy,-1,sizeof(dy));for(i=0;i<Nx;i++){if(Mx[i]==-1){Q.push(i);dx[i]=0;}}while(!Q.empty()){u=Q.front();Q.pop();if(dx[u]>dis) break;for(v=0;v<Ny;v++){if(g[u][v]&&dy[v]==-1){dy[v]=dx[u]+1;if(My[v]==-1) dis=dy[v];else{dx[My[v]]=dy[v]+1;Q.push(My[v]);}}}}return dis!=INF; }bool DFS(int u) {int v;for(v=0;v<Ny;v++){if(g[u][v]&&!used[v]&&dy[v]==dx[u]+1){used[v]=true;if(My[v]!=-1&&dy[v]==dis) continue;if(My[v]==-1||DFS(My[v])){My[v]=u;Mx[u]=v;return true;}}}return false; }int Hungary() {int u;int ret=0;memset(Mx,-1,sizeof(Mx));memset(My,-1,sizeof(My));while(searchP()){memset(used,false,sizeof(used));for(u=0;u<Nx;u++)if(Mx[u]==-1&&DFS(u)) ret++;}return ret; }int main() {int k,u,v;while(~scanf("%d",&k),k){scanf("%d%d",&Nx,&Ny);memset(g,0,sizeof(g));Ny=Nx>Ny? Nx:Ny;while(k--){scanf("%d%d",&u,&v);u--;v--;g[u][v]=1;}int ans=Hungary();printf("%d\n",ans);}return 0; }
?
總結
以上是生活随笔為你收集整理的二分匹配的Hopcroft-Carp算法的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: POJ2594(二分匹配+Floyd求传
- 下一篇: POJ2154(Pólya定理与欧拉函数