HDU1540 (线段树区间合并)

线段树区间合并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;
}

有时间的话把图补一下。。。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值