51nod 雷达覆盖

现在有N个雷达。每个雷达的覆盖区域是一个正方体，区域用一个顶点座标P(x,y,z)和边长L表示，P是区域的一个顶点，整个区域是由P沿x,y,z轴正方向延展L得到的一个正方体。现在想知道哪些区域是被所有雷达覆盖到的，计算出这些区域的总体积。

 收起

输入

单组测试数据。
第一行有一个整数N(2<=N<=1000)，表示雷达的数目。
接下来有N行，每一行给出四个整数x, y, z （|x|, |y|, |z|<=10000 ）和L(1<=L<=10000)，表示一个雷达所覆盖区域的参数。

输出

输出一个整数表示答案。

输入样例

2
0 0 0 10
9 1 1 5

输出样例

25

这个题没有什么大的难点，琢磨了一下午好歹卡数据卡过了，刚开始没有用long long导致数据输出不出来，总而言之很不容易。

思路很简单找x，y，z的公共区域，直接把长度相乘就好，若得数为负则没有相交。

附代码

#include <cmath>
#include <cstdio>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
long long int f;
int main() {
long long int a,b,c,d,e,g1=-10001,h1=100001,x[2000]={0},y[2000]={0},z[2000]={0},l[2000];
long long int g2=-10001,h2=100001,g3=-10001,h3=100001;
    cin>>a;
    for(b=1;b<=a;++b)
     {cin>>x[b]>>y[b]>>z[b]>>l[b];
     g1=max(x[b],g1);
        g2=max(y[b],g2);
        g3=max(z[b],g3);
     h1=min(x[b]+l[b],h1);
     h2=min(y[b]+l[b],h2);
     h3=min(z[b]+l[b],h3);}
     c=(h1-g1);
     d=(h2-g2);
     e=(h3-g3);
   long long int  f=c*d*e;
if(f>=0) cout<<f;
else cout<<0;
    return 0;
}