看到讲课安排上 线段树有一节课"区间合并" 我是迷茫的 因为并没有见过
然后了解了一下题目 发现以前写过 还是很麻烦的树链剖分 大概是 解决带修改的区间查询"连续问题"
意思就是给一个数组 要对这个数组进行修改 然后进行区间查询 查询的一般是 l r 区间内的 连续xx 可能是LCIS 也可能只是连续的数字..
解题的方法比较固定 由于是 连续性 相关 所以对于每个数组 维护它的左边与右边与全局的性质 例如紧靠左(右)边连续xx的长度和这个区间内连续xx的最长长度
关键就在up函数上 不管是单点还是区间修改 都能通过up来做巧妙地构造
POJ 3667
每次清空一段 或 查询一段中最靠左的x长度并填补
维护每个区间的左右连续 最简单的区间合并
HDU 3308
单点修改 或 查询一段中的LCIS
与上一道题的差别就是up函数 和 query函数中 对于 一个区间的两个子区间合并 需要判断一下数字的大小关系
HDU 3397
01串 区间赋值 区间反转 查询 区间sum 和 区间连续1的maxlen
同时维护左右全局的0和1的长度 在做反转的时候直接交换就可以了
在down的时候 如果down下去的是一个反转的历史 它可能会和上一次的历史发生化合 上一次的历史是赋值0 那这一次的历史就覆盖成赋值1而不是反转
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<map>
#include<vector>
#include<queue>
using namespace std;
const int maxn = 100050 ;
int n , m ;
struct node {
int rl1, ll1, tl1 ;
int rl0, ll0, tl0 ;
int sum ;
int l , r ;
int len;
}a[maxn * 4];
int mark[maxn * 4] ; /// 0 -> 0 1 -> 0 -1 notdid 2 -> ycy
int b[maxn] ;
void up(int p) {
if(a[p].l == a[p].r )return ;
int pl = p*2 , pr = p*2+1 ;
a[p].tl0 = max(a[pl].tl0 , a[pr].tl0) ;
a[p].tl0 = max(a[p].tl0 , a[pl].rl0 + a[pr].ll0) ;
a[p].tl1 = max(a[pl].tl1 , a[pr].tl1) ;
a[p].tl1 = max(a[p].tl1 , a[pl].rl1 + a[pr].ll1) ;
a[p].ll0 = a[pl].ll0;
if(a[pl].sum == 0 ) {
a[p].ll0 = max(a[p].ll0 , a[pl].len + a[pr].ll0) ;
}
a[p].rl0 = a[pr].rl0;
if(a[pr].sum == 0 ) {
a[p].rl0 = max(a[p].rl0 , a[pr].len + a[pl].rl0) ;
}
a[p].ll1 = a[pl].ll1;
if(a[pl].sum == a[pl].len ) {
a[p].ll1 = max(a[p].ll1 , a[pl].len + a[pr].ll1) ;
}
a[p].rl1 = a[pr].rl1;
if(a[pr].sum == a[pr].len ) {
a[p].rl1 = max(a[p].rl1 , a[pr].len + a[pl].rl1) ;
}
a[p].sum = a[pl].sum + a[pr].sum ;
return ;
}
void build(int p,int l,int r){
mark[p] = -1 ;
a[p].l = l , a[p].r = r;
if(b[l] == 0) a[p].ll0 = 1 , a[p].ll1 = 0; else a[p].ll0 = 0 , a[p].ll1 = 1;
if(b[r] == 0) a[p].rl0 = 1 , a[p].rl1 = 0; else a[p].rl0 = 0 , a[p].rl1 = 1;
if(b[l] == 0) a[p].tl0 = 1 , a[p].tl1 = 0; else a[p].tl0 = 0 , a[p].tl1 = 1;
a[p].len = r - l + 1 ;
if(b[l] == 0)a[p].sum = 0 ; else a[p].sum = 1 ;
if(l == r)return ;
int mid = (l+r) / 2 ;
build(p*2,l,mid);
build(p*2+1,mid+1,r);
up(p);
}
void ycy(int p) {
swap(a[p].ll0 , a[p].ll1) ;
swap(a[p].rl0 , a[p].rl1) ;
swap(a[p].tl0 , a[p].tl1) ;
a[p].sum = a[p].len - a[p].sum ;
}
void down(int p) {
if(a[p].l == a[p].r)return ;
if(mark[p] == -1)return ;
int pl = p*2 , pr = p*2+1 ;
if(mark[p] == 0) {
mark[pl] = mark[pr] = mark[p] ;
a[pl].ll0 = a[pl].rl0 = a[pl].tl0 = a[pl].len ;
a[pl].ll1 = a[pl].rl1 = a[pl].tl1 = a[pl].sum = 0 ;
a[pr].ll0 = a[pr].rl0 = a[pr].tl0 = a[pr].len ;
a[pr].ll1 = a[pr].rl1 = a[pr].tl1 = a[pr].sum = 0 ;
}
else if(mark[p] == 1){
mark[pl] = mark[pr] = mark[p] ;
a[pl].ll0 = a[pl].rl0 = a[pl].tl0 = 0 ;
a[pl].ll1 = a[pl].rl1 = a[pl].tl1 = a[pl].sum = a[pl].len ;
a[pr].ll0 = a[pr].rl0 = a[pr].tl0 = 0 ;
a[pr].ll1 = a[pr].rl1 = a[pr].tl1 = a[pr].sum = a[pr].len ;
}
else {
if(mark[pl] == 0) mark[pl] = 1;
else if(mark[pl] == 1) mark[pl] = 0 ;
else if(mark[pl] == 2) mark[pl] = -1 ;
else mark[pl] = 2;
if(mark[pr] == 0) mark[pr] = 1;
else if(mark[pr] == 1) mark[pr] = 0 ;
else if(mark[pr] == 2) mark[pr] = -1 ;
else mark[pr] = 2;
ycy(pl) , ycy(pr) ;
}
mark[p] = -1 ;
}
void upda1(int p , int l , int r ){ /// 0
if(a[p].l >= l && a[p].r <= r) {
mark[p] = 0 ;
a[p].ll0 = a[p].rl0 = a[p].tl0 = a[p].len ;
a[p].ll1 = a[p].rl1 = a[p].tl1 = a[p].sum = 0 ;
return ;
}
if(a[p].l == a[p].r)return ;
down(p) ;
int mid = (a[p].l + a[p].r) / 2 ;
if(r <= mid) {
upda1(p*2,l,r) ;
}
else if(l >= mid+1) {
upda1(p*2+1,l,r);
}
else {
upda1(p*2,l,mid);
upda1(p*2+1,mid+1,r);
}
up(p);
}
void upda2(int p , int l , int r ){ /// 1
if(a[p].l >= l && a[p].r <= r) {
mark[p] = 1 ;
a[p].ll0 = a[p].rl0 = a[p].tl0 = 0 ;
a[p].ll1 = a[p].rl1 = a[p].tl1 = a[p].sum = a[p].len ;
return ;
}
if(a[p].l == a[p].r)return ;
down(p) ;
int mid = (a[p].l + a[p].r) / 2 ;
if(r <= mid) {
upda2(p*2,l,r) ;
}
else if(l >= mid+1) {
upda2(p*2+1,l,r);
}
else {
upda2(p*2,l,mid);
upda2(p*2+1,mid+1,r);
}
up(p);
}
void upda3(int p, int l , int r ){
down(p);
if(a[p].l >= l && a[p].r <= r) {
ycy(p) ;
mark[p] = 2 ;
return ;
}
if(a[p].l == a[p].r) return ;
int mid = (a[p].l + a[p].r) / 2;
if(r <= mid) {
upda3(p*2,l,r);
}
else if(l >= mid +1) {
upda3(p*2+1,l,r);
}
else {
upda3(p*2,l,mid);
upda3(p*2+1,mid+1,r);
}
up(p);
}
int query1(int p,int l,int r){
down(p);
if(a[p].l>=l&&a[p].r<=r) return a[p].sum;
int mid = (a[p].l+a[p].r)/2;
if(r<=mid)return query1(p*2,l,r);
else if(l>=mid+1)return query1(p*2+1,l,r);
else return query1(p*2,l,r) + query1(p*2+1,l,r);
}
int query2(int p,int l,int r){
down(p);
if(a[p].l>=l&&a[p].r<=r)return a[p].tl1 ;
if(a[p].l == a[p].r)return a[p].tl1;
int mid = (a[p].l + a[p].r) / 2 ;
if(r <= mid) {
return query2(p*2,l,r);
}
else if(l >= mid+1)
return query2(p*2+1,l,r);
else {
int res = 0;
int r1 = r ;
int r2 = a[p*2+1].l + a[p*2+1].ll1 - 1;
r1 = min(r2,r1) ;
int l1 = l ;
int l2 = a[p*2].r - a[p*2].rl1 + 1;
l1 = max(l1,l2) ;
if(r1 >= l1)
res = max(res , r1-l1+1);
res = max(query2(p*2,l,r),res) ;
res = max(query2(p*2+1,l,r),res);
return res ;
}
}
int main(){
int t;
scanf("%d",&t);
while(t-- ){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)scanf("%d",&b[i]) ;
build(1,1,n) ;
for(int i=1;i<=m;i++){
int op;
scanf("%d",&op);
int A,B;
scanf("%d%d",&A,&B);
A++,B++;
if(op==0){
upda1(1,A,B);
}
else if(op==1){
upda2(1,A,B);
}
else if(op==2){
upda3(1,A,B);
}
else if(op==3){
printf("%d
",query1(1,A,B));
}
else {
printf("%d
",query2(1,A,B));
}
}
}
}
感觉区间合并是一个线段树的思路
学习了kd-tree之后对线段树有了新的认识