//这题从十一点开始写了四十分钟 然后查错一小时+ 要吐了
这题题意是给很多矩形的左下角（x,y,z最小的那个角）和三边的长（不是x,y,z最大的那个角T-T）,为组成图形的面积与表面积（包在内部的之算体积不算表面积）
解法：离散化+bfs，先把范围扩大（相当于在周围加上空气），然后bfs，遇到表面积直接加入，遇到非长方体的部分也直接加入，最后用总体积减去空气的体积，这样就可以把内部的体积计算进来而不计算其表面积。因为坐标范围比较大，要先离散化。
//其实我对这题一直耿耿于怀，当年没进省队多少与这题有关
//昨天没出题 明天周五课少 出两题
1 #include<cstdio> 2 #include<iostream> 3 #include<cmath> 4 #include<algorithm> 5 #include<cstring> 6 #include<cstdlib> 7 #include<queue> 8 #include<vector> 9 #include<map> 10 #include<stack> 11 #include<string> 12 13 using namespace std; 14 15 const int u[6]={0,0,0,0,1,-1}; 16 const int v[6]={1,-1,0,0,0,0}; 17 const int w[6]={0,0,1,-1,0,0}; 18 19 struct rec{ 20 int x1,y1,z1,x2,y2,z2; 21 }; 22 23 struct POINT{ 24 int x,y,z; 25 }; 26 27 bool vis[103][103][103]; 28 bool loc[103][103][103]; 29 int n,T,ans1,ans2; 30 int x[103],y[103],z[103]; 31 rec p[103]; 32 int numx,numy,numz; 33 34 int IDX(int aim){ 35 return lower_bound(x,x+numx,aim)-x; 36 } 37 38 int IDY(int aim){ 39 return lower_bound(y,y+numy,aim)-y; 40 } 41 42 int IDZ(int aim){ 43 return lower_bound(z,z+numz,aim)-z; 44 } 45 46 int solve(int flag,int nowx,int nowy,int nowz){//计算表面积 47 if (flag==0) return (x[nowx+1]-x[nowx])*(z[nowz+1]-z[nowz]); 48 if (flag==1) return (x[nowx+1]-x[nowx])*(z[nowz+1]-z[nowz]); 49 if (flag==2) return (x[nowx+1]-x[nowx])*(y[nowy+1]-y[nowy]); 50 if (flag==3) return (x[nowx+1]-x[nowx])*(y[nowy+1]-y[nowy]); 51 if (flag==4) return (z[nowz+1]-z[nowz])*(y[nowy+1]-y[nowy]); 52 if (flag==5) return (z[nowz+1]-z[nowz])*(y[nowy+1]-y[nowy]); 53 return -1; 54 } 55 56 void bfs(){ 57 queue<POINT> q; 58 while (!q.empty()) q.pop(); 59 q.push((POINT){0,0,0}); 60 ans2=(x[1]-x[0])*(y[1]-y[0])*(z[1]-z[0]); 61 vis[0][0][0]=1; 62 while (!q.empty()){ 63 POINT now=q.front(); 64 q.pop(); 65 for (int i=0;i<6;i++){ 66 int tx=now.x+u[i]; 67 int ty=now.y+v[i]; 68 int tz=now.z+w[i]; 69 if (tx<0 || tx>=numx-1 || ty<0 || ty>=numy-1 || tz<0 || tz>=numz-1 || vis[tx][ty][tz]) continue; 70 if (loc[tx][ty][tz]){ 71 ans1+=solve(i,now.x,now.y,now.z); 72 } 73 else{ 74 ans2+=(x[tx+1]-x[tx])*(y[ty+1]-y[ty])*(z[tz+1]-z[tz]); 75 vis[tx][ty][tz]=1; 76 q.push((POINT){tx,ty,tz}); 77 } 78 } 79 } 80 } 81 82 int main(){ 83 scanf("%d",&T); 84 for (int cas=1;cas<=T;cas++){ 85 scanf("%d",&n); 86 for (int i=0;i<n;i++){ 87 scanf("%d%d%d%d%d%d",&p[i].x1,&p[i].y1,&p[i].z1,&p[i].x2,&p[i].y2,&p[i].z2); 88 p[i].x2+=p[i].x1; 89 p[i].y2+=p[i].y1; 90 p[i].z2+=p[i].z1; 91 x[2*i+1]=p[i].x1; 92 x[2*i+2]=p[i].x2; 93 y[2*i+1]=p[i].y1; 94 y[2*i+2]=p[i].y2; 95 z[2*i+1]=p[i].z1; 96 z[2*i+2]=p[i].z2;//先把坐算出来 97 } 98 x[0]=0; 99 y[0]=0; 100 z[0]=0; 101 x[2*n+1]=1001; 102 y[2*n+1]=1001; 103 z[2*n+1]=1001;//拓展范围 104 sort(x,x+2*n+2); 105 sort(y,y+2*n+2); 106 sort(z,z+2*n+2); 107 numx=unique(x,x+2*n+2)-x; 108 numy=unique(y,y+2*n+2)-y; 109 numz=unique(z,z+2*n+2)-z;//离散化 110 memset(loc,0,sizeof(loc)); 111 memset(vis,0,sizeof(vis)); 112 for (int now=0;now<n;now++){ 113 for (int i=IDX(p[now].x1);i<IDX(p[now].x2);i++){ 114 for (int j=IDY(p[now].y1);j<IDY(p[now].y2);j++){ 115 for (int k=IDZ(p[now].z1);k<IDZ(p[now].z2);k++){ 116 loc[i][j][k]=1;//记录矩形位置 117 } 118 } 119 } 120 } 121 ans1=0;//^2 122 ans2=0;//^3 123 bfs(); 124 ans2=x[numx-1]*y[numy-1]*z[numz-1]-ans2;//总体积减去空气体积 125 printf("%d %d ",ans1,ans2); 126 } 127 return 0; 128 } 129 /* 130 1 131 2 132 1 2 3 3 4 5 133 6 2 3 3 4 5 134 135 1 136 7 137 1 1 1 5 5 1 138 1 1 10 5 5 1 139 1 1 2 1 4 8 140 2 1 2 4 1 8 141 5 2 2 1 4 8 142 1 5 2 4 1 8 143 3 3 4 1 1 1 144 145 2 146 2 147 1 2 3 3 4 5 148 6 2 3 3 4 5 149 7 150 1 1 1 5 5 1 151 1 1 10 5 5 1 152 1 1 2 1 4 8 153 2 1 2 4 1 8 154 5 2 2 1 4 8 155 1 5 2 4 1 8 156 3 3 4 1 1 1 157 */