线段树区间合并dalao总结:ino
(以下以求01序列中最长连续1为例)
tree[cur].left代表以区间左端点为起点的连续段的长度 tree[cur].right代表右边 tree[cur].all代表该区间内最长的连续段
通过这三个变量的组合 对于每一个区间一般有四个值供我们使用
1 tree[cur].left代表以该区间左端点为起点的连续段的长度(左连续段)
2 tree[cur].right代表以该区间右端点为终点的连续段的长度(右连续段)
3 tree[2cur].right+tree[2cur+1].left 代表包含该区间中点的连续段的长度 这是个隐含值(中间连续段)
4 tree[cur].all代表该区间内最长连续段
其值不止要从左右子节点的最长连续段中择最大 还要考虑该区间的中间连续段即
tree[cur].all=max(tree[2cur].right+tree[2cur+1].left,max(tree[2cur].all,tree[2cur+1].all));
对于tree[cur].all的取值可以这样理解 如果两个子区间合并后产生了新的连续段
那一定是tree[2cur].right+tree[2cur+1].left的值
即中间连续段 如果没有那就只能在两个子区间的最长连续段择优了
dalao讲解:%%%
这里说一下蒟蒻在看完dalao题解的理解:
在更新(insert)的时候首先递归到单一结点,然后再向上更新父节点,
父节点的区间最大连续区间就是在左,右儿子节点上取最大值,以及左儿子右端点和右儿子左端点的和之间取最大值。。。。
以及
左子树区间满了的话,父亲左区间要加上右孩子的左区间
if (a[2 * i].ls == a[2 * i].r - a[2 * i].l + 1)
a[i].ls += a[2 * i + 1].ls;
右儿子同理。。
查询的时候有可能
//因为t<=mid,看左子树,a[2×i].r-a[2×i].rs+1代表左子树右边连续区间的左边界值,
//如果t在左子树的右区间内,则要看右子树的左区间有多长并返回
代码是搬运的。。。
#include<iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#include <stdlib.h>
using namespace std;
const int maxn = 50000 + 10;
int n, m;
int s[maxn], top;//s为模拟栈
struct node
{
int l, r;
int ls, rs, ms;
//ls,左端最大连续区间,rs右端最大连续区间,ms区间内最大连续区间
} a[maxn << 2];
void init(int l, int r, int i)
{
a[i].l = l;
a[i].r = r;
a[i].ls = a[i].rs = a[i].ms = r - l + 1;
if (l != r)
{
int mid = (l + r) >> 1;
init(l, mid, i * 2);
init(mid + 1, r, 2 * i + 1);
}
}
void insert(int i, int t, int x)
{
if (a[i].l == a[i].r)
{
if (x == 1)
a[i].ls = a[i].rs = a[i].ms = 1;//修复
else
a[i].ls = a[i].rs = a[i].ms = 0;//破坏
return;
}
int mid = (a[i].l + a[i].r) >> 1;
if (t <= mid)
insert(2 * i, t, x);
else
insert(2 * i + 1, t, x);
a[i].ls = a[2 * i].ls;//左区间
a[i].rs = a[2 * i + 1].rs;//右区间
a[i].ms = max(max(a[2 * i].ms, a[2 * i + 1].ms), a[2 * i].rs + a[2 * i + 1].ls);
//父亲区间内的最大区间必定是,左子树最大区间,右子树最大区间,
//左右子树合并的中间区间,三者中最大的区间值
if (a[2 * i].ls == a[2 * i].r - a[2 * i].l + 1)
//左子树区间满了的话,父亲左区间要加上右孩子的左区间
a[i].ls += a[2 * i + 1].ls;
if (a[2 * i + 1].rs == a[2 * i + 1].r - a[2 * i + 1].l + 1)//同理
a[i].rs += a[2 * i].rs;
}
int query(int i, int t)
{
if (a[i].l == a[i].r || a[i].ms == 0 || a[i].ms == a[i].r - a[i].l + 1)
//到了叶子节点或者该访问区间为空或者已满都不必要往下走了
return a[i].ms;
int mid = (a[i].l + a[i].r) >> 1;
if (t <= mid)
{
if (t >= a[2 * i].r - a[2 * i].rs + 1)
//因为t<=mid,看左子树,a[2*i].r-a[2*i].rs+1代表左子树右边连续区间的左边界值,
//如果t在左子树的右区间内,则要看右子树的左区间有多长并返回
return query(2 * i, t) + query(2 * i + 1, mid + 1);
else
return query(2 * i, t);
//如果不在左子树的右边界区间内,则只需要看左子树
}
else
{
if (t <= a[2 * i + 1].l + a[2 * i + 1].ls - 1)//同理
return query(2 * i + 1, t) + query(2 * i, mid);
else
return query(2 * i + 1, t);
}
}
int main()
{
int i, j, x;
char ch[2];
while (~scanf("%d%d", &n, &m))
{
top = 0;
init(1, n, 1);
while (m--)
{
scanf("%s", ch);
if (ch[0] == 'D')
{
scanf("%d", &x);
s[top++] = x;
insert(1, x, 0);
}
else if (ch[0] == 'Q')
{
scanf("%d", &x);
printf("%d\n", query(1, x));
}
else
{
if (x>0)
{
x = s[--top];
insert(1, x, 1);
}
}
}
}
return 0;
}
有时间的话把图补一下。。。