HDU 4419 Colourful Rectangle（线段树+扫描线)

题目链接

主要是pushup的代码，其他和区间更新+扫描线差不多。

那个区间如果要再刷一层x，那么sum[x][rt] = que[r+1] - que[l];但是如果原本有颜色为i，颜色将会变成i|x，sum[x][rt] 要减去以前的i颜色的部分。sum[i|x][rt]要加上那部分。

这个题还可以用容斥，容斥的话，多次求面积并就可以了，代码直接是模版，对扫描线，还是不熟啊。

 #include <cstdio>
 #include <cstring>
 #include <string>
 #include <algorithm>
 using namespace std;
 #define LL __int64
 #define maxn 20100
 #define lson l , m, rt<<1
 #define rson m+1, r,rt<<1|1
 int que[*maxn];
 int sum[][*maxn];
 int cnt[*maxn][];
 LL ans[];
 struct node
 {
     int lx,rx,y;
     int s;
     node() {}
     node (int a,int b,int c,int d):lx(a),rx(b),y(c),s(d) {}
     bool operator < (const node &S)const
     {
         return y < S.y;
     }
 } mat[*maxn];
 int bin(int x,int n)
 {
     int str,end,mid;
     str = ;
     end = n;
     while(str <= end)
     {
         mid = (str + end)/;
         if(que[mid] == x)
             return mid;
         else if(que[mid] > x)
             end = mid - ;
         else
             str = mid + ;
     }
     return mid;
 }
 void pushup(int rt,int l,int r)
 {
     int i,x,t;
     x = ;
     for(i = ;i <= ;i ++)
     {
         if(cnt[rt][i] > )
         x |= (<<(i-));
     }
     if(x)
     {
         for(i = ;i < ;i ++)
         sum[i][rt] = ;
         sum[x][rt] = que[r+] - que[l];
         for(i = ; i < ; i ++)//注意这里
         {
             if(x != (i|x))
             {
                 t = sum[i][rt<<] + sum[i][rt<<|];
                 sum[x|i][rt] += t;
                 sum[x][rt] -= t;//减去原本的其他颜色
             }
         }
     }
     else if(l == r)
     {
         for(i = ; i < ; i ++)
             sum[i][rt] = ;
     }
     else
     {
         for(i = ; i < ; i ++)
             sum[i][rt] = sum[i][rt<<] + sum[i][rt<<|];
     }
 }
 void update(int L,int R,int c,int l,int r,int rt)
 {
     int m;
     if(l >= L&&r <= R)
     {
         c >  ? cnt[rt][c] ++:cnt[rt][-c] --;
         pushup(rt,l,r);
         return ;
     }
     m = (l+r)>>;
     if(L <= m)
         update(L,R,c,lson);
     if(R > m)
         update(L,R,c,rson);
     pushup(rt,l,r);
 }
 int main()
 {
     int n,cas = ,i,j,l,r,t,num;
     char ch[];
     int a,b,c,d,col;
     scanf("%d",&t);
     while(t --)
     {
         scanf("%d",&n);
         num = ;
         for(i = ; i <= n; i ++)
         {
             scanf("%s%d%d%d%d",ch,&a,&b,&c,&d);
             if(ch[] == 'R')
                 col = ;
             else if(ch[] == 'G')
                 col = ;
             else
                 col = ;
             mat[num] = node(a,c,b,col);
             que[num++] = a;
             mat[num] = node(a,c,d,-col);
             que[num++] = c;
         }
         sort(que,que+num);
         sort(mat,mat+num);
         int k = ;
         for(i = ; i < num; i ++)
         {
             if(que[i] != que[i-])
                 que[k++] = que[i];
         }
         memset(cnt,,sizeof(cnt));
         memset(sum,,sizeof(sum));
         memset(ans,,sizeof(ans));
         for(i = ; i < num-; i ++)
         {
             l = bin(mat[i].lx,k-);
             r = bin(mat[i].rx,k-) - ;
             if(l <= r)
                 update(l,r,mat[i].s,,k-,);
             for(j = ; j < ; j ++)
                 ans[j] += ((LL)mat[i+].y - (LL)mat[i].y)*(LL)sum[j][];
         }
         printf("Case %d:\n",cas++);
         printf("%I64d\n",ans[]);
         printf("%I64d\n",ans[]);
         printf("%I64d\n",ans[]);
         printf("%I64d\n",ans[]);
         printf("%I64d\n",ans[]);
         printf("%I64d\n",ans[]);
         printf("%I64d\n",ans[]);
     }
     return ;
 }