POJ 1379 Run Away 模拟退火

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <ctime>

using namespace std;
const int maxn = 1005;
const double eps = 1e-8;
const double r = 0.99;
const int dx[] = { 0, 0, 1, -1, 1, -1, 1, -1 };
const int dy[] = { 1, -1, 0, 0, 1, -1, -1, 1 };

int px[maxn], py[maxn];
int X, Y, M;

double sq(double x) {
	return x * x;
}

double dis(double x, double y) {
	double ret = 1e90;
	for (int i = 0; i < M; i++) {
		ret = min(ret, sqrt(sq(x - px[i]) + sq(y - py[i])));
	}
	return ret;
}

double solve(double &x, double &y) {
	double nowdis = dis(x, y);
	double step = max(X, Y);
	while (step > eps) {
		for (int i = 0; i < 8; i++) {
			double nx = x + dx[i] * step,
				ny = y + dy[i] * step,
				ndis = dis(nx, ny);
			if (nx > X || nx < 0 || ny > Y || ny < 0) continue;
			if (ndis > nowdis) {
				x = nx; y = ny; nowdis = ndis;
			}
		}
		step *= r;
	}
	return nowdis;
}

int main() {
	int T; scanf("%d", &T);
	srand(time(NULL));
	while (T--) {
		scanf("%d%d%d", &X, &Y, &M);
		for (int i = 0; i < M; i++) scanf("%d%d", &px[i], &py[i]);
		double nowdis = 0, x, y;
		for (int i = 0; i < 5; i++) {
			double tx = rand() % (X), ty = rand() % (Y);
			double ret = solve(tx, ty);
			if (ret > nowdis) {
				nowdis = ret;
				x = tx; y = ty;
			}
		}
		printf("The safest point is (%.1f, %.1f).\n",
			x, y);
	}
	return 0;
}