Atlantis 【POJ - 1151】【扫描线模板题+线段树更新】

题目链接

---

  是一道扫描线的模板题，也是我的第一道扫描线的题了，对扫描线也算是有了第一次的理解，无非就是更新新的向上的区间长度，然后去查询就是了，而查询是O(1)的，因为可以通过树的最上根节点得到的。

---

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
const int maxN = 500;
int N, num, flag[maxN<<1], shu;
double sum[maxN<<1], X[maxN<<1];
struct node
{
    double lx, rx, y;
    int s;
    node(double a=0, double b=0, double c=0, int d=0):lx(a), rx(b), y(c), s(d) {}
}line[maxN<<1];
bool cmp(node e1, node e2) { return e1.y < e2.y; }
inline void pushup(int rt, int l, int r)    //pushup()多用，可以直接赋值时候使用
{
    if(flag[rt]) sum[rt] = X[r+1] - X[l];   //如果“>0”就说明这个点是有边的
    else if(l == r) sum[rt] = 0;    //清零操作使用
    else sum[rt] = sum[rt<<1] + sum[rt<<1|1];   //求其下的边的总和
}
inline void update(int rt, int l, int r, int ql, int qr, int val)
{
    if(ql<=l && qr>=r)
    {
        flag[rt] += val;
        pushup(myself);
        return;
    }
    int mid = HalF;
    if(ql > mid) update(QR, val);
    else if(qr <= mid) update(QL, val);
    else { update(QL, val); update(QR, val); }
    pushup(myself);
}
inline void init()
{
    num = 0;    shu = 1;
    memset(flag, 0, sizeof(flag));
    memset(sum, 0, sizeof(sum));
}
int main()
{
    int Cas = 0;
    while(scanf("%d", &N) && N)
    {
        init();
        for(int i=1; i<=N; i++)
        {
            double lx, ly, rx, ry;
            scanf("%lf%lf%lf%lf", &lx, &ly, &rx, &ry);
            line[++num] = node(lx, rx, ly, 1);
            X[num] = lx;
            line[++num] = node(lx, rx, ry, -1);
            X[num] = rx;
        }
        sort(X+1, X+1+num);
        sort(line+1, line+num+1, cmp);
        for(int i=2; i<=num; i++) if(X[i] != X[i-1]) X[++shu] = X[i];
        double ans = 0;
        for(int i=1; i<num; i++)    //离散化更新的过程
        {
            int l = (int)(lower_bound(X+1, X+shu+1, line[i].lx) - X);
            int r = (int)(lower_bound(X+1, X+shu+1, line[i].rx) - X - 1);
            update(1, 1, shu, l, r, line[i].s);
            ans += sum[1] * (line[i+1].y - line[i].y);
        }
        printf("Test case #%d\n", ++Cas);
        printf("Total explored area: %.2lf\n\n", ans);
    }
    return 0;
}
/*
2
15 15 25 25.5
10 10 20 20
ans:180
 
5
1.37 0.10 2.75 0.20
2.96 2.01 3.06 2.75
2.65 1.05 4.97 1.15
2.33 0.74 4.02 0.84
2.96 0.10 4.02 1.47
ans:1.85

8
9.33 19.52 30.44 24.39
8.06 29.38 9.75 51.77
22.38 15.06 32.35 33.30
30.66 19.20 44.13 19.62
17.93 15.38 39.35 41.58
9.65 6.47 29.80 10.39
26.20 24.61 40.63 35.85
25.25 12.51 29.80 36.17
ans:751.10
*/