hdu 3642 长方体体积交....线段树

解题思路：

     一开始看题完全没思路，本来以为要建一个三维的线段树，以为暂时不能够动了，但是之后忍不住了，看了解题报告，原来是枚举z轴上面的值就可以了，这样猥琐过都可以。。。看来还是自己观察数据不够到位。。。

     然后针对每一层z轴，建立一个线段树，对覆盖三次以上的面积求和，然后再求体积。。。

  

#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
#include<string.h>
using std::sort;
using std::unique;
const int N = 1111;
struct node
{
    int x1,y1,z1;
    int x2,y2,z2;
}Vo[N];
struct line
{
    int y,x1,x2;
    int s;
    bool operator<(const line &tmp)const
    {
        return y<tmp.y;
    }
}Li[N<<2];
int topl;
int X[N<<2],Z[N<<2],topx,topz;
int cover[N<<2];
__int64 sum1[N<<2],sum2[N<<2],sum3[N<<2];
int Fin(double k,int len)
{
    int l=0,r=len;
    while(l<r)
    {
        int m=(l+r)>>1;
        if(X[m]==k)return m;
        if(X[m]>k)r=m-1;
        else l=m+1;
    }
    return l;
}
void PushUp(int t,int l,int r)
{
     if(cover[t]>=3)
     {
         sum3[t]=sum2[t]=sum1[t]=(__int64)X[r+1]-X[l];
     }
     else
       if(cover[t]==2)
       {

              sum2[t]=sum1[t]=(__int64)X[r+1]-X[l];
              if(l==r)sum3[t]=0;
              else
              sum3[t]=sum1[t<<1]+sum1[t<<1|1];
       }
       else
         if(cover[t]==1)
         {
                 sum1[t]=(__int64)X[r+1]-X[l];
                 if(l==r)sum2[t]=sum3[t]=0;
                 else
                 {
                   sum2[t]=sum1[t<<1]+sum1[t<<1|1];
                   sum3[t]=sum2[t<<1]+sum2[t<<1|1];
                 }
         }
         else
         {
                if(l==r)sum1[t]=sum2[t]=sum3[t]=0;
                else
                {
                  sum1[t]=sum1[t<<1]+sum1[t<<1|1];
                  sum2[t]=sum2[t<<1]+sum2[t<<1|1];
                  sum3[t]=sum3[t<<1]+sum3[t<<1|1];
                }
         }
}
void update(int t,int l,int r,int L,int R,int val)
{
     if(L<=l&&r<=R)
     {
        cover[t]+=val;
        PushUp(t,l,r);
        return ;
     }
     int m=(r+l)>>1;
     if(L<=m)update(t<<1,l,m,L,R,val);
     if(R>m)update(t<<1|1,m+1,r,L,R,val);
     PushUp(t,l,r);
}
int main()
{
    int t,n,T=0;
    scanf("%d",&t);
    while(t--)
    {
          scanf("%d",&n);
          topz=0;
          for(int i=0;i<n;i++)
          {
             scanf("%d%d%d%d%d%d",&Vo[i].x1,&Vo[i].y1,&Vo[i].z1,&Vo[i].x2,&Vo[i].y2,&Vo[i].z2);
             Z[topz++]=Vo[i].z1;
             Z[topz++]=Vo[i].z2;
          }
          sort(Z,Z+topz);
          topz=unique(Z,Z+topz)-Z;
          __int64 ans=0,ret;
          for(int i=0;i<topz-1;i++)
          {
             topl=topx=0;
             for(int j=0;j<n;j++)
             {
                 if(Vo[j].z1<=Z[i]&&Vo[j].z2>Z[i])
                 {
                    X[topx++]=Vo[j].x1;
                    X[topx++]=Vo[j].x2;
                    Li[topl].y=Vo[j].y1;Li[topl].x1=Vo[j].x1;Li[topl].x2=Vo[j].x2;Li[topl].s=1;topl++;
                    Li[topl].y=Vo[j].y2;Li[topl].x1=Vo[j].x1;Li[topl].x2=Vo[j].x2;Li[topl].s=-1;topl++;
                 }
             }
             sort(X,X+topx);
             sort(Li,Li+topl);
             topx=unique(X,X+topx)-X;
             memset(sum1,0,sizeof(sum1));
             memset(sum2,0,sizeof(sum2));
             memset(sum3,0,sizeof(sum3));
             memset(cover,0,sizeof(cover));
             ret=0;
             for(int k=0;k<topl-1;k++)
             {
                 int l=Fin(Li[k].x1,topx-1);
                 int r=Fin(Li[k].x2,topx-1)-1;
                 if(l<=r)update(1,0,topx-1,l,r,Li[k].s);
                 ret+=sum3[1]*(__int64)(Li[k+1].y-Li[k].y);
             }
             ans+=ret*(__int64)(Z[i+1]-Z[i]);
          }
          printf("Case %d: %I64d\n",++T,ans);
    }
    return 0;
}