码迷,mamicode.com
首页 > 其他好文 > 详细

HDU 6183 Color it cdq分治 + 线段树 + 状态压缩

时间:2017-09-01 11:52:20      阅读:291      评论:0      收藏:0      [点我收藏+]

标签:log   line   orm   define   fine   oid   point   颜色   占用   

Color it

Time Limit: 20000/10000 MS (Java/Others)    Memory Limit: 132768/132768 K (Java/Others)


Problem Description
Do you like painting? Little D doesn‘t like painting, especially messy color paintings. Now Little B is painting. To prevent him from drawing messy painting, Little D asks you to write a program to maintain following operations. The specific format of these operations is as follows.

0 : clear all the points.

1 x y c : add a point which color is c at point (x,y).

2 x y1 y2 : count how many different colors in the square (1,y1) and (x,y2). That is to say, if there is a point (a,b) colored c, that 1ax and y1by2, then the color c should be counted.

3 : exit.
 
Input
The input contains many lines. 

Each line contains a operation. It may be ‘0‘, ‘1 x y c‘ ( 1x,y106,0c50 ), ‘2 x y1 y2‘ (1x,y1,y2106 ) or ‘3‘. 

x,y,c,y1,y2 are all integers.

Assume the last operation is 3 and it appears only once.

There are at most 150000 continuous operations of operation 1 and operation 2. 

There are at most 10 operation 0. 

 
Output
For each operation 2, output an integer means the answer .
 
Sample Input
0
1 1000000 1000000 50 
1 1000000 999999 0
1 1000000 999999 0
1 1000000 1000000 49
2 1000000 1000000 1000000
2 1000000 1 1000000 0
1 1 1 1
2 1 1 2
1 1 2 2
2 1 1 2
1 2 2 2
2 1 1 2
1 2 1 3
2 2 1 2
2 10 1 2
2 10 2 2
0 1 1 1
1 2 1 1
1 1 1 2
1 2 1 1
2 1 2 2
1 2 1 1
2 1 2 1
1 2 2 1
2 2 10
1 2 2 10
2 2
3
 
Sample Output
2 3 1 2 2 3 3 1 1 1 1 1 1 1
 
题解:
 
  对于加入的点,我把第一维, x 进行排序, cdq分治优化时间这一维,其余部分用线段树
  因为只有50中颜色,我将每一位颜色进行二进制压缩,当作一个数存在线段树里
  利用线段树查询一个区间有多少不同的颜色(位运算)
  
#include <bits/stdc++.h>

inline int read(){int x=0,f=1;char ch=getchar();while(ch<0||ch>9){if(ch==-)f=-1;ch=getchar();}while(ch>=0&&ch<=9){x=x*10+ch-0;ch=getchar();}return x*f;}
using namespace std;
#define ls i<<1
#define rs ls | 1
#define mid ((ll+rr)>>1)
const int N = 1e6 + 10;

namespace IO {
    const int MX = 4e7; //1e7占用内存11000kb
    char buf[MX]; int c, sz;
    void begin() {
        c = 0;
        sz = fread(buf, 1, MX, stdin);
    }
    inline bool read(int &t) {
        while(c < sz && buf[c] != - && (buf[c] < 0 || buf[c] > 9)) c++;
        if(c >= sz) return false;
        bool flag = 0; if(buf[c] == -) flag = 1, c++;
        for(t = 0; c < sz && 0 <= buf[c] && buf[c] <= 9; c++) t = t * 10 + buf[c] - 0;
        if(flag) t = -t;
        return true;
    }
}

struct ss{
    int op,x,y,z,id;
    long long ans;
    ss(int op = 0,int x = 0,int y = 0,int z = 0,int id = 0,long long ans = 0) : op(op), x(x), y(y), z(z), id(id), ans(ans) {}
}Q[N],t[N];

bool cmp(ss s1,ss s2) { if(s1.x == s2.x) return s1.op < s2.op;else return s1.x < s2.x; }

int n,mx,san[N];
void init() {n = 0;}
long long v[N * 20];

void update(int i,int ll,int rr,int x,long long c,int ff) {
    if(ll == rr && x == ll) {if(!ff)v[i] |= 1LL<<c;else v[i] = c;return ;}
    if(x <= mid) update(ls,ll,mid,x,c,ff);
    else update(rs,mid+1,rr,x,c,ff);
    v[i] = v[ls] | v[rs];
}
long long ask(int i,int ll,int rr,int x,int y) {
    if(x > y) return 0;
    if(ll == x && rr == y) return v[i];
    if(y <= mid) return ask(ls,ll,mid,x,y);
    else if(x > mid) return ask(rs,mid+1,rr,x,y);
    else return (ask(ls,ll,mid,x,mid) | ask(rs,mid+1,rr,mid+1,y));
}
void cdq(int ll,int rr) {
    if(ll == rr) return ;
    for(int i = ll; i <= rr; ++i) {
        if(Q[i].id <= mid && Q[i].op == 1)
            update(1,1,mx,Q[i].y,Q[i].z,0);
        else if(Q[i].id > mid && Q[i].op == 2)
            Q[i].ans |= ask(1,1,mx,Q[i].y,Q[i].z);
    }
    for(int i = ll; i <= rr; ++i) {
        if(Q[i].id <= mid && Q[i].op == 1)
            update(1,1,mx,Q[i].y,0,1);
    }
    int L1 = ll, R1 = mid+1;
    for(int i = ll; i <= rr; ++i) {
        if(Q[i].id <= mid) t[L1++] = Q[i];
        else t[R1++] = Q[i];
    }
    for(int i = ll; i <= rr; ++i) Q[i] = t[i];
    cdq(ll,mid);cdq(mid+1,rr);
}
void solve() {
    if(n == 0) return ;
    int cny = 0;
    for(int i = 1; i <= n; ++i) {
        if(Q[i].op == 1) san[++cny] = Q[i].y;
        else {
            san[++cny] = Q[i].y;
            san[++cny] = Q[i].z;
        }
    }
    sort(san+1,san+cny+1);
    int SA = unique(san+1,san+cny+1) - san - 1;
    for(int i = 1; i <= n; ++i) {
        if(Q[i].op == 1) Q[i].y = lower_bound(san+1,san+SA+1,Q[i].y) - san;
        else {
            Q[i].y = lower_bound(san+1,san+SA+1,Q[i].y) - san;
            Q[i].z = lower_bound(san+1,san+SA+1,Q[i].z) - san;
        }
    }
    mx = SA;
    sort(Q+1,Q+n+1,cmp);
    cdq(1,n);
    for(int i = 1; i <= n; ++i) {
        if(Q[i].op == 2) {
            int sum = 0;
            for(int j = 0; j <= 50; ++j)
                if((Q[i].ans >> j) & 1) sum++;
            printf("%d\n",sum);
        }
    }
}
int op,x,z,y;
int main() {
    IO::begin();
    while(1) {
        IO::read(op);
        if(op == 3 || op == 0) {
            solve();
            init();
            if(op == 3) return 0;
            continue;
        }
        IO::read(x);
        IO::read(y);
        IO::read(z);
        Q[++n] = ss(op,x,y,z,n,0);
    }
    return 0;
}

 

 

HDU 6183 Color it cdq分治 + 线段树 + 状态压缩

标签:log   line   orm   define   fine   oid   point   颜色   占用   

原文地址:http://www.cnblogs.com/zxhl/p/7462417.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!