ACM 数据结构 线段树 HDU 1540 Tunnel Warfare

1. Problem Description
   During the War of Resistance Against Japan, tunnel warfare was carried out extensively in the vast areas of north China Plain. Generally speaking, villages connected by tunnels lay in a line. Except the two at the ends, every village was directly connected with two neighboring ones.
   
   Frequently the invaders launched attack on some of the villages and destroyed the parts of tunnels in them. The Eighth Route Army commanders requested the latest connection state of the tunnels and villages. If some villages are severely isolated, restoration of connection must be done immediately!
   
    
   
   
   Input
   The first line of the input contains two positive integers n and m (n, m ≤ 50,000) indicating the number of villages and events. Each of the next m lines describes an event.
   
   There are three different events described in different format shown below:
   
   D x: The x-th village was destroyed.
   
   Q x: The Army commands requested the number of villages that x-th village was directly or indirectly connected with including itself.
   
   R: The village destroyed last was rebuilt.
   
    
   
   
   Output
   Output the answer to each of the Army commanders’ request in order on a separate line.
   
    
   
   
   Sample Input

       
       
   
        
        7 9
   D 3
   D 6
   D 5
   Q 4
   Q 5
   R
   Q 4
   R
   Q 4
       
       

    
   
   
   Sample Output

       
       
   
        
        1
   0
   2
   4
       
       

   
   
2. 
   
3. #include <stdio.h>  
4. #include <string.h>  
5. #include <algorithm>  
6. #include <math.h>  
7. #include <stdlib.h>  
8. using namespace std;  
9.   
10. const int maxn = 50000+10;  
11.   
12. int n,m;  
13. int s[maxn],top;//s为模拟栈  
14.   
15. struct node  
16. {  
17.     int l,r;  
18.     int ls,rs,ms;//ls,左端最大连续区间，rs右端最大连续区间，ms区间内最大连续区间  
19. } a[maxn<<2];  
20.   
21. void init(int l,int r,int i)  
22. {  
23.     a[i].l = l;  
24.     a[i].r = r;  
25.     a[i].ls = a[i].rs = a[i].ms = r-l+1;  
26.     if(l!=r)  
27.     {  
28.         int mid = (l+r)>>1;  
29.         init(l,mid,i*2);  
30.         init(mid+1,r,2*i+1);  
31.     }  
32. }  
33.   
34. void insert(int i,int t,int x)  
35. {  
36.     if(a[i].l == a[i].r)  
37.     {  
38.         if(x==1)  
39.             a[i].ls = a[i].rs = a[i].ms = 1;//修复  
40.         else  
41.             a[i].ls = a[i].rs = a[i].ms = 0;//破坏  
42.         return ;  
43.     }  
44.     int mid = (a[i].l+a[i].r)>>1;  
45.     if(t<=mid)  
46.         insert(2*i,t,x);  
47.     else  
48.         insert(2*i+1,t,x);  
49.     a[i].ls = a[2*i].ls;//左区间  
50.     a[i].rs = a[2*i+1].rs;//右区间  
51.     a[i].ms = max(max(a[2*i].ms,a[2*i+1].ms),a[2*i].rs+a[2*i+1].ls);//父亲区间内的最大区间必定是，左子树最大区间，右子树最大区间，左右子树合并的中间区间，三者中最大的区间值  
52.     if(a[2*i].ls == a[2*i].r-a[2*i].l+1)//左子树区间满了的话，父亲左区间要加上右孩子的左区间  
53.         a[i].ls += a[2*i+1].ls;  
54.     if(a[2*i+1].rs == a[2*i+1].r-a[2*i+1].l+1)//同理  
55.         a[i].rs += a[2*i].rs;  
56. }  
57.   
58. int query(int i,int t)  
59. {  
60.     if(a[i].l == a[i].r || a[i].ms == 0 || a[i].ms == a[i].r-a[i].l+1)//到了叶子节点或者该访问区间为空或者已满都不必要往下走了  
61.         return a[i].ms;  
62.     int mid = (a[i].l+a[i].r)>>1;  
63.     if(t<=mid)  
64.     {  
65.         if(t>=a[2*i].r-a[2*i].rs+1)//因为t<=mid，看左子树，a[2*i].r-a[2*i].rs+1代表左子树右边连续区间的左边界值，如果t在左子树的右区间内，则要看右子树的左区间有多长并返回  
66.             return query(2*i,t)+query(2*i+1,mid+1);  
67.         else  
68.             return query(2*i,t);//如果不在左子树的右边界区间内，则只需要看左子树  
69.     }  
70.     else  
71.     {  
72.         if(t<=a[2*i+1].l+a[2*i+1].ls-1)//同理  
73.             return query(2*i+1,t)+query(2*i,mid);  
74.         else  
75.             return query(2*i+1,t);  
76.     }  
77. }  
78.   
79. int main()  
80. {  
81.     int i,j,x;  
82.     char ch[2];  
83.     while(~scanf("%d%d",&n,&m))  
84.     {  
85.         top = 0;  
86.         init(1,n,1);  
87.         while(m--)  
88.         {  
89.             scanf("%s",ch);  
90.             if(ch[0] == 'D')  
91.             {  
92.                 scanf("%d",&x);  
93.                 s[top++] = x;  
94.                 insert(1,x,0);  
95.             }  
96.             else if(ch[0] == 'Q')  
97.             {  
98.                 scanf("%d",&x);  
99.                 printf("%d\n",query(1,x));  
100.             }  
101.             else  
102.             {  
103.                 if(x>0)  
104.                 {  
105.                     x = s[--top];  
106.                     insert(1,x,1);  
107.                 }  
108.             }  
109.         }  
110.     }  
111.   
112.     return 0;  
113. }