HDU-2255
题目意思转化之后就是,给你一个二分图(也称 二部图) ,要求选择一些边让左边的点都对应左边的某一个点!该问题也叫做二分图最大匹配。所以可以用KM算法来做这道题。KM前提你要理解匈牙利算法,最大二分匹配问题。
所以先简单阐述一下KM 算法过程:定义连个点集合A,B两个集合 定义A与B之间的边为 E
(1) 初始化顶标 lx[] , ly[] 两个数组;
● lx[i]初始化为A集合中 i 点能到B集合某一点的最大权值,
● ly[i] 初始化0;
(2) 用匈牙利找最大匹配;
(3) 如果找到了 进行步骤(4);
否则 扩边操作,回到步骤(2);
(4) 根据匈牙利算法中的二分图连接,求出最大权值!
注释:为什么要用 lx,ly数组? KM有着贪心的思想,一开始最大匹配时都找每个点的最大权值边,不行再把要求放稍微低一点,再来进行最大匹配!所以,才需要lx,ly数组,以及扩边操作。lx 我们暂且叫做期望权值
KM算法中 在匈牙利算法那一部分加了一个条件 假设x→y ,则需要 (lx[x]+ly[y])==(value:x→y)
因为一开始lx都是最大权边,ly为0,在进行匈牙利算法之中选择时候,就只会选择指定的边,如果不能够 ,这样做的母的是筛选出边,哪一些对于左图点权值较大的边,这样第一选择一定都是自己点所能到达的最大权值。如果匹配中断了,有点无法匹配,就需要下降要求,就是将 lx 减去某个数(这个数字,是稍微降低已经匹配好的点的期望权值,稍微降低的意思就是:使得减少量尽可能最少);
如果还是不好理解,推荐这一篇博客;
#include <cstdio>
#include <cstring>
#include <cctype>
#include <cmath>
#include <set>
#include <map>
#include <list>
#include <queue>
#include <deque>
#include <stack>
#include <string>
#include <vector>
#include <iostream>
#include <algorithm>
#include <stdlib.h>
#include <time.h>
using namespace std;
typedef long long LL;
const int INF=0x3f3f3f3f;
const int MOD=1e9+7;
const int MAXSIZE=1e6+5;
const double eps=0.0000000001;
void fre()
{
freopen("in.txt","r",stdin);
freopen("out.txt","w",stdout);
}
#define memst(a,b) memset(a,b,sizeof(a))
#define fr(i,a,n) for(int i=a;i<n;i++)
const int MAXN=305;
int adj[MAXN][MAXN],n;
int lx[MAXN],ly[MAXN],link[MAXN];
bool visx[MAXN],visy[MAXN];
bool dfs(int x) // 匈牙利算法部分
{
visx[x]=1;
for(int i=1;i<=n;i++)
{
if(adj[x][i]==lx[x]+ly[i]&&visy[i]==0)
{
visy[i]=1;
if(link[i]==-1||dfs(link[i]))
{
link[i]=x;
return true;
}
}
}
return false;
}
int KM()
{
memset(link,-1,sizeof(link));
for(int i=1;i<=n;i++)
{
lx[i]=-INF,ly[i]=0;
for(int j=1;j<=n;j++) lx[i]=max(lx[i],adj[i][j]);
}
for(int k=1;k<=n;k++)//匈牙利算法
{
while(1)//对于每个左图的点,进行不断的查找最大权边,一旦发现有
// 某一个左图点不能匹配,要求下降一点;继续
{
memst(visy,0);
memst(visx,0);
if(dfs(k)) break;
int minval=INF;
for(int i=1;i<=n;i++) if(visx[i]) for(int j=1;j<=n;j++) if(!visy[j]) minval=min(minval,lx[i]+ly[j]-adj[i][j]);
if(minval==INF) return -1;
for(int i=1;i<=n;i++) if(visx[i]) lx[i]-=minval;
for(int i=1;i<=n;i++) if(visy[i]) ly[i]+=minval;
}
}
int res=0;
for(int i=1;i<=n;i++)
{
if(link[i]!=-1) res+=adj[link[i]][i];
}
return res;
}
int main()
{
while(scanf("%d",&n)+1)
{
memset(adj,-1,sizeof(adj));
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
scanf("%d",&adj[i][j]);
printf("%d\n",KM());
}
return 0;
}
/**************************************************/
/** Copyright Notice **/
/** writer: wurong **/
/** school: nyist **/
/** blog : http://blog.csdn.net/wr_technology **/
/**************************************************/