hdu 3642 Get The Treasure

　　对于这题，一开始看到z值较小，就像枚举z坐标，然后对于一个单独的z平面做面积交来累加答案。。(三个以上柱体的体积交)

　　搜题解学完up函数的姿势后感觉自己写的还是蛮不错的嘛。。

　　线段树每个节点包含当前节店的覆盖次数cnt，表示覆盖次数-长度映射的数组，由于只要cnt>=3就可以统计，故数组只开到3。。

　　由于n只有1000，所以平面交的时候先对于x坐标进行离散化，然后就是扫描线。。

　　一开始自己的姿势是对于每个z建立一个线段树，但是会mle。。。后来才改成使用一个重复利用。。然后就是没注意到cnt会大于三可能会re到无法自拔（比如我）。。

　　

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define mod 1000000007
#define LL 2*n
#define RR 2*n+1
#define lson l,mid,2*n
#define rson mid+1,r,2*n+1
#define mlt map<ll,ll>::iterator
#define mit map<int,int>::iterator
#define vi vector<int>
#define mii map<int,int>
#define msi map<string,int>
#define si set<int>
#define sit set<int>::iterator
#define slt set<l>::iterator
#define mp make_pair
#define pb push_back
#define pii pair<int,int>
#define pdd pair<double,double>
#define pi acos(-1)
const int maxn =2222;
struct Seg{
    int l , r , h , s;
    Seg() {}  
    Seg(ll a,ll b,ll c,ll d):l(a) , r(b) , h(c) , s(d) {}
    bool operator < (const Seg &cmp) const {
    if (h == cmp.h) return s > cmp.s;
    return h < cmp.h;
    }
}ss[maxn<<2];
int x[maxn<<2];
struct node
{
    int cnt,a[5];
}tr[maxn<<2];
int bin(int z,int n)
{
    int l=0,r=n;
    while(r-l>1)
    {
        int mid=(l+r)>>1;
        if(x[mid]>=z)
        r=mid;
        else
        l=mid;
    }
    return r;
}
void up(int l,int r,int n)
{
    memset(tr[n].a,0,sizeof(tr[n].a));
    int cnt=min(3,tr[n].cnt);
    if(l==r)
    {
        tr[n].a[cnt]=x[r+1]-x[l];
        return;
    }
    int tmp=x[r+1]-x[l];
    for(int i=0;i<=3;i++)
    {
        tr[n].a[min(cnt+i,3)]+=tr[LL].a[i]+tr[RR].a[i];
    }
    for(int i=0;i<=3;i++)
    tmp-=tr[n].a[i];
    tr[n].a[cnt]+=tmp;
}
void update(int l,int r,int n,int left,int right,int targ)
{
    if(left<=l&&r<=right)
    {
        tr[n].cnt+=targ;
        up(l,r,n);
        return ;
    }
    int mid=(l+r)>>1;
    if(left<=mid)
    update(lson,left,right,targ);
    if(right>mid)
    update(rson,left,right,targ);
    up(l,r,n);
}
int a[5000][20];
int main()
{
    int t,cas=1;
   //freopen("C:\Users\Administrator\Documents\duipai\data.txt","r",stdin);freopen("C:\Users\Administrator\Documents\duipai\out1.txt","w",stdout);
    scanf("%d",&t);
    while(t--)
    {
        int n;
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=6;j++)
            scanf("%d",&a[i][j]);
        }
        ll ans=0;
        for(int i=-500;i<=500;i++)
        {
            int num1=0,num=0;
            for(int k=1;k<=n;k++)
            {
                    if(i>=a[k][3]&&i<a[k][6])
                    {
                        ss[++num]=Seg(a[k][1],a[k][4],a[k][2],1);
                        ss[++num]=Seg(a[k][1],a[k][4],a[k][5],-1);
                    }
            }
            sort(ss+1,ss+1+num);
            for(int j=1;j<=num;j++)
            {
                x[++num1]=ss[j].l;
                x[++num1]=ss[j].r;
            }
            sort(x+1,x+1+num1);
            int tmp=1;
            for(int j=2;j<=num1;j++)
            {
                if(x[j]!=x[j-1])
                x[++tmp]=x[j];
            }
            num1=tmp;
            for(int j=1;j<=num;j++)
            {
                int l=bin(ss[j].l,num1);
                int r=bin(ss[j].r,num1)-1;
                if(l<=r)
                update(1,num1-1,1,l,r,ss[j].s);
                ans+=1ll*tr[1].a[3]*(ss[j+1].h-ss[j].h);
            }
        }
        printf("Case %d: %I64d
",cas++,ans);
    }
    return 0;
 }