ZOJ - 3324 Machine

线段树还真是纠结了我好几天。由于几乎0基础，所以先研究别人AC的代码。

链接如下：http://hi.baidu.com/forverlin1204/blog/item/5b0b533110d9d595a8018e5b.html

 

开始不懂：为什么一个区间的高度不是0，就直接不去判断它的子区间中是否有满足题意的连续块了。

 

后来才发现原来是题意理解错了，一个push对应一个release。所以如果父区间的height不为0的话，子区间的height就不可能为0。

 

最后自己跌跌撞撞得写好了代码， 发现提交WA。对比那AC的代码发现，自己的代码就差：对于离散化后的点，如果两点间的间距大于1的

 

话，那么要在这两个点之间加入这两个点的终点。

 

表示无法理解为什么这么做，感觉离散化后的点就能正确地表示每个区间了的。所以为什么要加个点呢？后来自己在纸上模拟了下两种不同表示时候的插入情况。才发现如果不加这个点的话， 那么离散化后，原先不连续的两点就被当成连续的了。

 

比如输入的坐标包含1, 2, 4, 5。

 

转化后为0(1)， 1(2)， 2(4)， 3(5)。于是乎，原先的2,4点变成1,2点后被当做连续的了。所以我们只需这两个点之间随便地插入一个点。

 

下面说一下解题的细节：

 

一、首先树的构造：

 

struct Node{

int height, right, left, cnt, r, l;

};

height: 区间的高度。push则加1，release则减一。(从图像上看好像有点有违常理)

right, left:记录区间的最右边（左边）是不是属于某个连续块的。

cnt: 记录这个区间里有多少个高度为0的连续块。

r, l: 区间的左右界限。

 

 

二、接着是如何标记

 

我们只需标记找到的能覆盖的最大区间就行。

由于一个区间的height为0的话肯定没有符合题意的连续块，所以push或release一个区间后，如果heigh不为0，则将它的cnt以及right和

 

left都赋值为0。

 

三、如何更新（传递标记）

 

做法总是递归更新，问题就在于如何根据子区间的状态来更新父区间。

 

首先，如果父区间的高度不为0的话，那么是不用考虑子区间的具体情况的。cnt总为0。

 

如果父区间高度为0，那么就要考虑子区间了。cnt等于两个子区间的cnt相加。还要注意减去重复计算的（如果有的话）。

 

还有就是要把子区间的right，left标记传上来。 最后，更新时还要留意区间是不是叶子节点的情况。

 

 

ps.这回终于体会到什么是真正的离散化了：只记录需要的。 之前做的那个求线段并的，只是离散化初体验，感觉不到明显的离散—_—

 

 

#include <cstdio> #include <algorithm> using namespace std; #define MAXN 40010 struct Op { int l, r; char c; }; struct Node{ int height, right, left, cnt, r, l; }; Node tree[MAXN*4]; Op p[MAXN]; int t[MAXN], tt[MAXN]; void build(int k, int l, int r) { tree[k].right = tree[k].left = 1; tree[k].l = l; tree[k].r = r; tree[k].height = 0; tree[k].cnt = 1; if(l == r) return; int mid = (l+r)/2;; build(k<<1, l, mid); build(k<<1|1, mid+1, r); } //当区间的高度值为0时，才更新该区间。 void update(int k) { if(tree[k].l == tree[k].r) { tree[k].cnt = 1; tree[k].left = tree[k].right = 1; return; } int lson = k<<1, rson = k<<1|1; tree[k].cnt = tree[lson].cnt + tree[rson].cnt; if(tree[lson].right && tree[rson].left) { tree[k].cnt--; } //根据子区间是否连续，更新该区间的连续性 tree[k].left = tree[lson].left; tree[k].right = tree[rson].right; } void add(int k, int l, int r) { int lson = k<<1, rson = k<<1|1; if(l<=tree[k].l && tree[k].r<=r) { tree[k].height++; //if(tree[k].height == 0) { // update(k); //} //else { tree[k].left = tree[k].right = 0; tree[k].cnt = 0; //} return; } int mid = (tree[k].l+tree[k].r)>>1; if(l <= mid) { add(lson, l, r); } if(r > mid) { add(rson, l, r); } if(tree[k].height == 0) { //当tree[k].height值不为0是，那么第k个区间肯定是不连续的。 update(k); } } void del(int k, int l, int r) { int lson = k<<1, rson = k<<1|1; if(l<=tree[k].l && tree[k].r<=r) { tree[k].height--; if(tree[k].height == 0) { update(k); } else { tree[k].left = tree[k].right = 0; tree[k].cnt = 0; } return; } int mid = (tree[k].l+tree[k].r)>>1; if(l <= mid) { del(lson, l, r); } if(r > mid) { del(rson, l, r); } if(tree[k].height == 0) { update(k); } } int find(int x, int l, int r) { while(l < r) { int mid = (l+r)>>1; if(tt[mid] == x) { return mid; } else if(tt[mid] < x){ l = mid+1; } else { r = mid; } } return -1; } int main() { int T, cas = 1, tlen, ttlen, i; char op[4]; scanf("%d", &T); while(T--) { printf("Case #%d:/n", cas++); int n, m; tlen = 0; scanf("%d%d", &n, &m); for(i = 0; i < m; i++) { scanf("%s%d%d", op, &p[i].l, &p[i].r); p[i].c = op[0]; t[tlen++] = p[i].l; t[tlen++] = p[i].r; } t[tlen++] = 0; t[tlen++] = n-1; sort(t, t+tlen); ttlen = 1; for(i = 1; i < tlen; i++) { if(t[i] != t[i-1]) { t[ttlen++] = t[i]; } } tlen = ttlen; ttlen = 0; tt[ttlen++] = t[0]; for(i = 1; i < tlen; i++) { //如果不加的话，那么原先序列1，2， 3，5转化为0, 1, 2, 3后。 //2和3 原先对应的3 和 5就被当成连续块了 if(t[i] > t[i-1]+1) tt[ttlen++] = (t[i]+t[i-1])>>1; tt[ttlen++] = t[i]; } build(1, 0, ttlen-1); for(i = 0; i < m; i++) { int l = find(p[i].l, 0, ttlen); int r = find(p[i].r, 0, ttlen); if(p[i].c == 'p') { add(1, l, r); } else { del(1, l, r); } printf("%d/n", tree[1].cnt); } } return 0; }