【图论】hdu-2255 奔小康赚大钱

这段时间要沉迷刷题一段时间了，就让CSDN陪我一起吧！

一、题目大意

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

输入数据包含多组测试用例，每组数据的第一行输入n,表示房子的数量(也是老百姓家的数量)，接下来有n行,每行n个数表示第i个村名对第j间房出的价格(n<=300)。

请对每组数据输出最大的收入值，每组的输出占一行。

二、题目思路以及AC代码

这题吧，说难也难，说简单也简单。因为这题一看完题目就知道是一道要用KM算法求解的，要求二分图完备匹配下的最大权匹配。但难就难在，这个KM算法不好理解呀！

因为我也是KM算法的初学者，所以也就不献丑了，不过我的解题代码大家还是可以看看的，我为了方便大家理解，在上面进行了详细的标注，也是给自己理清一下思路。如果大家要详细了解KM算法的话，我可以给大家推荐一个博客，可以作为参考。

KM算法详解+模板

下面给出我注释的详细代码，当然也是hdu的AC代码啦！希望大家可以在看完之前那个博客，对KM算法有所了解之后再阅读我的代码及注释，不然我相信就算有注释，你也是看不懂的0.0

#include <iostream>
#define MAXN 305
#define INF 10000000
using namespace std;

int N;

int lines[MAXN][MAXN];			// 记录二分图的边权
int f_people[MAXN];				// 各个村民的期望值(标杆)，初始为连接边中的最大值
int f_house[MAXN];				// 各个房屋的期望值(标杆)，初始为0
bool vis_people[MAXN];			// 标记在此次匹配中是否被访问过
bool vis_house[MAXN];			// 同上
int match[MAXN];				// match[i]记录了i号房屋所归属的村民号
								// i号房屋已经归属于match[i]，为0表示无归属
int slack[MAXN];				// 记录L(x)+L(y)-lines[x][y] 用于计算减少的最小权值

int Min(int a, int b) {			// 最小值函数
	return a < b ? a : b;
}

int Max(int a, int b) {			// 最大值函数
	return a > b ? a : b;
}


// 初始化函数
void init() {			
	for (int i = 0; i < MAXN; i++) {
		for (int j = 0; j < MAXN; j++) {
			lines[i][j] = 0;
		}
		f_people[i] = 0;
		f_house[i] = 0;
		vis_house[i] = false;
		vis_people[i] = false;
		match[i] = 0;
		slack[i] = INF;
	}
}

/*
 * 匹配函数
 * 输入参数： 
 *			x：表示为x号村民查找匹配房屋
 * 返回值：为x号村民匹配成功返回true，失败返回false
 */
bool find(int x) {
	// 标记x号村民已访问，用于记录增广路径，从而更新期望值(标杆值)f_people
	vis_people[x] = true;

	// 遍历每一个房屋，找满足要求的，即权值损失最少
	for (int i = 1; i <= N; i++) {
		// 如果已经遍历过，跳过，表明在此次匹配中已经有村民预定了这个房屋
		if (vis_house[i]) continue;
		// 计算优化值 L(x)+L(y)-lines[x][y]
		int gap = f_people[x] + f_house[i] - lines[x][i];

		// gap等于0，说明此时是权值损失最小的时候
		if (gap == 0) {
			// 预定该房屋，防止其他村民购买
			vis_house[i] = true;

			// 第一个条件是如果该房屋暂时没有主人
			// 第二个条件是如果该房屋的主人可以找到别的房子，把这个房子让出来
			if (match[i] == 0 || find(match[i])) {
				// 那就占有这个房子，匹配成功
				// (此时该村民之前的村民会有房屋更换的情况，但也是匹配好房屋的)
				match[i] = x;
				return true;
			}
		}
		// 如果没满足条件，记录最优值，为后续求最小值准备
		// 这里slack[i]，是i号房子经历所有后的最优值最小值，外面再对所有房屋求最小就是全部的最小
		// 可以理解一下，不论怎样，slack最终求得的肯定是全部的最小
		else {
			slack[i] = Min(gap, slack[i]);
		}
	}
	return false;
}

/*
 * KM算法
 * 输入参数：
 *			None
 * 返回值：二分图完备匹配的最大权值
 */
int KM() {

	// 初始化村民的期望值(标杆值)
	for (int i = 1; i <= N; i++) {
		f_people[i] = lines[i][0];
		for (int j = 1; j <= N; j++) {
			f_people[i] = Max(f_people[i], lines[i][j]);
		}
	}

	// 初始化房屋的期望值(标杆值)，全0
	for (int i = 1; i <= N; i++) {
		f_house[i] = 0;
	}

	// 为每一个村民匹配房屋，从1 ~ N
	for (int i = 1; i <= N; i++) {
		
		// 因为后来要求最小值，所以置为INF
		// 如果不每次重置的话，前面的匹配遗留会影响后面
		for (int j = 1; j <= N; j++) {
			slack[j] = INF;
		}

		while (1) {
			// 每次匹配都要重置，因为每次要寻找增广路径
			for (int j = 0; j < MAXN; j++) {
				vis_house[j] = false;
				vis_people[j] = false;
			}

			// 如果匹配成功，则跳出，开始下一个村民
			if (find(i)) break;

			// 如果匹配没成功，则降低期望值
			int d = INF;
			for (int j = 1; j <= N; j++) {
				// 这里与上述求最小值对应，求全部的最小d
				if(!vis_house[j]) d = Min(d, slack[j]);
			}

			// 对于在增广路径上的村民点，减少d
			for (int j = 1; j <= N; j++) {
				if (vis_people[j]) f_people[j] -= d;
			}

			// 对于在增广路径上的房屋点，增加d
			// 对于不在增广路上的房屋，由于村民点减少d，所以其需要的差值也减少d
			for (int j = 1; j <= N; j++) {
				if (vis_house[j]) f_house[j] += d;
				else slack[j] -= d;
			}
		}
	}

	// 求解最大权值
	// 此时以为每个房屋i分配了归属match[i]
	int res = 0;
	for (int i = 1; i <= N; i++) {
		res += lines[match[i]][i];
	}

	return res;
}

int main() {
	
	while (scanf("%d", &N) != EOF) {
		init();

		for (int i = 1; i <= N; i++) {
			for (int j = 1; j <= N; j++) {
				scanf("%d", &lines[i][j]);
			}
		}

		printf("%d\n", KM());
	}

	return 0;
}