hdu5126 stars

题目描述

题解：

和二维的比起来差不多。

但是这是四维偏序。

所以搞一下CDQ套CDQ。

CDQ是维度a已经有序，按维度b排序，然后将维度c存入一维数据结构。

所以我们在第一层CDQ中分治处理，将合法的前一半打标记。

然后进入第二层CDQ，处理打标记的点对没打标记的点的影响。

可以说是将两维压成一维。

想法仍然是容斥。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 50050
#define ll long long
inline int rd()
{
    int f=,c=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){c=*c+ch-'';ch=getchar();}
    return f*c;
}
struct node
{
    int x,y,z,t,w,k,id;
    int typ;
    node(){}
    node(int x,int y,int t,int w,int k,int i):x(x),y(y),t(t),w(w),k(k),id(i){}
}p[N<<],tmp[N<<],sec[N<<];
bool cmpx(node a,node b)
{
    return a.x<b.x;
}
bool cmpid(node a,node b)
{
    return a.id<b.id;
}
int cnt;
struct PR
{
    int x,id;
    PR(){}
    PR(int x,int i):x(x),id(i){}
}Z[N<<];
bool cmp(PR a,PR b)
{
    return a.x<b.x;
}
int T,Q,typ,tim;
int kkk[N<<],ct;
ll ans[N];
struct BIT
{
    ll v[N<<];
    void up(int x,ll d)
    {
        if(!x)return ;
        while(x<(N<<))
            v[x]+=d,x+=(x&-x);
    }
    ll down(int x)
    {
        if(!x)return ;
        ll ret = 0ll;
        while(x)
            ret+=v[x],x-=(x&-x);
        return ret;
    }
}tr;
void init()
{
    cnt=ct=;
    tim=;
}
void ap(int x,int y,int z,int t,int w,int k)
{
    cnt++;
    p[cnt] = node(x,y,t,w,k,cnt);
    Z[cnt] = PR(z,cnt);
}
void Sort(int l,int r)
{
    int mid = (l+r)>>;
    int i=l,j=mid+,k=l;
    while(i<=mid&&j<=r)
    {
        while(i<=mid&&tmp[i].y<=tmp[j].y)
        {
            sec[k]=tmp[i];
            i++,k++;
        }
        while(j<=r&&tmp[i].y>tmp[j].y)
        {
            sec[k]=tmp[j];
            j++,k++;
        }
    }
    while(i<=mid)
    {
        sec[k]=tmp[i];
        i++,k++;
    }
    while(j<=r)
    {
        sec[k]=tmp[j];
        j++,k++;
    }
    for(i=l;i<=r;i++)
        tmp[i]=sec[i];
}
void cdq2(int l,int r)
{
    if(l==r)return ;
    int mid = (l+r)>>;
    cdq2(l,mid);cdq2(mid+,r);
    Sort(l,mid);Sort(mid+,r);
    int i,j;
    for(i=l,j=mid+;j<=r;j++)
    {
        while(i<=mid&&tmp[i].y<=tmp[j].y)
        {
            if(tmp[i].typ&&tmp[i].w)tr.up(tmp[i].z,tmp[i].w);
            i++;
        }
        if(kkk[tmp[j].t]&&!tmp[j].typ)ans[kkk[tmp[j].t]]+=1ll*tmp[j].k*tr.down(tmp[j].z);
    }
    for(i=i-;i>=l;i--)
        if(tmp[i].typ&&tmp[i].w)tr.up(tmp[i].z,-tmp[i].w);
}
void cdq1(int l,int r)
{
    if(l==r)return ;
    int mid = (l+r)>>;
    cdq1(l,mid),cdq1(mid+,r);
    sort(p+l,p+mid+,cmpx);
    sort(p+mid+,p+r+,cmpx);
    int i=l,j=mid+,k=l;
    while(i<=mid&&j<=r)
    {
        while(i<=mid&&p[i].x<=p[j].x)
        {
            tmp[k]=p[i],tmp[k].typ=;
            i++,k++;
        }
        while(j<=r&&p[i].x>p[j].x)
        {
            tmp[k]=p[j],tmp[k].typ=;
            j++,k++;
        }
    }
    while(i<=mid)
    {
        tmp[k]=p[i],tmp[k].typ=;
        i++,k++;
    }
    while(j<=r)
    {
        tmp[k]=p[j],tmp[k].typ=;
        j++,k++;
    }
    cdq2(l,r);
}
int main()
{
    T=rd();
    while(T--)
    {
        init();
        Q = rd();
        for(int ax,ay,az,bx,by,bz,i=;i<=Q;i++)
        {
            typ=rd();
            if(typ==)
            {
                ax=rd(),ay=rd(),az=rd();
                ap(ax,ay,az,tim,,);
            }else
            {
                ax=rd(),ay=rd(),az=rd(),bx=rd(),by=rd(),bz=rd();
                tim++,ct++;
                kkk[tim]=ct;
                ax--,ay--,az--;
                ap(ax,ay,az,tim,,-);
                ap(ax,ay,bz,tim,,);
                ap(ax,by,az,tim,,);
                ap(ax,by,bz,tim,,-);
                ap(bx,ay,az,tim,,);
                ap(bx,ay,bz,tim,,-);
                ap(bx,by,az,tim,,-);
                ap(bx,by,bz,tim,,);
                tim++;
            }
        }
        sort(Z+,Z++cnt,cmp);
        for(int las=-,k=,i=;i<=cnt;i++)
        {
            if(las!=Z[i].x)
            {
                las = Z[i].x;
                k++;
            }
            p[Z[i].id].z = k;
        }
        cdq1(,cnt);
        for(int i=;i<=ct;i++)
            printf("%I64d\n",ans[i]);
        for(int i=;i<=ct;i++)
            ans[i]=;
    }
    return ;
}