HDU-2255(KM算法)

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  **//**************************************************/