HDU1542 Atlantis

Posted 2022-07-14 onionqaq

线段树扫描线模板

#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
#define full(a, b) memset(a, b, sizeof a)
#define __fastIn ios::sync_with_stdio(false), cin.tie(0)
#define pb push_back
using namespace std;
typedef long long LL;
inline int lowbit(int x) return x & (-x); 
inline int read()
    int ret = 0, w = 0; char ch = 0;
    while(!isdigit(ch))
        w |= ch == '-', ch = getchar();
    
    while(isdigit(ch))
        ret = (ret << 3) + (ret << 1) + (ch ^ 48);
        ch = getchar();
    
    return w ? -ret : ret;

template <typename A>
inline A __lcm(A a, A b) return a / __gcd(a, b) * b; 
template <typename A, typename B, typename C>
inline A fpow(A x, B p, C lyd)
    A ans = 1;
    for(; p; p >>= 1, x = 1LL * x * x % lyd)if(p & 1)ans = 1LL * x * ans % lyd;
    return ans;

const int N = 200;
int n, k, cnt[N<<3], cs;
double b[N<<2], sum[N<<3];
struct Seg
    double l, r, h;
    int tag;
    Seg()
    Seg(double l, double r, double h, int tag): l(l), r(r), h(h), tag(tag)
    bool operator < (const Seg &rhs) const 
        return h < rhs.h;
    
;
vector<Seg> v;

void push_up(int rt, int l, int r)
    if(cnt[rt]) sum[rt] = b[r + 1] - b[l];
    else if(l == r) sum[rt] = 0;
    else sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];


void buildTree(int rt, int l, int r)
    if(l == r)
        cnt[rt] = sum[rt] = 0;
        return;
    
    int mid = (l + r) >> 1;
    buildTree(rt << 1, l, mid);
    buildTree(rt << 1 | 1, mid + 1, r);


void insert(int rt, int l, int r, int il, int ir, int val)
    if(l == il && r == ir)
        cnt[rt] += val;
        push_up(rt, l, r);
        return;
    
    int mid = (l + r) >> 1;
    if(ir <= mid) insert(rt << 1, l, mid, il, ir, val);
    else if(il > mid) insert(rt << 1 | 1, mid + 1, r, il, ir, val);
    else insert(rt << 1, l, mid, il, mid, val), insert(rt << 1 | 1, mid + 1, r, mid + 1, ir, val);
    push_up(rt, l, r);


void init()
    k = 0;
    v.clear(), full(b, 0);
    full(sum, 0), full(cnt, 0);


int main()

    double x1, y1, x2, y2;
    while(~scanf("%d", &n) && n)
        init();
        for(int i = 1; i <= n; i ++)
            scanf("%lf%lf%lf%lf", &x1, &y1, &x2, &y2);
            v.emplace_back(x1, x2, y1, 1);
            v.emplace_back(x1, x2, y2, -1);
            b[++k] = x1, b[++k] = x2;
        
        sort(b + 1, b + k + 1);
        k = (int)(unique(b + 1, b + k + 1) - b - 1);
        sort(v.begin(), v.end());
        buildTree(1, 1, k - 1);
        double ans = 0;
        for(int i = 0; i < v.size() - 1; i ++)
            int l = (int)(lower_bound(b + 1, b + k + 1, v[i].l) - b);
            int r = (int)(lower_bound(b + 1, b + k + 1, v[i].r) - b);
            r --;
            insert(1, 1, k - 1, l, r, v[i].tag);
            ans += sum[1] * (v[i + 1].h - v[i].h);
        
        printf("Test case #%d\nTotal explored area: %.2lf\n\n", ++ cs, ans);
    
    return 0;