アトランティス - スキャン

いくつかの古代ギリシャ語の本の中で伝説のアトランティス島の記述が含まれていますがあります。
これらの島々やマップのも、一部の一部。
しかし残念ながら、これらのマップは、アトランティスの異なる領域を記述する。
あなたの友人ビルはマップの総面積を知っている必要があります。
あなたは、総面積を計算するプログラムを書くために志願しました。
入力形式の
入力はテストケースを複数含みます。
各試験のために、最初の行は、整数n、機能の総数を含んでいます。
次のn行は、各マップが示され、各列は、4つの番号X1、Y1、X2、y2x1、含ま Y1、X2、Y2（ 必ずしも整数）、（X1、Y1）、（ X1、Y1） と（X2を、Y2）、（X2、Y2）が左上隅と地図位置の右下隅です。
座標軸x軸から下方に延び、y軸は左から右に向かって延びていることに留意されたいです。
実施形態では、n = 0を入力すると、入力終了、アクションは不要です。このユースケースを表します。
出力形式
各試験2本の出力線。
kは1つの開始からテストケースの数であり、出力「テストケース#K」の最初の行。
「エリアを探検合計：」出力の2行目、マップの総面積である（注複数のマップを含む領域であれば、矩形の全領域におけるこの試験で、すなわちとは、総面積を算出する際には一度のみカウント） 、小数点の後に2桁の精度。
空白行の各テスト出力後。
データ範囲
1≤n≤1000001≤n≤100000、

0≤x1<x2≤1000000≤x1<x2≤100000、

0≤y1<y2≤1000000≤y1<y2≤100000

この質問の上限は十億万に強化nnはあることに注意してください。
サンプル入力：
2
10 10 20 20である
15 15 25 25.5
0

サンプル出力：
テストケース＃1
総探検面積：180.00
あまりにもハード、アイデアを書き、確かに混乱の書き込みはありません
、次のようにACコードは次のとおりです。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;
const int N = 100010;
int n;
struct Segment
{
    double x, y1, y2;
    int k;
    bool operator< (const Segment &t)const
    {
        return x < t.x;
    }
}seg[N * 2];
struct Node
{
    int l, r;
    int cnt;
    double len;
}tr[N * 8];
vector<double> ys;
int find(double y)
{
    return lower_bound(ys.begin(), ys.end(), y) - ys.begin();
}
void pushup(int u)
{
    if (tr[u].cnt) tr[u].len = ys[tr[u].r + 1] - ys[tr[u].l];
    else if (tr[u].l != tr[u].r)
    {
        tr[u].len = tr[u << 1].len + tr[u << 1 | 1].len;
    }
    else tr[u].len = 0;
}
void build(int u, int l, int r)
{
    tr[u] = {l, r, 0, 0};
    if (l != r)
    {
        int mid = l + r >> 1;
        build(u << 1, l, mid), build(u << 1 | 1, mid + 1, r);
    }
}
void modify(int u, int l, int r, int k)
{
    if (tr[u].l >= l && tr[u].r <= r)
    {
        tr[u].cnt += k;
        pushup(u);
    }
    else
    {
        int mid = tr[u].l + tr[u].r >> 1;
        if (l <= mid) modify(u << 1, l, r, k);
        if (r > mid) modify(u << 1 | 1, l, r, k);
        pushup(u);
    }
}
int main()
{
    int T = 1;
    while (scanf("%d", &n), n)
    {
        ys.clear();
        for (int i = 0, j = 0; i < n; i ++ )
        {
            double x1, y1, x2, y2;
            scanf("%lf%lf%lf%lf", &x1, &y1, &x2, &y2);
            seg[j ++ ] = {x1, y1, y2, 1};
            seg[j ++ ] = {x2, y1, y2, -1};
            ys.push_back(y1), ys.push_back(y2);
        }
sort(ys.begin(), ys.end());
        ys.erase(unique(ys.begin(), ys.end()), ys.end());
build(1 ,0, ys.size() - 2);
sort(seg, seg + n * 2);
 double res = 0;
        for (int i = 0; i < n * 2; i ++ )
        {
            if (i > 0) res += tr[1].len * (seg[i].x - seg[i - 1].x);
            modify(1, find(seg[i].y1), find(seg[i].y2) - 1, seg[i].k);
        }   
    printf("Test case #%d\n", T ++ );
        printf("Total explored area: %.2lf\n\n", res);
    }
 return 0;
}