线段树->面积并 Atlantis HDU - 1542

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <string>
 4 #include <cstring>
 5 #include <algorithm>
 6 #include <cmath>
 7 #include <vector>
 8 #include <queue>
 9 #include <stack>
10 #include <ctime>
11 #define ll long long
12 # define lson l,m,rt<<1
13 # define rson m+1,r,rt<<1|1
14 using namespace std;
15 const int maxn = 600+100;
16 # define inf 0x3f3f3f3f
17 struct node
18 {
19     double l,r,h;
20     int d;
21     node() {}
22     node(double xx,double yy,double zz,int t)
23     {
24         l=xx;
25         r=yy;
26        h=zz;
27         d=t;
28     }
29     bool friend operator < (node t1,node t2)
30     {
31         return t1.h<t2.h;
32     }
33 } q[maxn];
34 double tree[maxn<<2];
35 double Hash[maxn<<2];
36 int lazy[maxn<<2];
37 void up(int l,int r,int rt)
38 {
39     if(lazy[rt])
40         tree[rt]=Hash[r+1]-Hash[l];//注意这个地方，我们需要计算的是区域，而不是一个点。
41     else if(l==r)
42         tree[rt]=0;
43     else
44         tree[rt]=tree[rt<<1]+tree[rt<<1|1];
45 }
46 void update(int l,int r,int rt,int L,int R,int pos)
47 {
48     if(L<=l&&R>=r)
49     {
50         lazy[rt]+=pos;
51         up(l,r,rt);
52         return ;
53     }
54     int m=(l+r)>>1;
55     if(L<=m)
56         update(lson,L,R,pos);
57     if(R>m)
58         update(rson,L,R,pos);
59     up(l,r,rt);
60 }
61 int main()
62 {
63     int n,Case=0;
64     while(scanf("%d",&n)&&n)
65     {
66         int num=0;
67         double x1,y1,x2,y2;
68         for(int i=1; i<=n; i++)
69         {
70             scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2);
71             q[num]= {x1,x2,y1,1};
72             Hash[num++]=x1;
73             q[num]= {x1,x2,y2,-1};
74             Hash[num++]=x2;
75         }
76         sort(Hash,Hash+num);
77         sort(q,q+num);
78         int m=1;
79         for(int i=1; i<num; i++)
80         {
81             if(Hash[i-1]==Hash[i])
82                 continue;
83             Hash[m++]=Hash[i];
84         }
85         double ans=0;
86         memset(tree,0,sizeof(tree));
87         memset(lazy,0,sizeof(lazy));
88         for(int i=0; i<num; i++)
89         {
90             int l=lower_bound(Hash,Hash+m,q[i].l)-Hash;
91             int r=lower_bound(Hash,Hash+m,q[i].r)-Hash-1;
92             update(0,m-1,1,l,r,q[i].d);
93             ans+=tree[1]*(q[i+1].h-q[i].h);
94         }
95         printf("Test case #%d\n",++Case);
96         printf("Total explored area: %.2lf\n\n",ans);
97     }
98     return 0;
99 }