As chief engineer of the Starship Interprize, the task of repairing the hyperstellar, cubic, transwarped-out software has fallen on your shoulders. Simply put, you must compute the volume of the intersection of anywhere from 2 to 1000 cubes.
Input
The input data file consists of several sets of cubes for which the volume of their intersections must be computed. The first line of the data file contains a number (from 2 to 1000) which indicates the number of cubes which follow, one cube per line. Each line which describes a cube contains four integers. The first three integers are the x, y, and z coordinates of the corner of a cube, and the fourth integer is the positive distance which the cube extends in each of the three directions (parallel to the x, y, and z axes) from that corner.
Output
Following the data for the first set of cubes will be a number which indicates how many cubes are in a second set, followed by the cube descriptions for the second set, again one per line. Following this will be a third set, and so on. Your program should continue to process sets of cubes, outputting the volume of their intersections to the output file, one set per line, until a zero is read for the number of cubes.
Note that the data file will always contain at least one set of cubes, and every set will contain at least 2 and at most 1000 cubes. For any given set of cubes, the volume of their intersections will not exceed 1,000,000 units.
Sample Input
2 0 0 0 10 9 1 1 5 3 0 0 0 10 9 1 1 5 8 2 2 3 0
Sample Output
25 9
题意:给出三维空间的一个点以及从这个点开始延长成的长方体的相交体积
分析:通过对样例的分析不难得出,只需求出长宽高的公共长度,相乘即可;注意没有交点或只有一个交点的情况
代码:
#include<iostream>
#include<algorithm>
#include<cstring>
#include<string.h>
using namespace std;
const int inf=999999999;
int mapp[1000+10][5];
int main()
{
int n;
while(scanf("%d",&n)!=EOF)
{
if(n==0)
break;
for(int i=1;i<=n;i++)
cin>>mapp[i][1]>>mapp[i][2]>>mapp[i][3]>>mapp[i][4];
//for(int i=1;i<=n;i++)
// cout<<mapp[i][1]<<" "<<mapp[i][2]<<" "<<mapp[i][3]<<" "<<mapp[i][4];
int a1=0,a2=0,a3=0;
int b1=inf,b2=inf,b3=inf;
for(int i=1;i<=n;i++)
{
a1=max(a1,mapp[i][1]);
a2=max(a2,mapp[i][2]);
a3=max(a3,mapp[i][3]);
b1=min(b1,mapp[i][1]+mapp[i][4]);
b2=min(b2,mapp[i][2]+mapp[i][4]);
b3=min(b3,mapp[i][3]+mapp[i][4]);
}
if(b1<a1||b2<a2||b3<a3)
cout<<"0"<<endl;
else
cout<<(b1-a1)*(b2-a2)*(b3-a3)<<endl;
}
return 0;
}