Hdu 4419 Colourful Rectangle(矩形颜色交)

传送门:Hdu 4419 Colourful Rectangle

---

题意:给你10000个三种颜色的矩形，不同颜色的矩形相互覆盖会形成不同的颜色，问形成的七种颜色的面积是多少。

---

思路:很显然,一种颜色的时候是很容易解决的(面积并)
两种颜色或以上的时候,我们可以在线段树的每个节点开三个标记(表示对应区间的某一种颜色是否出现)
利用二进制状态维护线段树中的面积交

后面利用容斥减去多余部分就可以了
R记为0,G记为1,B记为2

收获:一次扫描线只能维护一种二进制状态
思路二:用0,1,…7八个数记录节点的不同颜色的覆盖情况，每次上传时根据子节点的信息及父节点的覆盖情况算出父节点的信息即可。

---

#include<bits/stdc++.h>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
const int N=1e4+100;
struct node{
    long long L,R,h;
    int flag,d;
}a[N*2],b[2*N];
char s[2];
int mark[8*N][3],t[2*N],m,Sta[8*N];
long long area[8*N][8],Ans[8];

bool cmp(const node& u,const node& v){
    if(u.h==v.h)
        return u.d>v.d;
    return u.h<v.h;
}

int L,R,d,flag,sta;

void pushup(int rt,int l,int r){
    if(Sta[rt]==sta){
        for(int i=1;i<=sta;i++)
            area[rt][i]=t[r+1]-t[l];
    }
    else if(l==r){
        int j=Sta[rt];
        memset(area[rt],0,sizeof(area[rt]));
        do{ //枚举子集
            area[rt][j]=t[r+1]-t[l];
            j=(j-1)&Sta[rt];
        }while(Sta[rt]!=j);
        area[rt][0]=t[r+1]-t[l];
    }
    else{   //其它状态
        for(int i=1;i<=sta;i++) //想得到哪一种状态
            area[rt][i]=area[rt<<1][i-(Sta[rt]&i)]+area[rt<<1|1][i-(Sta[rt]&i)];
    }
}

void update(int l,int r,int rt){
    if(L<=l&&R>=r){
        if(mark[rt][flag]==0)   //从0->1
            Sta[rt]^=(1<<flag);
        mark[rt][flag]+=d;
        if(mark[rt][flag]==0)   //从1->0
            Sta[rt]^=(1<<flag);
        pushup(rt,l,r);
        return ;
    }
    int m=(l+r)>>1;
    if(L<=m)
        update(lson);
    if(R>m)
        update(rson);
    pushup(rt,l,r);
}

void build(int l,int r,int rt){
    for(int i=1;i<=sta;i++)
        area[rt][i]=0;
    area[rt][0]=t[r+1]-t[l];
    memset(mark[rt],0,sizeof(mark[rt]));
    Sta[rt]=0;
    if(l==r)    return ;
    int m=l+r>>1;
    build(lson);
    build(rson);
}

long long solve(int n){
    build(1,m,1);
    long long ans=0,cnt=0;
    for(int i=1;i<=2*n;i++)
        if((1<<a[i].flag)&sta)
            b[++cnt]=a[i];
    for(int i=1;i<cnt;i++){
        L=lower_bound(t+1,t+m+1,b[i].L)-t,R=lower_bound(t+1,t+m+1,b[i].R)-t-1,d=b[i].d,flag=b[i].flag;
        if(L<=R)
            update(1,m,1);
        ans+=area[1][sta]*(b[i+1].h-b[i].h);
    }
    return ans;
}

int main(){
    int _,n;
    scanf("%d",&_);
    for(int case1=1;case1<=_;case1++){
        scanf("%d",&n);
        int flag,x1,x2,y1,y2;
        for(int i=1;i<=n;i++){
            scanf("%s",s);
            if(s[0]=='R')   flag=0;
            if(s[0]=='G')   flag=1;
            if(s[0]=='B')   flag=2;
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            t[i*2-1]=x1,t[i*2]=x2;
            a[i*2-1]=(node){x1,x2,y1,flag,1};
            a[i*2]=(node){x1,x2,y2,flag,-1};
        }
        sort(t+1,t+2*n+1);
        m=unique(t+1,t+2*n+1)-t-1;
        t[m+1]=t[m];
        sort(a+1,a+2*n+1,cmp);
        printf("Case %d:\n",case1);
        for(int i=1;i<=7;i++)
            sta=i,Ans[i]=solve(n);
        Ans[3]-=Ans[7],Ans[5]-=Ans[7],Ans[6]-=Ans[7];
        Ans[1]-=(Ans[3]+Ans[5]+Ans[7]);
        Ans[2]-=(Ans[3]+Ans[6]+Ans[7]);
        Ans[4]-=(Ans[6]+Ans[5]+Ans[7]);
        printf("%I64d\n%I64d\n%I64d\n%I64d\n%I64d\n%I64d\n%I64d\n",Ans[1],Ans[2],Ans[4],Ans[3],Ans[5],Ans[6],Ans[7]);
    }
    return 0;
}