Hdu 4419 Colourful Rectangle(矩形颜色交)

关注

原创最新推荐文章于 2026-02-03 21:28:05 发布 576 阅读 AI 写同款· GEO 优化 ›

MathType 7.4保姆级安装教程：从下载到Word/WPS完美加载（附常见问题解决）本文提供MathType 7.4的详细安装教程，涵盖从下载到Word/WPS完美加载的全过程，并附常见问题解决方案。特别针对科研写作中的公式编辑需求，介绍如何避免安装陷阱、配置Office集成及高效排版技巧，帮助用户快速掌握这一公式编辑神器。阅读详情

传送门: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;
}