uva 12171 sculpture （超级好题）——yhx
1 #include<cstdio> 2 #include<algorithm> 3 #include<map> 4 #include<cstring> 5 #include<queue> 6 using namespace std; 7 int a[110][110][110],x_map[110],y_map[110],z_map[110], 8 x1_ori[110],y1_ori[60],z1_ori[60],x2_ori[60],y2_ori[60],z2_ori[60], 9 x_move[6]={1,-1,0,0,0,0},y_move[6]={0,0,1,-1,0,0},z_move[6]={0,0,0,0,1,-1}, 10 x_cnt,y_cnt,z_cnt; 11 long long ans_s,ans_v; 12 map<int,int> x_rev,y_rev,z_rev; 13 struct cond 14 { 15 int x,y,z; 16 }c1,c2; 17 queue<cond> q; 18 int main() 19 { 20 int i,j,k,l,m,n,p,x,y,z,t,xx,yy,zz,temp; 21 scanf("%d",&t); 22 while (t--) 23 { 24 memset(a,0,sizeof(a)); 25 memset(x_map,0,sizeof(x_map)); 26 memset(y_map,0,sizeof(y_map)); 27 memset(z_map,0,sizeof(z_map)); 28 memset(x1_ori,0,sizeof(x1_ori)); 29 memset(y1_ori,0,sizeof(y1_ori)); 30 memset(z1_ori,0,sizeof(z1_ori)); 31 memset(x2_ori,0,sizeof(x2_ori)); 32 memset(y2_ori,0,sizeof(y2_ori)); 33 memset(z2_ori,0,sizeof(z2_ori)); 34 x_rev.clear(); 35 y_rev.clear(); 36 z_rev.clear(); 37 scanf("%d",&n); 38 for (i=1;i<=n;i++) 39 { 40 scanf("%d%d%d%d%d%d",&x1_ori[i],&y1_ori[i],&z1_ori[i],&x,&y,&z); 41 x_map[i]=x1_ori[i]; 42 y_map[i]=y1_ori[i]; 43 z_map[i]=z1_ori[i]; 44 x_map[i+n]=x2_ori[i]=x1_ori[i]+x; 45 y_map[i+n]=y2_ori[i]=y1_ori[i]+y; 46 z_map[i+n]=z2_ori[i]=z1_ori[i]+z; 47 } 48 sort(x_map+1,x_map+2*n+1); 49 sort(y_map+1,y_map+2*n+1); 50 sort(z_map+1,z_map+2*n+1); 51 x_cnt=unique(x_map+1,x_map+2*n+1)-(x_map+1); 52 y_cnt=unique(y_map+1,y_map+2*n+1)-(y_map+1); 53 z_cnt=unique(z_map+1,z_map+2*n+1)-(z_map+1); 54 for (i=1;i<=x_cnt;i++) 55 x_rev[x_map[i]]=i; 56 for (i=1;i<=y_cnt;i++) 57 y_rev[y_map[i]]=i; 58 for (i=1;i<=z_cnt;i++) 59 z_rev[z_map[i]]=i; //以上为离散化 60 x_map[0]=y_map[0]=z_map[0]=0; 61 x_map[x_cnt+1]=y_map[y_cnt+1]=z_map[z_cnt+1]=1005; //边界，要在外面多围一圈空气。 62 for (i=1;i<=n;i++) 63 for (x=x_rev[x1_ori[i]];x<x_rev[x2_ori[i]];x++) 64 for (y=y_rev[y1_ori[i]];y<y_rev[y2_ori[i]];y++) 65 for (z=z_rev[z1_ori[i]];z<z_rev[z2_ori[i]];z++) 66 a[x][y][z]=1; //填充 67 ans_v=1005*1005*1005; 68 ans_s=0; 69 c1.x=c1.y=c1.z=0; 70 q.push(c1); 71 a[0][0][0]=2; //2表示已经入过队 72 while (!q.empty()) 73 { 74 c1=q.front(); 75 q.pop(); 76 xx=c1.x; 77 yy=c1.y; 78 zz=c1.z; 79 ans_v-=(x_map[xx+1]-x_map[xx])*(y_map[yy+1]-y_map[yy])*(z_map[zz+1]-z_map[zz]); 80 for (i=0;i<=5;i++) 81 { 82 x=xx+x_move[i]; 83 y=yy+y_move[i]; 84 z=zz+z_move[i]; 85 if (x>=0&&x<=x_cnt&&y>=0&&y<=y_cnt&&z>=0&&z<=z_cnt) 86 { 87 if (a[x][y][z]==0) 88 { 89 c2.x=x; 90 c2.y=y; 91 c2.z=z; 92 a[x][y][z]=2; //要在入队时标记，不能在取出时。 93 q.push(c2); 94 } 95 if (a[x][y][z]==1) 96 { 97 temp=1; 98 if (x==xx) temp*=(x_map[xx+1]-x_map[xx]); 99 if (y==yy) temp*=(y_map[yy+1]-y_map[yy]); 100 if (z==zz) temp*=(z_map[zz+1]-z_map[zz]); 101 ans_s+=temp; //算面积 102 } 103 } 104 } 105 } 106 printf("%lld %lld\n",ans_s,ans_v); 107 } 108 }
离散化+floodfill。
最开始用的dfs结果RE想到可能是堆栈溢出，改用bfs后AC。
由于体积和面积都是从外面看的，可以从外面一圈“空气”开始floodfill，每次遇到雕塑就加上表面积，因为每一块表面积都会且只会和一块空气接触。体积用总体积减去遇到的空气体积即可。
由于坐标较大，需要离散化。
