Atlantis HDU - 1542 （扫描线，线段树）

2021-10-10 10:22:58

扫描线的模板题，先把信息接收，然后排序，记录下上边和下边，然后用一条虚拟的线从下往上扫。如果我扫到的是下边，那么久用线段树在这个区间内加上1，表示这个区间现在是有的，等我扫描到上边的时候在加上-1，把之前的消掉，然后线段树维护区间内的长度，这里不是直接用下标维护，而是需要另一个数组来维护，每次记录我当前的下标在原本的图中的长度。

在update部分用一个点表示我这个点以后的一个长度为1的区间，然后这样算出来的sum[1]就是我现在线段树里包括的区间大小，然后用现在的区间大小去乘上这一部分的高度，然后每次相加，就可以了

#include<map>

#include<set>

#include<ctime>

#include<cmath>

#include<stack>

#include<queue>

#include<string>

#include<vector>

#include<cstdio>

#include<cstdlib>

#include<cstring>

#include<iostream>

#include<algorithm>

#define lowbit(x) (x & (-x))

typedef unsigned long long int ull;

typedef long long int ll;

const double pi = 4.0*atan(1.0);

const int inf = 0x3f3f3f3f;

const int maxn = ;

const int maxm = ;

using namespace std;

int n, m, tol, T;

struct Node {

    double l, r, h;

    int f;

    bool operator <(Node a) const {

        return h < a.h;

    }

};

Node node[maxn];

double a[maxn];

double sum[maxn << ];

int cnt[maxn << ];

void init() {

    memset(sum, , sizeof sum);

    memset(cnt, , sizeof cnt);

}

void pushup(int left, int right, int root) {

    if(cnt[root] != )        sum[root] = a[right + ] - a[left];

    else if(left == right)    sum[root] = ;

    else    sum[root] = sum[root << ] + sum[root <<  | ];

}

void update(int left, int right, int prel, int prer, int val, int root) {

    if(prel <= left && right <= prer) {

        cnt[root] += val;

        pushup(left, right, root);

        return ;

    }

    int mid = (left + right) >> ;

    if(prel <= mid)    update(left, mid, prel, prer, val, root << );

    if(prer > mid)    update(mid+, right, prel, prer, val, root <<  | );

    pushup(left, right, root);

}

int main() {

    int cas = ;

    while(scanf("%d", &n)==  && n) {

        init();

        double x1, y1, x2, y2;

        for(int i=; i<=n; i++) {

            scanf("%lf%lf%lf%lf", &x1, &y1, &x2, &y2);

            node[*i-].l = x1;

            node[*i].l = x1;

            node[*i-].r = x2;

            node[*i].r = x2;

            node[*i-].h = y1;

            node[*i].h = y2;

            node[*i-].f = ;

            node[*i].f = -;

            a[*i-] = x1;

            a[*i] = x2;

        }

        n = *n;

        sort(node+, node++n);

        sort(a+, a++n);

        int nn = unique(a+, a++n) - (a+);

        double ans = ;

        for(int i=; i<n; i++) {

            int l = lower_bound(a+, a++nn, node[i].l) - a;

            int r = lower_bound(a+, a++nn, node[i].r) - a;

            update(, nn, l, r-, node[i].f, );

            ans += (node[i+].h - node[i].h) * sum[];

        }

        printf("Test case #%d\n", cas++);

        printf("Total explored area: %0.2f\n", ans);

        printf("\n");

    }

    return ;

}

码农公寓

相关文章