码迷,mamicode.com
首页 > 其他好文 > 详细

HDU 1542 Atlantis(矩形面积并)

时间:2014-10-14 14:06:38      阅读:119      评论:0      收藏:0      [点我收藏+]

标签:style   blog   http   color   io   os   ar   for   sp   

HDU 1542 Atlantis

题目链接

题意:给定一些矩形,求面积并

思路:利用扫描线,由于这题矩形个数不多,直接暴力扫就可以了,如果数据大,就要用线段树

代码:

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;

const int N = 205;
const int M = 100005;
const double eps = 1e-8;

int n, vis[N], hn;
double hash[N];

struct Line {
	double l, r, y;
	int flag;
	Line() {}
	Line(double l, double r, double y, int flag) {
		this->l = l;
		this->r = r;
		this->y = y;
		this->flag = flag;
	}
} line[N];

bool cmp(Line a, Line b) {
	return a.y < b.y;
}

int get(double x) {
	return lower_bound(hash, hash + hn, x) - hash;
}

int main() {
	int cas = 0;
	while (~scanf("%d", &n) && n) {
		double x1, x2, y1, y2; hn = 0;
		memset(vis, 0, sizeof(vis));
		for (int i = 0; i < n; i++) {
			scanf("%lf%lf%lf%lf", &x1, &y1, &x2, &y2);
			line[i * 2] = Line(x1, x2, y1, 1);
			line[i * 2 + 1] = Line(x1, x2, y2, -1);
			hash[hn++] = x1; hash[hn++] = x2;
		}
		n *= 2;
		sort(line, line + n, cmp);
		sort(hash, hash + hn);
		hn = 1;
		for (int i = 1; i < n; i++) {
			if (fabs(hash[i] - hash[i - 1]) < eps) continue;
			hash[hn++] = hash[i];
		}
		double ans = 0;
		for (int i = 0; i < n; i++) {
			int l = get(line[i].l), r = get(line[i].r);
			double len = 0;
			for (int j = 0; j < hn - 1; j++) if (vis[j] > 0) len += (hash[j + 1] - hash[j]);
			if (i) ans += len * (line[i].y - line[i - 1].y);
			for (int j = l; j < r; j++) vis[j] += line[i].flag;
		}
		printf("Test case #%d\n", ++cas);
		printf("Total explored area: %.2lf\n\n", ans);
	}
	return 0;
}


HDU 1542 Atlantis(矩形面积并)

标签:style   blog   http   color   io   os   ar   for   sp   

原文地址:http://blog.csdn.net/accelerator_/article/details/40074903

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!