HDU2295 Radar (DLX)

 

下面的代码99％参考了这个网站http://www.cnblogs.com/183zyz/archive/2011/08/07/2130193.html

人生的第一道DLX肯定是需要作一些参考的啦。

题意：给你N个城市，M个雷达，你要在其中选K个，问当半径最小是多少的时候可以覆盖到所有的N个城市。

做法：二分所需要的半径r,然后处理出这时雷达能够覆盖到哪些位置。然后就转化成了一个重复覆盖问题，跑一下DLX即可。

花了两天的时间学习了一下DLX这种神奇的数据结构,它在解决精确覆盖问题上尤其有帮助，因为随着每次选的一行，对应的包含了这个行的列的行也不会再用到，所以整个图会很快的变得非常小，从而可以极大得加快搜索速度。但是重复覆盖就没有缩小得那么快，包含了对应的列的行是不会消除掉的。

下面对模板的一些数组说下理解  L,R,U,D就是链表的四个方向,其中H是行的头指针，当该行不存在东西的时候H[]=-1,而col[size]表示结点编号为size所在的是哪一列，在精确覆盖的时候还要多维护一个row[size]去表示是在哪一行。对应的link函数是用来加条件的，然后dance就是递归，remove和resume则是对应的删除和恢复，在重复覆盖和精确覆盖问题上这两个操作要做相应的修改。

题目了用了个A*，多亏我有数据结构大作业的基础，A*就是在搜的时候多了个估价函数，但是这个估价函数必须保证不能高估了可能的解，譬如解是6的话你不能估成7，这样的估价函数是不行的。下面给的h函数是一定估小的，虽然我觉得可能优化的不是特别的大。

细节等别的地方后面再慢慢学习，先照着别人的拍一题加深一下认识。。

#pragma warning(disable:4996)

#include<iostream>

#include<cstring>

#include<string>

#include<vector>

#include<cstdio>

#include<cmath>

#include<algorithm>

#include<queue>

using namespace std;



#define maxnode 4000

#define maxn 80

#define eps 1e-8

int L[maxnode], R[maxnode], U[maxnode], D[maxnode];

int col[maxnode]; 

int S[maxn], H[maxn];

int size;

int n, m, K;



int ak; // currently best dep



void link(int r, int c)

{

	S[c]++; col[size] = c;

	U[size] = U[c]; D[U[c]] = size;

	D[size] = c; U[c] = size;



	if (H[r] == -1) H[r] = L[size] = R[size] = size;

	else{

		L[size] = L[H[r]]; R[L[H[r]]] = size;

		R[size] = H[r]; L[H[r]] = size;

	}

	size++;

}



void remove(int c)

{

	for (int i = D[c]; i != c; i = D[i]){

		L[R[i]] = L[i]; R[L[i]] = R[i];

	}

}



void resume(int c)

{

	for (int i = U[c]; i != c; i = U[i]){

		L[R[i]] = R[L[i]] = i;

	}

}

// evaluating function

int hfunc()

{

	bool vis[maxn];

	memset(vis, 0, sizeof(vis));

	int cnt = 0;

	for (int i = R[0]; i; i = R[i]){

		if (vis[i]) continue;

		cnt++; vis[i] = true;

		for (int j = D[i]; j != i; j = D[j]){

			for (int k = R[j]; k != j; k = R[k]){

				vis[col[k]] = true;

			}

		}

	}

	return cnt;

}



void dance(int dep)

{

	int ans = hfunc();

	if (ans + dep > K || ans + dep >= ak) return;

	if (R[0] == 0){

		if (dep < ak) ak = dep; return;

	}

	int minv = maxn; int c;

	for (int i = R[0]; i; i = R[i]){

		if (S[i] < minv) minv = S[i], c = i;

	}

	//枚举删除掉的每一行

	for (int i = D[c]; i != c; i = D[i]){

		remove(i);

		// 将该行上的对应的列也删掉

		for (int j = R[i]; j != i; j = R[j]){

			remove(j);

		}

		dance(dep + 1);

		for (int j = L[i]; j != i; j = L[j]){

			resume(j);

		}

		resume(i);

	}

	return;

}



double x[100], y[100];

double xx[100], yy[100];



int dcmp(double x){

	return (x > eps) - (x < -eps);

}



double dist(double x1, double y1, double x2, double y2){

	return sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2));

}



int main()

{

	int T; cin >> T;

	while (T--)

	{

		scanf("%d%d%d", &n, &m, &K);

		for (int i = 1; i <= n; i++){

			scanf("%lf%lf", &x[i], &y[i]);

		}

		for (int i = 1; i <= m; i++){

			scanf("%lf%lf", &xx[i], &yy[i]);

		}

		double l = 0, r = 1500;

		while (dcmp(r - l) > 0){

			for (int i = 0; i <= n; i++){

				S[i] = 0; U[i] = D[i] = i;

				L[i + 1] = i; R[i] = i + 1;

			}

			R[n] = 0;

			memset(H, -1, sizeof(H)); size = n + 1;



			double mid = (l + r) / 2;

			for (int i = 1; i <= m; i++){

				for (int j = 1; j <= n; j++){

					if (mid>=dist(xx[i], yy[i], x[j], y[j])){

						link(i, j);

					}

				}

			}

			ak = maxn;

			dance(0);

			if (ak <= K) r = mid;

			else l = mid;

		}

		printf("%.6lf\n", l);

	}

	return 0;

}