hdu3642扫描线 长方体

立方体交，自己写的莫名其妙MLE了，不知道为什么

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define maxn 2010
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define ll long long
struct cube{
    int x1,y1,z1,x2,y2,z2;
    cube(){}
    cube(int a,int b,int c,int d,int e,int f):
        x1(a),y1(b),z1(c),x2(d),y2(e),z2(f){}
}cubes[maxn];
struct seg{//水平横线段
    int l,r,h,d;
    seg(){}
    seg(int l,int r,int h,int d):l(l),r(r),h(h),d(d){}
    bool operator<(const seg & a)const
    { return h<a.h; }
}segs[maxn];
int tots;
int y[maxn],z[maxn],tot,toty,totz;//y轴，z轴上的数据
int len1[maxn<<],len2[maxn<<],len3[maxn<<];//线段树区间被覆盖了1，2，3次的长度
int cnt[maxn<<];//区间被完全覆盖的次数
inline void pushup(int rt,int l,int r){
    if(cnt[rt]>=){
        len3[rt]=y[r+]-y[l];
        len2[rt]=len1[rt]=;
    }
    else if(cnt[rt]==){
        len1[rt]=;
        len2[rt]=y[r+]-y[l];
        if(l==r) len3[rt]=;//没有子区间的情况
        else {
            len3[rt]=len3[rt<<]+len3[rt<<|]+len2[rt<<]+len2[rt<<|]+len1[rt<<]+len1[rt<<|];
            len2[rt]-=len3[rt];
        }
    }
    else if(cnt[rt]==){
        len1[rt]=y[r+]-y[l];
        if(l==r) len3[rt]=len2[rt]=;//没有子区间的情况
        else {
            len3[rt]=len3[rt<<]+len3[rt<<|]+len2[rt<<]+len2[rt<<|];
            len2[rt]=len1[rt<<]+len1[rt<<|];
            len1[rt]-=len2[rt]+len3[rt];
        }
    }
    else {//cnt[rt]==0
        if(l==r) len3[rt]=len2[rt]=len1[rt]=;
        else {
            len3[rt]=len3[rt<<]+len3[rt<<|];
            len2[rt]=len2[rt<<]+len2[rt<<|];
            len1[rt]=len1[rt<<]+len1[rt<<|];
        }
    }
}
void update(int L,int R,int c,int l,int r,int rt){
    if(L<=l && R>=r){
        cnt[rt]+=c;
        pushup(rt,l,r);
        return;
    }
    int m=l+r>>;
    if(L<=m) update(L,R,c,lson);
    if(R>m) update(L,R,c,rson);
    pushup(rt,l,r);
}
void init(){
    tot=toty=totz=tots=;
    memset(y,,sizeof y);
    memset(z,,sizeof z);
    memset(len1,,sizeof len1);
    memset(len2,,sizeof len2);
    memset(len3,,sizeof len3);
    memset(cnt,,sizeof cnt);
}
int main(){
    int T,n;
    scanf("%d",&T);
    for(int tt=;tt<=T;tt++){
        init();
        scanf("%d",&n);
        for(int i=;i<=n;i++){
            int a,b,c,d,e,f;
            scanf("%d%d%d%d%d%d",&a,&b,&c,&d,&e,&f);
            cubes[i]=cube(a,b,c,d,e,f);
            y[tot]=b;z[tot++]=c;
            y[tot]=e;z[tot++]=f;
        }
        if(n<) {puts("");continue;}

        sort(y,y+tot);toty=unique(y,y+tot)-y;
        sort(z,z+tot);totz=unique(z,z+tot)-z;
        ll res=;
        for(int i=;i<totz-;i++){
            tots=;
            for(int j=;j<=n;j++)
                if(cubes[j].z1<=z[i] && cubes[j].z2>=z[i+]){
                    segs[tots++]=seg(cubes[j].x1,cubes[j].x2,cubes[j].y1,);
                    segs[tots++]=seg(cubes[j].x1,cubes[j].x2,cubes[j].y2,-);
                }
            sort(segs,segs+tots);//将水平横线段排序
            for(int j=;j<tots;j++){//将这些边更新到线段树中
                if(j!=) res+=(z[i+]-z[i])*(segs[j].h-segs[j-].h)*len3[];
                int L=lower_bound(y,y+toty,segs[j].l)-y;
                int R=lower_bound(y,y+toty,segs[j].r)-y-;
                update(L,R,segs[j].d,,toty,);//为什么要计算在前更新在后？
            }
        }
        printf("%lld\n",res);
    }
    return ;
}

kuangbin的板子是可以过的。。