hdu 2255奔小康赚大钱 KM算法模板

题目链接:http://acm.hdu.edu.cn/showproblem.php?

pid=2255

一,KM算法:(借助这个题写一下个人对km的理解与km模板)

KM算法主要是用来求解图的最优匹配的。

1。带权二分图:  在二分图中每一条边(x。y)相应一个权值Wi这样的二分图叫带权二分图。

一个匹配的权值就是该匹配中全部边的权值之和。

2,最优匹配:

权值最大的一个完美匹配。叫做最优匹配。

     
《km算法思想》

对于一个带权全然二分图:G(V。E),对于当中每一条边(x。y)边都相应一个权值W(x,y)。

因为是二分图所以节点集合能够分解成两个集合(X。Y)。X={x1。x2。,,。。xn},Y={y1,y2。,。。yn}。

可行性标签:我们将在全部顶点上定义一个实函数F,且对全部定点满足:F(x)+ F(y)>=W(x,y),

则我们称F为图G的一个可行性标签。

(这里能够理解成F是对顶点加权值F(x))

平庸标签:对于每个图都有存在这样一个标签 l:

称该标签为平庸标签

相等子图:对于每个标签都相应原图的一个子图:G’=(V’,E’)  ,当中V’=V。E’属于E,

且E’中全部边满足:F(x) + F(y)=W(x,y) 。该子图称为可行性标签的F的相等子图。

定理:(KM最重要的定理)

对于一个可行性标签L,其相等子图G’,若存在完美匹配M。则该匹配必是原图G的完美匹配。

由于G中与G’中的节点是全然一样的,且对于该匹配的权值为:

所以对于最大的可行性标签L,若它的相等子图G‘含有完美匹配M,则M就是原图的最优匹配。

注意:上面的平庸匹配就是全部匹配中最大的匹配。

匈牙利树:

匈牙利树。事实上是一颗DFS搜索树。

可是在搜索是有条件限制。

(1)。该树必需要从未匹配点開始搜索。即,该树必须以未匹配点为根。

(2)。并且在搜索过程中仅仅能走叫错路。

(3)。叶子节点必须是匹配点

注意:若有叶子节点为未匹配点则从树根到该节点一定有一条增广路。

事实上在匈牙利算法中,若我们从一个未匹配的点。利用DFS找增广路。若到最后我们没找到增广路。

则我们DFS搜索过程的路线就是一颗DFS树且是一颗匈牙利树,且在匈牙利算法中树我们能够永久

性的地把匈牙利树从图中删去,而不影响结果。即:这颗树对树以外的节点的匹配没有影响。

在上图中(2)不是匈牙利树由于有一个叶子节点7,非常明显有一条2到7的增广路。

(3)是一颗匈牙利树。

更新可行性标签:

当我们在当前大的可行性标签l下。无法在L相等子图中找到完美匹配,

这时我们要更新可行性标签,

事实上就是在原来的可行性标签l下降低最小的一个值(key)得到还有一个可行性标签L。

更新方程为:

       《KM算法》

如果二分图为:G=(V,E),(X。Y)= V。X={x1。x2。,,。,xn}。Y={y1,y2,,,,yn}。

1。初始化可行性标签L为平庸标签(最大标签)。

2。在可行性标签相等子图G’=(V。E’)中用匈牙利算法(DFS)对X集合中节点一个个进行匹配。

若在匹配过程中,若有X集合中的点没匹配成功。转向(3)

(3)匹配失败这时DFS索搜出的路径必然是一颗匈牙利树。

由匈牙利树性质可知匈牙利树中的节点与边与匈牙利树外部是相互独立的,所以我们改变匈牙利

树种节点的可行新标签就可以。且不会对其他树外有影响,于是我们利用上面标签转化公式改变

匈牙利中节点的可行性标签,然后再转向(2)。

Km 例题:

令:图节点为:V1={x1,x2,x3,x4,x5}  。V2={y1。y2,y3。y4,y5}。边矩阵为:

初始化可行性标签(平庸标签)

L(y1)= L(y2) = L(y3) = L(y4) = L(y5) = 0,

L(x1)=max(3,5,5,4,1)=5 , L(x2)=max(2,2,0,2,2) =2  ,L(x3)=max(2,4,4,1,0) =4

L(x4)=max(0,1,1,0,0) =1  。L(x4)=max(0,1,1,0,0) =1 。

L的相等子图为:

L相等子图的匹配(匹配到x4点失败)

x4搜索出的匈牙利树为:

则 S={x4,x1,x3}。T={y2,y3} 利用可行性标签变换公式为L’:

可得到:key=1;改变:

L’(x1)=4。L’(x2)= 2, L’(x4)=0 ;   L’(y2)=1。L’(y3)=1。

然后x4继续匹配:能够匹配完:

这上图L’完美匹配也就是原图最优匹配。

W(M)=12。

代码:

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
using namespace std;
const int Max=310;
const int Inf=(1<<31)-1;
int N;
int graph[Max][Max];//图
int match[Max];//行匹配
int l_x[Max];//x可行性标签
int l_y[Max];//y可行性标签
int usedx[Max];//标记数组
int usedy[Max];

bool Hungary(int x)
{//匈牙利DFS
    usedx[x]=true;
    for(int i=0;i<N;i++)
    {
        if(!usedy[i] && l_x[x]+l_y[i]==graph[x][i])
        {
            usedy[i]=true;
            if(match[i]==-1 || Hungary(match[i]))
            {
                match[i]=x;
                return true;
            }
        }
    }
    return false;
}

void KM()
{
    memset(l_x,0,sizeof(l_x));
    memset(l_y,0,sizeof(l_y));
    for(int i=0; i<N; i++)
    {//初始化可信性标签
        for(int j=0; j<N; j++)
        {
            if(l_x[i]<graph[i][j])
                l_x[i]=graph[i][j];
        }
    }
    for(int i=0;i<N;i++)
    {//匹配V1所以点
        while(1)
        {
            memset(usedx,false,sizeof(usedx));
            memset(usedy,false,sizeof(usedy));
            if(Hungary(i))
                break;//匹配成功就跳出
            else
            {//匹配失败 匈牙利树节点更新
                int tem=Inf;//求最小的 key
                for(int j=0;j<N;j++)
                {//用到那公式
                    if(usedx[j])
                    {
                        for(int k=0;k<N;k++)
                        {
                            if(!usedy[k] && l_x[j]+l_y[k]-graph[j][k]<tem)
                                 tem= l_x[j]+l_y[k]-graph[j][k];
                        }
                    }
                }
                for(int j=0;j<N;j++)
                {//更新匈牙利树可行性标签
                    if(usedx[j])
                        l_x[j]-=tem;
                    if(usedy[j])
                        l_y[j]+=tem;
                }
            }
        }
    }
}

int main()
{
    while(scanf("%d",&N)!=EOF)
    {
        for(int i=0; i<N; i++)
            for(int j=0; j<N; j++)
                scanf("%d",&graph[i][j]);
        memset(match,-1,sizeof(match));
       KM();
       int total=0;
       for(int i=0;i<N;i++)
       {
           total+=l_x[i];
           total+=l_y[i];
       }
       cout<<total<<endl;
    }
    return 0;
}

时间: 2024-10-14 07:43:05

hdu 2255奔小康赚大钱 KM算法模板的相关文章

hdu 2255 奔小康赚大钱 (km算法模板)

题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=2255 解题思路: 了解km算法以后,就可以直接套用km算法,km算法:完备匹配下的最大权匹配, 代码: 1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <algorithm> 5 using namespace std; 6 7 #define maxn 310

HDU 2255 奔小康赚大钱 KM算法的简单解释

KM算法一般用来寻找二分图的最优匹配. 步骤: 1.初始化可行标杆 2.对新加入的点用匈牙利算法进行判断 3.若无法加入新编,修改可行标杆 4.重复2.3操作直到找到相等子图的完全匹配. 各步骤简述: 1.根据二分图建立2个可行标杆; lx为x的可行标杆,初始化lx[i]为与i点相连的最大边 ly为y的可行标杆,初始化为0. 可行性的判断条件应为lx[x]+ly[y] >= Map[x][y]. 2.对于新加入的点用匈牙利算法经行判断,确定改点能否加入旧子图中,形成新子图. 对于该加入的点有两种

HDU 2255 奔小康赚大钱(KM算法求完备匹配)

Problem Description: 传说在遥远的地方有一个非常富裕的村落,有一天,村长决定进行制度改革:重新分配房子.这可是一件大事,关系到人民的住房问题啊.村里共有n间房间,刚好有n家老百姓,考虑到每家都要有房住(如果有老百姓没房子住的话,容易引起不安定因素),每家必须分配到一间房子且只能得到一间房子.另一方面,村长和另外的村领导希望得到最大的效益,这样村里的机构才会有钱.由于老百姓都比较富裕,他们都能对每一间房子在他们的经济范围内出一定的价格,比如有3间房子,一家老百姓可以对第一间出1

HDU 2255 奔小康赚大钱 KM算法题解

KM算法求的是完备匹配下的最大权匹配,是Hungary算法的进一步,由于Hungary算法是最大匹配的算法,不带权. 经典算法,想不出来的了,要參考别人的.然后消化吸收吧. 由于真的非常复杂的算法. 我理解算法匹配思想: 1 開始的时候,全部边都记录自己的最优匹配,无论有没有冲突 2 递归循环的时候.假设找不到自己的最优匹配,那么就找次要匹配 3 次要匹配不行,继续找下一个次优匹配,全部点都一定要找到解 难点: 怎样记录这些信息,以便循环处理全部点. 牵涉到非常多什么增广路,交错树之类的,名词,

HDU 2255 奔小康赚大钱 KM裸题

#include <stdio.h> #include <string.h> #define M 310 #define inf 0x3f3f3f3f int n,nx,ny; int link[M],lx[M],ly[M],slack[M]; //lx,ly为顶标,nx,ny分别为x点集y点集的个数 int visx[M],visy[M],w[M][M]; int DFS(int x) { visx[x] = 1; for (int y = 1;y <= ny;y ++)

hdu 2255 奔小康赚大钱(KM算法)

奔小康赚大钱 Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 3248    Accepted Submission(s): 1413 Problem Description 传说在遥远的地方有一个非常富裕的村落,有一天,村长决定进行制度改革:重新分配房子. 这可是一件大事,关系到人民的住房问题啊.村里共有n间房间,刚好有n家老百姓,考

hdu2255 奔小康赚大钱 km算法解决最优匹配(最大权完美匹配)

/** 题目:hdu2255 奔小康赚大钱 km算法 链接:http://acm.hdu.edu.cn/showproblem.php?pid=2255 题意:lv 思路:最优匹配(最大权完美匹配) km算法 模板来自:http://www.cnblogs.com/wenruo/p/5264235.html 如果是求最小权完美匹配,那么将所有权值取相反数,然后求得最大权,输出最大权的相反数即可. */ #include <iostream> #include <cstring> #

HDU 2255 ——奔小康赚大钱——————【KM算法裸题】

奔小康赚大钱 Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u Submit Status Practice HDU 2255 Description 传说在遥远的地方有一个非常富裕的村落,有一天,村长决定进行制度改革:重新分配房子. 这可是一件大事,关系到人民的住房问题啊.村里共有n间房间,刚好有n家老百姓,考虑到每家都要有房住(如果有老百姓没房子住的话,容易引起不安定因素),每家必须分配到

[ACM] HDU 2255 奔小康赚大钱 (二分图最大权匹配,KM算法)

奔小康赚大钱 Problem Description 传说在遥远的地方有一个很富裕的村落,有一天,村长决定进行制度改革:又一次分配房子. 这但是一件大事,关系到人民的住房问题啊. 村里共同拥有n间房间,刚好有n家老百姓,考虑到每家都要有房住(假设有老百姓没房子住的话.easy引起不安定因素),每家必须分配到一间房子且仅仅能得到一间房子. 还有一方面,村长和另外的村领导希望得到最大的效益,这样村里的机构才会有钱.因为老百姓都比較富裕,他们都能对每一间房子在他们的经济范围内出一定的价格,比方有3间房