POJ 1151 矩形面积并

题解：
把横坐标离散化，然后把矩形拆成上下两条边，上边标成-1，下边标成1，然后按情况增减面积即可~
View Code
 1 #include <iostream>
 2 #include <cstring>
 3 #include <cstdio>
 4 #include <cstdlib>
 5 #include <algorithm>
 6 
 7 #define N 222
 8 
 9 using namespace std;
10 
11 struct LI
12 {
13     double y,sx,tx;
14     int nsx,ntx,fg;
15 }li[N<<1];
16 
17 int gs,cs,cnt,n,cas;
18 int c[N<<1];
19 double a[N<<1],ans;
20 
21 inline bool cmp(const LI &a,const LI &b)
22 {
23     return a.y<b.y;
24 }
25 
26 inline void read()
27 {
28     gs=cs=0;double tx,ty,sx,sy;
29     for(int i=1;i<=n;i++)
30     {
31         scanf("%lf%lf%lf%lf",&sx,&sy,&tx,&ty);
32         a[++gs]=sx; a[++gs]=tx;
33         li[++cs].y=sy; li[cs].sx=sx; li[cs].tx=tx; li[cs].fg=-1;
34         li[++cs].y=ty; li[cs].sx=sx; li[cs].tx=tx; li[cs].fg=1;
35     }
36     sort(li+1,li+1+cs,cmp);
37     sort(a+1,a+1+gs);
38     cnt=unique(a+1,a+1+gs)-a-1;
39     for(int i=1;i<=cs;i++)
40     {
41         li[i].nsx=lower_bound(a+1,a+1+cnt,li[i].sx)-a;
42         li[i].ntx=lower_bound(a+1,a+1+cnt,li[i].tx)-a;
43     }
44 }
45 
46 inline void go()
47 {
48     memset(c,0,sizeof c);
49     ans=0.0;
50     for(int i=1;i<=cs;i++)
51         for(int j=li[i].nsx;j<li[i].ntx;j++)
52         {
53             if(li[i].fg==-1)
54             {
55                 c[j]--;
56                 if(c[j]==-1) ans-=(a[j+1]-a[j])*li[i].y;
57             }
58             else
59             {
60                 c[j]++;
61                 if(c[j]==0) ans+=(a[j+1]-a[j])*li[i].y;
62             }
63         }
64     printf("Test case #%d\n",++cas);
65     printf("Total explored area: %.2lf\n\n",ans);
66 }
67 
68 int main()
69 {
70     while(scanf("%d",&n),n) read(),go();
71     return 0;
72 }