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

题目链接:https://cn.vjudge.net/problem/HDU-1542

题目大意:求面积并

具体思路:我们首先把矩形分割成一横条一横条的,然后对于每一个我们给定的矩形,我们将储存两个点,一个是左下角,一个是右上角。对于左下角的点,我们将他标记为-1,对于右上角的点,我们把它标记成1,-1代表这块区域的面积是需要减去的,1代表这块区域的面积是需要加上的,然后我们通过扫描线的形式,从y轴从下往上扫就可以了。

离散化的过程,数组还是尽量从0开始吧,因为在更新的过程中会出现0的情况,如果下标从1开始存储的话,会mle的,,,无限递归?

ps:这个线段树的区间分开的时候和之前的不太一样,这个线段树{1,5}分开的是{1,3}和{4,5}

注意这个线段树每一个节点代表的是区间,节点1-3代表的是区间1-4,也就是 节点1代表 区间1-2 ,然后我们查询1-4 的时候,直接查询节点1-3就可以了,更新区间1-3的时候,更新节点1-2就可以了。

AC代码:

 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 }

 

posted @ 2019-01-29 00:05  Let_Life_Stop  阅读(186)  评论(0编辑  收藏  举报