POJ 1151，扫描线过样例了，但不知道哪里有问题/kk

#include <iostream>
#include <algorithm>
#include <map>

using namespace std;

typedef double db;

const int MAXN = 1e6 + 10;

struct Line {
	db l, r, h;
	int mark;
	bool operator < (const Line &p) const { return h < p.h; }
} line[MAXN << 1]; 

int tot;

db num[MAXN << 1];

struct segtree {
	int l, r, cnt;
	db len;
} t[MAXN << 2];

inline 
void pushup(int p) {
	if (t[p].cnt) t[p].len = num[t[p].r + 1] - num[t[p].l];
	else t[p].len = t[p << 1].len + t[p << 1 | 1].len;
}

void build(int l, int r, int p) {
	t[p].l = l, t[p].r = r;
	if (l == r) return ;
	int mid = l + r >> 1;
	build(l, mid, p << 1), build(mid + 1, r, p << 1 | 1);
}

void add(db l, db r, int k, int p) {
	if (num[t[p].r + 1] <= l || r <= num[t[p].l]) return ;
	if (l <= num[t[p].l] && num[t[p].r + 1] <= r) return t[p].cnt += k, pushup(p);
	add(l, r, k, p << 1), add(l, r, k, p << 1 | 1);
	pushup(p);
}

int n;

db ax, ay, bx, by, ans;

int main() {
	for (int c = 1; scanf("%d", &n), n; c++) {
		tot = 0, ans = 0;
		for (int i = 1; i <= n; i++) {
			scanf("%lf%lf%lf%lf", &ax, &ay, &bx, &by);
			num[++tot] = ax, line[tot] = (Line){ ax, bx, ay, 1 };
			num[++tot] = bx, line[tot] = (Line){ ax, bx, by, -1 };
		}
		sort(line + 1, line + tot + 1);
		sort(num + 1, num + tot + 1);
		n = tot, tot = unique(num + 1, num + tot + 1) - num - 1;
		build(1, tot - 1, 1);
		for (int i = 1; i < n; i++) {
			add(line[i].l, line[i].r, line[i].mark, 1);
			ans += t[1].len * (line[i + 1].h - line[i].h);
		}
		printf("Test case #%d\nTotal explored area: %.2lf\n\n", c, ans);
	}
}