线段树学习之：用线段数统计矩形的周长hdu(1828)

    学习如何用线段树来统计矩形的并面积后，接下来我学习一下如何用线段树来统计矩形周长的并！尽管其基本思路跟统计矩形的面积差不多，但估计是因为没有完全明白线段树的思想，我还是没有想出如何用线段树统计线段树的面积！

   在学习的过程中参考了以下两个博客的代码，最终以博客二的代码作为研究对象！

博客一：http://www.notonlysuccess.com/index.php/divide-tree/

博客二：http://hi.baidu.com/lw0090s/item/159f6ad6cab61b9c260ae731

   个人觉得博客二的代码更适合于像我这样的菜鸟！更容易理解！把代码贴出来！(注代码来源于博客二，不是本人写出来的)

1. #include<iostream> 
2. #include<algorithm> 
3. #include<stdio.h> 
4. using namespace std; 
5. struct Line 
6. { 
7.     int x,y1,y2; 
8.     int flag; 
9. } line[5005*2]; 
10. int y[5005*2]; 
11. struct tree 
12. { 
13.     int lines,len,l,r,c,lb,rb; //lb表示线段左边端点是否被覆盖为0，1，rb同理右边 
14. //lines表示当前线段下面覆盖的间隔的元线段的条数，len表示覆盖的长度 
15. } tree[5005*2*5]; 
16. int cmp(Line a,Line b) 
17. { 
18.     if(a.x!=b.x) 
19.         return a.x<b.x; 
20.     else 
21.         return a.flag>b.flag; 
22. } 
23. /* 
24.  建树不用说，很直观明白！就是以元线段来建树的 而且用到的离散化 
25. */ 
26. void build(int l,int r,int num) 
27. { 
28.     tree[num].l=l; 
29.     tree[num].r=r; 
30.     tree[num].lines=0; 
31.     tree[num].len=0; 
32.     tree[num].c=0; 
33.     tree[num].lb=0; 
34.     tree[num].rb=0; 
35.     if(l==r-1) 
36.         return; 
37.     int mid=(l+r)/2; 
38.     build(l,mid,2*num); 
39.     build(mid,r,2*num+1); 
40. } 
41. /* 
42.  个人感觉这里写得特别清晰，计算长度也很简单 
43. */ 
44. void calen(int num) 
45. { 
46.     if(tree[num].c>0) 
47.         tree[num].len=y[tree[num].r]-y[tree[num].l];//这里典型的以下标来建树 
48.     else 
49.     { 
50.         if(tree[num].l==tree[num].r-1)//元线段 
51.             tree[num].len=0; 
52.         else 
53.             tree[num].len=tree[num*2].len+tree[num*2+1].len; 
54.     } 
55.     return; 
56. } 
57. /* 
58.   这里统计线段的条数不好理解，我看第一种代码的时候就没有看懂 lb与rb的意思！ 
59.   在这里依然有点朦胧，感觉没有理解透！ 
60. */ 
61. void caline(int num) 
62. { 
63.     tree[num].lines=0; 
64.     tree[num].lb=0; 
65.     tree[num].rb=0; 
66.     if(tree[num].c>0) 
67.     { 
68.         tree[num].lines=2; 
69.         tree[num].rb=1; 
70.         tree[num].lb=1; 
71.     } 
72.     else 
73.     { 
74.         if(tree[num].l==tree[num].r-1)//元线段 
75.             tree[num].lines=0; 
76.         else 
77.         { 
78.             tree[num].lines=tree[num*2].lines+tree[2*num+1].lines-2*(tree[num*2].rb&tree[num*2+1].lb); 
79.             //这里就是lb和rb标记的原因，在合并时使用 
80.             tree[num].lb=tree[num*2].lb; 
81.             tree[num].rb=tree[num*2+1].rb; 
82.         } 
83.     } 
84.     return; 
85. } 
86. void update(Line e,int num) 
87. { 
88.     if(e.y1==y[tree[num].l]&&e.y2==y[tree[num].r])//与交和并的区别是这里插入到底了， 
89. //因为这里要计算下面的元线段得条数，所以必须插入到底 
90.     { 
91.         tree[num].c+=e.flag; 
92.         tree[num].lb=1; 
93.         tree[num].rb=1;//标记当前线段的端点的覆盖情况 
94.     } 
95.     else if(e.y2<=y[tree[num*2].r]) 
96.         update(e,num*2); 
97.     else if(e.y1>=y[tree[num*2+1].l]) 
98.         update(e,num*2+1); 
99.     else 
100.     { 
101.         Line tmp=e; 
102.         int mid=(tree[num].l+tree[num].r)/2; 
103.         tmp.y2=y[mid]; 
104.         update(tmp,num*2); 
105.         tmp=e; 
106.         tmp.y1=y[mid]; 
107.         update(tmp,num*2+1); 
108.     } 
109.     calen(num);//计算长度 
110.     caline(num); //计算元线段的条数 
111. } 
112. int main() 
113. { 
114.     int i,j,n; 
115.     int x1,y1,x2,y2,t; 
116.     int lenth,lines,ans; 
117.     while(scanf("%d",&n)!=EOF) 
118.     { 
119.         t=1; 
120.         for(i=1; i<=n; i++) 
121.         { 
122.             scanf("%d%d%d%d",&x1,&y1,&x2,&y2); 
123.             line[t].x=x1; 
124.             line[t].y1=y1; 
125.             line[t].y2=y2; 
126.             line[t].flag=1;//建立标记要看x的坐标大小 
127.             line[t+1].x=x2; 
128.             line[t+1].y1=y1; 
129.             line[t+1].y2=y2; 
130.             line[t+1].flag=-1; 
131.             y[t]=y1; 
132.             y[t+1]=y2; 
133.             t=t+2; 
134.         } 
135.         sort(line+1,line+t,cmp); 
136.         sort(y+1,y+t); 
137.         int l=1; 
138.         for(i=2; i<t; i++)//线段的离散化 
139.         { 
140.             if(y[i]!=y[i-1]) 
141.                 y[++l]=y[i];//这里因为要看元线段的条数所以将相同的点去掉 
142.         } 
143.         build(1,l,1); 
144.         lenth=0; 
145.         lines=0; 
146.         ans=0; 
147.         for(i=1; i<t; i++) 
148.         { 
149.             update(line[i],1); 
150.             ans+=abs(lenth-tree[1].len)+lines*(line[i].x-line[i-1].x); 
151.             lenth=tree[1].len; 
152.             lines=tree[1].lines; 
153.         } 
154.         printf("%d\n",ans); 
155.  
156.     } 
157.   //  system("pause"); 
158. } 

先把代码收藏着吧！以后学习的时候慢慢来参悟！