HDOJ 3642 Get The Treasury(扫描线 + 线段树 + 离散化 立方体的并)
题意:
给一些长方体,求这些长方体相交至少3次的体积的并。
思路:
1. 先注意到 z 的变化范围 0 - 500,于是可以先把 z 离散化再枚举 z 坐标的高度。
2. 对于每一个 z[i] - z[i+1] 的高度,找出跨越这个长的立方体,然后根据面积的方法求出扫描线 >= 3 的个数。
3. 因为题目最终要求至少相交 3 次体积的并,对于每个固定 z 区间,要求 x-y 平面面积重合度 >= 3 即可,设置了几组标记:
a. cnt[] 表示整个区间已经被扫描线覆盖的次数: cnt[rt] = 1 表示区间被扫描线完整的覆盖了 1 次。 等于 2 和等于 3 都是一样的道理。
b. once[] 表示区间内仅仅被 1 次覆盖的扫描线的长度,twice[] 表示区间内仅仅被 2 次覆盖的扫描线的长度。
c. sum[] 表示区间最终被覆盖 >= 3 的扫描线长度,对于不同的 cnt[],sum[] once[] twice[] 都有不同的计算方法,具体见代码注释。
#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]; // 表示整个区间已经被扫描线覆盖 3 次,sum[] 自然为区间长度once[rt] = twice[rt] = 0; // once[], twice[] 都为 0 表示区间内没有仅覆盖 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]; // 因为整个区间已经被扫描线覆盖 2 次,所以 once[] 一定要为 0 sum[rt] += twice[rt << 1] + twice[rt << 1 | 1]; // 由于线段树中 cnt[] 并不是更新到底的,所以 sum[] 要加上左右子区间 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; }
转载于:https://www.cnblogs.com/kedebug/archive/2013/02/02/2890373.html
HDOJ 3642 Get The Treasury(扫描线 + 线段树 + 离散化 立方体的并)相关推荐
- 【POJ 2482】 Stars in Your Window(线段树+离散化+扫描线)
[POJ 2482] Stars in Your Window(线段树+离散化+扫描线) Time Limit: 1000MS Memory Limit: 65536K Total Submiss ...
- 【HDU5091】Beam Cannon,扫描线+线段树
传送门 思路: 扫描线的经典问题 然而并不是很会做-- 对x坐标差分一下(例如(x,y)类型为1,(x+w,y)的类型就是-1,以便于之后扫描可以去除该点),然后对y坐标转化为正数(便于操作为线段树上 ...
- 【BZOJ3958】[WF2011]Mummy Madness 二分+扫描线+线段树
[BZOJ3958][WF2011]Mummy Madness Description 在2011年ACM-ICPC World Finals上的一次游览中,你碰到了一个埃及古墓. 不幸的是,你打开了 ...
- 扫描线+线段树简介 AcWing 248窗内的星星题解
----出自南昌理工学院ACM集训队 这周学习了线段树和扫描线的解题方法,下面由小菜鸡简介一下: 一般扫描线的题目最简单的便是扫描线裸模板(一般来说的话:数据范围小),其次的话便是进行拓展成线段树+扫 ...
- HDOJ 2492 Ping pong 线段树+离散化
//2492 Ping pong 线段树+离散化 /* 题意: 有一陀人从左到右排成一排,每个人有一个唯一的技能值,每个人都找其他人比赛, 比赛前要再找一个人做裁判,裁判的技能值不能比这两个人都高,也 ...
- poj2528贴海报(线段树离散化)
//poj2528贴海报(线段树离散化) #include<cstring> #include<iostream> #include<cstdio> #includ ...
- poj 2528 Mayor's posters(线段树+离散化)
1 /* 2 poj 2528 Mayor's posters 3 线段树 + 离散化 4 5 离散化的理解: 6 给你一系列的正整数, 例如 1, 4 , 100, 1000000000, 如果利用 ...
- POJ 1177 Picture [离散化+扫描线+线段树]
http://poj.org/problem?id=1177 给若干矩形,求被覆盖的区域的周长. 将 y 坐标离散化后,按 x 坐标进行扫描.用线段树维护两个东西,当前竖线的叠加长度 len 和 条数 ...
- hdu 4419 Colourful Rectangle (离散化扫描线线段树)
Problem - 4419 题意不难,红绿蓝三种颜色覆盖在平面上,不同颜色的区域相交会产生新的颜色,求每一种颜色的面积大小. 比较明显,这题要从矩形面积并的方向出发.如果做过矩形面积并的题,用线段树 ...
最新文章
- Java项目命名规范
- 开放源码,华为鸿蒙HarmonyOS 2.0来了
- 李宏毅强化学习完整笔记!开源项目《LeeDeepRL-Notes》发布
- python map lambda 分割字符串_Python特殊语法:filter、map、reduce、lambda [转]
- 会计的思考(17):还原会计报表的企业个性之一
- ElementUI:项目中引入自己的ICON
- 并查集(ACWING三道题)
- 给linux默认mysql设置root密码
- 用C语言实现俄罗斯方块游戏
- miniprograme-wsx
- c语言string函数的用法_C语言让电脑关机?system函数功能够大够硬
- android模拟器 分辨率,Android模拟器各个皮肤的分辨率
- 多目标、多阶段、多层次的强化学习合作方法
- 哪有什么太迟,世界多的是大器晚成
- 台式计算机连wifi,台式电脑怎么连wifi
- 关于java 操作word的几种方式
- 致敬不凡·最美的星火:国产飞腾CPU研发力量
- java计算机毕业设计技术的游戏交易平台源代码+数据库+系统+lw文档
- web前段网图分类规划
- (转载)终极解密---房价,物价飞涨--分析得不错的文章
热门文章
- JQuery 的each方法
- [导入]Don's Loose
- Selenium常用方法及函数
- Android 自带图标库 android.R.drawable
- oracle 里面定时执行任务,比如存储过程内容等
- Lazarus IOCP 移植
- MYSQL5.5 YUM更新安装
- mysql 5.6 初始化_MySQL 5.6 关于登陆的初始化设置
- python数据整理代码_熬夜整理的资料:分享Python数据可视化图表代码和案例给大家...
- python中ndarray如何生成单独的记事本文件_NumPy-快速处理数据--ndarray对象--数组的创建和存取...