題意：

給一些長(zhǎng)方體，求這些長(zhǎng)方體相交至少3次的體積的并。

思路：

1. 先注意到 z 的變化范圍 0 - 500，于是可以先把 z 離散化再枚舉 z 坐標(biāo)的高度。

2. 對(duì)于每一個(gè) z[i] - z[i+1] 的高度，找出跨越這個(gè)長(zhǎng)的立方體，然后根據(jù)面積的方法求出掃描線 >= 3 的個(gè)數(shù)。

3. 因?yàn)轭}目最終要求至少相交 3 次體積的并，對(duì)于每個(gè)固定 z 區(qū)間，要求 x-y 平面面積重合度 >= 3 即可，設(shè)置了幾組標(biāo)記：

? ?a. cnt[] 表示整個(gè)區(qū)間已經(jīng)被掃描線覆蓋的次數(shù)： cnt[rt] = 1 表示區(qū)間被掃描線完整的覆蓋了 1 次。 等于 2 和等于 3 都是一樣的道理。

? ?b. once[] 表示區(qū)間內(nèi)僅僅被 1 次覆蓋的掃描線的長(zhǎng)度，twice[] 表示區(qū)間內(nèi)僅僅被 2 次覆蓋的掃描線的長(zhǎng)度。

? ?c. sum[] 表示區(qū)間最終被覆蓋 >= 3 的掃描線長(zhǎng)度，對(duì)于不同的 cnt[]，sum[] once[] twice[] 都有不同的計(jì)算方法，具體見(jiàn)代碼注釋。

?

#include <iostream> #include <algorithm> using namespace std;#define lhs l, m, rt << 1 #define rhs m + 1, r, rt << 1 | 1const int maxn = 2222; int sum[maxn << 2], cnt[maxn << 2]; int once[maxn << 2], twice[maxn << 2]; int xcord[maxn], zcord[maxn];struct Cube {int x1, y1, z1;int x2, y2, z2; } cube[maxn] ;struct Segment {int l, r, h, v;Segment() { } Segment(int _l, int _r, int _h, int _v): l(_l), r(_r), h(_h), v(_v) { }bool operator < (const Segment& other){if (h == other.h)return v > other.v;elsereturn h < other.h;} } seg[maxn] ; void pushUp(int l, int r, int rt) {if (cnt[rt] >= 3){sum[rt] = xcord[r + 1] - xcord[l];　　　　　　　　// 表示整個(gè)區(qū)間已經(jīng)被掃描線覆蓋 3 次，sum[] 自然為區(qū)間長(zhǎng)度once[rt] = twice[rt] = 0;　　　　　　　　　　　　　// once[], twice[] 都為 0 表示區(qū)間內(nèi)沒(méi)有僅覆蓋 1 次或 2 次的掃描線}else if (cnt[rt] == 2){sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];　　　　sum[rt] += once[rt << 1] + once[rt << 1 | 1];　　// 因?yàn)檎麄€(gè)區(qū)間已經(jīng)被掃描線覆蓋 2 次，所以 once[] 一定要為 0 sum[rt] += twice[rt << 1] + twice[rt << 1 | 1]; // 由于線段樹(shù)中 cnt[] 并不是更新到底的，所以 sum[] 要加上左右子區(qū)間 once[], twice[] 的值once[rt] = 0;twice[rt] = xcord[r + 1] - xcord[l] - sum[rt];}else if (cnt[rt] == 1){sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];sum[rt] += twice[rt << 1] + twice[rt << 1 | 1];twice[rt] = once[rt << 1] + once[rt << 1 | 1];once[rt] = xcord[r + 1] - xcord[l] - sum[rt] - twice[rt];}else if (cnt[rt] == 0){sum[rt] = sum[rt << 1] + sum[rt << 1 | 1];once[rt] = once[rt << 1] + once[rt << 1 | 1];twice[rt] = twice[rt << 1] + twice[rt << 1 | 1];} }void update(int beg, int end, int value, int l, int r, int rt) {if (beg <= l && r <= end){cnt[rt] += value;pushUp(l, r, rt);return ;}int m = (l + r) >> 1;if (beg <= m)update(beg, end, value, lhs);if (end > m)update(beg, end, value, rhs);pushUp(l, r, rt); }__int64 solve(int n) {__int64 ret = 0;int m = 0;for (int i = 0; i < n; ++i){xcord[m] = cube[i].x1;zcord[m++] = cube[i].z1;xcord[m] = cube[i].x2;zcord[m++] = cube[i].z2; }sort(xcord, xcord + m);sort(zcord, zcord + m);int xnum = unique(xcord, xcord + m) - xcord;int znum = unique(zcord, zcord + m) - zcord;for (int z = 0; z < znum - 1; ++z){m = 0;for (int i = 0; i < n; ++i){if (cube[i].z1 <= zcord[z] && zcord[z] < cube[i].z2){seg[m++] = Segment(cube[i].x1, cube[i].x2, cube[i].y1, 1);seg[m++] = Segment(cube[i].x1, cube[i].x2, cube[i].y2, -1);}}sort(seg, seg + m);memset(sum, 0, sizeof(sum));memset(cnt, 0, sizeof(cnt));memset(once, 0, sizeof(once));memset(twice, 0, sizeof(twice));__int64 area = 0;for (int i = 0; i < m - 1; ++i){int beg = lower_bound(xcord, xcord + xnum, seg[i].l) - xcord;int end = lower_bound(xcord, xcord + xnum, seg[i].r) - xcord;if (beg < end)update(beg, end - 1, seg[i].v, 0, xnum - 1, 1);area += (__int64)sum[1] * (seg[i+1].h - seg[i].h);}ret += area * (__int64)(zcord[z+1] - zcord[z]);}return ret; }int main() {int n, cases, cc = 0;scanf("%d", &cases);while (cases--){scanf("%d", &n);for (int i = 0; i < n; ++i){scanf("%d %d %d", &cube[i].x1, &cube[i].y1, &cube[i].z1);scanf("%d %d %d", &cube[i].x2, &cube[i].y2, &cube[i].z2);}__int64 ret;if (n < 3)ret = 0;elseret = solve(n);printf("Case %d: %I64d\n", ++cc, ret);}return 0; }

轉(zhuǎn)載于:https://www.cnblogs.com/kedebug/archive/2013/02/02/2890373.html

總結(jié)

以上是生活随笔為你收集整理的HDOJ 3642 Get The Treasury（扫描线 + 线段树 + 离散化 立方体的并）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。