hdu1542线段树 扫描线 求面积并

各变量名称严格按照算法竞赛进阶指南的起名习惯,可结合起来研究。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <cmath>

using namespace std;
#define MID ((l(p) + r(p)) >> 1)

const int N = 1e5 + 10;

struct Side {
    double x, y1, y2;
    int k;
    Side(double a, double b, double c, int d) : x(a), y1(b), y2(c), k(d) {}
    Side() {}
} side[N];

struct Node {
    int l, r, cnt;
    double len;
#define l(x) c[x].l
#define r(x) c[x].r
#define len(x) c[x].len
#define cnt(x) c[x].cnt
} c[N];

double raw[N];
int ccnt = 0;

bool cmp(Side a, Side b) {
    return a.x < b.x;
}

bool cmp1(double a, double b) {
    return (fabs(a - b) < 1e-8);
}

int val(double d) {
    return lower_bound(raw, raw + ccnt, d) - raw;
}

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

void build(int p, int l, int r) {
    l(p) = l, r(p) = r, cnt(p) = 0, len(p) = 0.0;
    if (l == r) return;
    build(p << 1, l, MID);
    build(p << 1 | 1, MID + 1, r);
}

void change(int p, int l, int r, int k) {
    if (l <= l(p) && r(p) <= r) {
        cnt(p) += k, pushup(p);
        return;
    }
    if (l <= MID) change(p << 1, l, r, k);
    if (r > MID) change(p << 1 | 1, l, r, k);
    pushup(p);
}

double solve(int n) {
    ccnt = 0;
    for (int i = 0; i < n; i++) {
        double x1, x2, y1, y2;
        scanf("%lf%lf%lf%lf", &x1, &y1, &x2, &y2);
        raw[ccnt++] = y1, raw[ccnt++] = y2;
        side[i] = Side(x1, y1, y2, 1);
        side[i + n] = Side(x2, y1, y2, -1);
    }
    sort(raw, raw + ccnt);
    unique(raw, raw + ccnt, cmp1);
    sort(side, side + 2 * n, cmp);
    build(1, 0, ccnt);
    double ans = 0.0;
    for (int i = 0; i < 2 * n; i++) {
        change(1, val(side[i].y1), val(side[i].y2) - 1, side[i].k);
        ans += len(1) * (side[i + 1].x - side[i].x);
    }
    return ans;
}

int main() {
    int t, cas = 0;
    while (cin >> t, t) {
        printf("Test case #%d\nTotal explored area: %.2f\n\n", ++cas, solve(t));
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值