Uva11383 二分圖的完美匹配(深入理解KM演算法)
https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2378
/**
Uva11383 二分圖的完美匹配(深入理解KM演算法)
劉汝佳大白書P351 少林決勝
*/
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <iostream>
using namespace std;
const int maxn=510;
const int inf=0x3f3f3f3f;
int n,nx,ny;
int link[maxn],lx[maxn],ly[maxn],slack[maxn];
int visx[maxn],visy[maxn],w[maxn][maxn];
int DFS(int x)
{
visx[x]=1;
for(int y=1;y<=ny;y++)
{
if(visy[y])continue;
int t=lx[x]+ly[y]-w[x][y];
if(t==0)
{
visy[y]=1;
if(link[y]==-1||DFS(link[y]))
{
link[y]=x;
return 1;
}
}
else if(slack[y]>t)
{
slack[y]=t;
}
}
return 0;
}
int KM()
{
memset(link,-1,sizeof(link));
memset(ly,0,sizeof(ly));
for(int i=1;i<=nx;i++)
{
lx[i]=-inf;
for(int j=1;j<=ny;j++)
{
if(w[i][j]>lx[i])
lx[i]=w[i][j];
}
}
for(int x=1;x<=nx;x++)
{
for(int i=1;i<=ny;i++)
slack[i]=inf;
while(1)
{
memset(visx,0,sizeof(visx));
memset(visy,0,sizeof(visy));
if(DFS(x))break;
int d=inf;
for(int i=1;i<=ny;i++)
{
if(!visy[i]&&d>slack[i])
d=slack[i];
}
for(int i=1;i<=nx;i++)
{
if(visx[i])
lx[i]-=d;
}
for(int i=1;i<=ny;i++)
{
if(visy[i])
ly[i]+=d;
else
slack[i]-=d;
}
}
}
int res=0;
for(int i=1;i<=ny;i++)
{
if(link[i]>-1)
res+=w[link[i]][i];
}
return res;
}
int main()
{
while(~scanf("%d",&n))
{
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
scanf("%d",&w[i][j]);
}
}
nx=n,ny=n;
int ans=KM();
// printf("KM=%d\n",ans);
int sum=0;
for(int i=1;i<=n;i++)
{
printf("%d%c",lx[i]," \n"[i==n]);
sum+=lx[i];
}
for(int i=1;i<=n;i++)
{
printf("%d%c",ly[i]," \n"[i==n]);
sum+=ly[i];
}
printf("%d\n",sum);
}
return 0;
}
相關文章
- poj2400 KM演算法二分圖的完美匹配演算法
- KM演算法——二分圖的最佳匹配演算法
- 二分圖的最大匹配、完美匹配和匈牙利演算法演算法
- hdu2255 二分圖的最佳匹配 KM演算法演算法
- POJ 3565 Ants (最小權完美匹配 KM演算法)演算法
- 二分圖最大權完美匹配
- 圖論-二分圖匹配匈牙利演算法圖論演算法
- 匈牙利演算法--二分圖的最大匹配演算法
- 對匈牙利演算法理解——對二分圖進行最大匹配的演算法演算法
- 二分圖的最大匹配的匈牙利演算法演算法
- 二分圖最大匹配(匈牙利演算法)演算法
- 二分圖匹配
- 二分圖的最大匹配(匈牙利演算法)程式碼演算法
- 匹配字串之——KMP演算法深入理解字串KMP演算法
- 詳解匈牙利演算法與二分圖匹配演算法
- 二分圖最大匹配問題匈牙利演算法演算法
- 對KM演算法暫時性的理解演算法
- 《啊哈!演算法》我要做月老 ——二分圖最大匹配演算法
- HDU 2063 匈牙利演算法二分圖的最大匹配演算法
- POJ 1325-Machine Schedule(二分圖匹配-匈牙利演算法)Mac演算法
- UVA 11383 - Golden Tiger Claw【二分圖km原理】Go
- 演算法學習之路|二分圖的最大匹配—匈牙利演算法(Dfs實現)演算法
- 二分圖最小點覆蓋等於二分圖最大匹配
- KM演算法演算法
- POJ 1469-COURSES(二分圖匹配入門-匈牙利演算法)演算法
- hdu5090 匈牙利演算法二分圖最大匹配問題演算法
- HDU 2255-奔小康賺大錢(Kuhn-Munkras演算法/KM演算法-完備匹配下的最大權匹配)演算法
- 深入理解二分的兩種實現
- POJ - 3041 Asteroids 【二分圖匹配】AST
- POJ 3041-Asteroids(二分圖匹配)AST
- KM演算法小記演算法
- KMP字串匹配演算法 通俗理解KMP字串匹配演算法
- KMP字串匹配演算法理解(轉)KMP字串匹配演算法
- 匈牙利演算法模板(二分圖)演算法
- 學習筆記----KM演算法筆記演算法
- 深入理解React:diff 演算法React演算法
- 深入理解React Diff演算法React演算法
- 地圖匹配演算法實踐地圖演算法