以前对离散化的理解不够,所以把端点和区间区分来考虑但是做完这题以后有了新的认识:

先来看一个问题:给你以下的网格,你需要多少空间去存储红点区间的信息呢?

只需要图上所示的1,2,3,4个点就足够表示红点所在区间了,为什么不是一个区间的第一个红点和最后一个红点呢?(如果这样记录的话则必须加一区间点,记录区间内部信息,因为端点可能是两个区间的交集而区间内可能只被操作了一次)这样做的好处是空白区域的长度也能轻易计算出来。

因此离散化的核心在于以点代表区间

为了计算总区间两端空白的长度,增加A和B点。

原数据找离散后的值直接二分,没必要建map,因为map的查询也是logN级别的

知道了这些,存储问题就解决了。

/*
Created by Rey Chen on 2015.6.25
*/

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
//#define local
using namespace std;
const int maxn = + ;
const int max2n = + ;
const int maxc = + ;//coordinate;
typedef int ll;
// original data
int n, x0[maxn], y0[maxn], z0[maxn], x1[maxn], y1[maxn], z1[maxn];
// discretized data
int xs[max2n], ys[max2n], zs[max2n];
int xn, yn, zn;//the numbers of each coordinate //1 means solid, 2 means visited air
int color[max2n][max2n][max2n]; void discretize(int* x, int& n) {
sort(x, x+n);
n = unique(x, x+n) - x;
} struct Cell
{
int x, y, z;
Cell(int x,int y,int z):x(x),y(y),z(z){}
};
//related to flood fill
const int dx[] = {,-,,,,};
const int dy[] = {,,,-,,};
const int dz[] = {,,,,,-}; #define vis(u) (color[v.x][v.y][v.z] == 2)
#define setVis(u) color[u.x][u.y][u.z] = 2//solid..
#define uNeigh(i) Cell v(u.x+dx[i], u.y+dy[i], u.z+dz[i])
#define inside(v) (v.x >=0 && v.x < xn && v.y >= 0 && v.y < yn && v.z >= 0 && v.z < zn)
#define issolid(v) (color[v.x][v.y][v.z] == 1)
#define dis(u,x) (x##s[u.x+1]-x##s[u.x])
#define vol(u) (ll)(xs[u.x+1]-xs[u.x])*(ys[u.y+1]-ys[u.y])*(zs[u.z+1]-zs[u.z])//为了让原本最后一层空气计算时不下标不越界 还要再加一层 inline ll area(Cell& v, int dir){
if(dx[dir]) return (ll)dis(v,y)*dis(v,z);
if(dy[dir]) return (ll)dis(v,x)*dis(v,z);
return (ll)dis(v,x)*dis(v,y);
} //dfs might stack overflow
void floodfill(ll &V,ll &s)
{
queue<Cell> q;
Cell u(,,);// Cell u = Cell() 1,2
setVis(u);
q.push(u);
V = s = ;
while(!q.empty()) {
u = q.front();q.pop();
V += vol(u); //calculate air volume
for(int i = ; i < ; i++) {
uNeigh(i);
if(!inside(v)) continue;//最后一层不加
if(issolid(v)) s += area(v,i); //ensure air vis once meet solid add area
else if(!vis(v)) {
setVis(v);
q.push(v);
}
}
}
V = maxc*maxc*maxc - V;
} #define add(x) x##1[i] += x##0[i]
#define dadd(a) xs[t] = x##a[i]; ys[t] = y##a[i]; zs[t++] = z##a[i]
#define ID(x,n,val) (lower_bound(x,x+n,val) - x)
#define GetId(X,x) X##0 = ID(x##s, x##n, x##0[i]);X##1 = ID(x##s, x##n, x##1[i])
//X0 = ID(xs,xn,x0[i]);
int main()
{
#ifdef local
freopen("in.txt","r",stdin);
#endif // local int T;
scanf("%d",&T);
while(T--) {
scanf("%d",&n); int t = ;
*xs = *ys = *zs = ;//add air
xs[] = ys[] = zs[] = maxc;//最后面有两层空气
for(int i = ; i < n; i++) {
scanf("%d%d%d%d%d%d", x0+i, y0+i, z0+i, x1+i, y1+i, z1+i);
add(x);add(y);add(z);
dadd(); dadd();
}
xn = yn = zn = t ;
discretize(xs,xn);
discretize(ys,yn);
discretize(zs,zn);
xn--;yn--;zn--;//不访问最后一层空气,计算总体积时要加上最后一层空气 memset(color,,sizeof(color));
for(int i = ; i < n; i++) { //solid set
int X0, X1, Y0, Y1, Z0, Z1;
GetId(X,x); GetId(Y,y); GetId(Z,z);
for(int X = X0; X < X1; X++)
for(int Y = Y0; Y < Y1; Y++)
for(int Z = Z0; Z < Z1; Z++)
color[X][Y][Z] = ;
}
ll v,s;
floodfill(v,s);
printf("%d %d\n",s,v);
}
return ;
}

UVA 12171 (hdu 2771)sculptrue(离散化)的更多相关文章

  1. uva 12171 hdu 1771 Sculpture

    //这题从十一点开始写了四十分钟 然后查错一小时+ 要吐了 这题题意是给很多矩形的左下角(x,y,z最小的那个角)和三边的长(不是x,y,z最大的那个角T-T),为组成图形的面积与表面积(包在内部的之 ...

  2. hdu 2771(uva 12171) Sculpture bfs+离散化

    题意: 给出一些边平行于坐标轴的长方体,这些长方体可能相交.也可能相互嵌套.这些长方体形成了一个雕塑,求这个雕塑的整体积和表面积. 题解: 最easy想到直接进行bfs或者dfs统计,但此题的麻烦之处 ...

  3. Uva 12171 Sculpture - 离散化 + floodfill

    题目连接:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem ...

  4. UVa 12171 (离散化 floodfill) Sculpture

    题意: 三维空间中有n个长方体组成的雕塑,求表面积和体积. 分析: 我们可以在最外边加一圈“空气”,然后求空气的连通块的体积,最后用总体积减去即是雕塑的体积. 还有一个很“严重”的问题就是5003所占 ...

  5. UVa 221城市正视图(离散化)

    https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem& ...

  6. UVA 10816 + HDU 1839 Dijstra + 二分 (待研究)

    UVA 题意:两个绿洲之间是沙漠,沙漠的温度不同,告诉起点,终点,求使得从起点到终点的最高温度最小的路径,如果有多条,输出长度最短的路径: 思路:用最小费用(最短路径)最大流(最小温度)也能搞吧,但因 ...

  7. HDU - 1255 扫描线+离散化进阶

    这道题最开始我以为和HDU - 1542 那道题一样,只需要把cover次数改成2次即可,但是后面仔细一想,我们需要求的是覆盖次数大于等于2次的,这样的话,我们需要维护两个长度,HDU-1542 由于 ...

  8. HDU 5925 Coconuts 离散化

    题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5925 Coconuts Time Limit: 9000/4500 MS (Java/Others) ...

  9. hdu 5303 DP(离散化,环形)+贪心

    题目无法正常粘贴,地址:http://acm.hdu.edu.cn/showproblem.php?pid=5303 大意是给出一个环形公路,和它的长度,给出若干颗果树的位置以及树上的果子个数. 起点 ...

随机推荐

  1. linux下使用c判断文件夹是否为空的小程序

    自己写了一个 判断文件夹是否为空的小代码 //文件夹操作相关的函数的帮助$: man 3 readdir #include <stdio.h> #include <sys/types ...

  2. Lightoj1028【计算约数个数】

    思路: 最终就是求一个数的约数(除了1)对吧. 然后想要枚举sqrt(N)受阻,枚举素数数组受阻,加上prime[i]*prime[i]<=n就好了?那就好了吧. #include <bi ...

  3. Node-SASS安装

    鉴于国内的环境,node-sass实在是太难安装了,可以直接通过淘宝的npm镜像来安装. 1.安装cnpm(https://npm.taobao.org/) npm install -g cnpm - ...

  4. git常用命令和linux下安装git教程

    Debian/Ubuntu下使用apt-get安装git #安装libz-dev libssl-dev $ apt-get install libcurl4-gnutls-dev libexpat1- ...

  5. 黑马学习Ajax 概念和基本使用

  6. 单例设计模式singleton

    简单引入 单例设计模式作为最简单,最常用的设计模式.一般是这两中写法,这两种写法教科书所谓的标准写法,但是实际上存在不少问题.后面介绍标准写法,以规避这些问题. 1.懒汉式: /** * 问题在于,当 ...

  7. docker 在Windows下使用遇到的坑

    1.大部分系统不支持直接安装docker for windows,只能使用docker toolbox,相当于在Windows上安装了一个linux的虚拟机 2.启动docker toolbox的时候 ...

  8. Codeforces Round #375 (Div. 2) D. Lakes in Berland 并查集

    http://codeforces.com/contest/723/problem/D 这题是只能把小河填了,题目那里有写,其实如果读懂题这题是挺简单的,预处理出每一块的大小,排好序,从小到大填就行了 ...

  9. 解决win10下python3和python2共存pip的问题

    经过在查阅网友的各种经验,发现仍然解决不了问题,python2和python3在win10下的安装就不再演示了,直接在python的官网下载就好,我机器上使用的是python2.7.15和python ...

  10. Spark Mllib里如何将数据取消暂存在内存中(图文详解)

    不多说,直接上干货! 具体,见 Hadoop+Spark大数据巨量分析与机器学习整合开发实战的第18章 决策树回归分类Bike Sharing数据集