（转自一位神牛，收藏起来以后留着看）hdu 1540 Tunnel Warfare（线段树区间合并）

题意：D代表破坏村庄，R代表修复最后被破坏的那个村庄，Q代表询问包括x在内的最大连续区间是多少

思路：代码中的注释已经比较详细了，所以不多解释什么，在线段树的区间内，我们要用三个变量记录左边连续区间，右边连续区间和最大连续区间

代码：

#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;  
}