HDU 1255 覆盖的面积 （扫描线 线段树 离散化 矩形面积并）

题目链接

题意：中文题意。

分析：纯手敲，与上一道题目很相似，但是刚开始我以为只是把cnt》=0改成cnt>=2就行了，、

但是后来发现当当前加入的线段的范围之前 还有线段的时候就不行了，因为虽然现在都不等于

2，但是之前的那个线段加上现在的已经覆盖2次了。

 #include <iostream>
 #include <cstdio>
 #include <vector>
 #include <cstring>
 #include <cstdlib>
 #include <algorithm>
 #define LL __int64
 #define lson l, mid, 2*rt
 #define rson mid+1, r, 2*rt+1
 const int maxn = +;
 using namespace std;
 int n;
 double y[maxn];
 struct node
 {
     int l, r, c;
     double cnt, lf, rf, more; //cnt还是代表覆盖的长度，增加了more代表两次及以上覆盖的长度
 }tr[*maxn];
 struct Line
 {
     double x, y1, y2;
     int f;
 }line[maxn];
 bool cmp(Line a, Line b)
 {
     return a.x < b.x;
 }
 void build(int l, int r, int rt)
 {
     tr[rt].l = l; tr[rt].r = r;
     tr[rt].c = ;  tr[rt].cnt = ;
     tr[rt].more = ;
     tr[rt].rf = y[r];  tr[rt].lf = y[l];
     if(l+==r) return;
     int mid = (l+r)/;
     build(l, mid, *rt);
     build(mid, r, *rt+);
 }
 void calen(int rt)
 {
     if(tr[rt].c==)
     {
         if(tr[rt].l+==tr[rt].r)
         {
             tr[rt].cnt = ; tr[rt].more = ;
         }
         else
         {
             tr[rt].cnt = tr[*rt].cnt+tr[*rt+].cnt;
             tr[rt].more = tr[*rt].more+tr[*rt+].more;
         }
     }
     if(tr[rt].c==)  //注意这一步是关键
     {
         tr[rt].cnt = tr[rt].rf-tr[rt].lf;
         if(tr[rt].l+==tr[rt].r)  //因为没有注意是否到最后，错了一遍
         tr[rt].more = ;
         else
         tr[rt].more = tr[*rt].cnt + tr[*rt+].cnt; //为1的时候如果下面也有就加上
     }
     if(tr[rt].c>=)
     {
         tr[rt].more = tr[rt].rf-tr[rt].lf;
         tr[rt].cnt = tr[rt].more;
     }
 }
 void update(int rt, Line e)
 {
     if(e.y1==tr[rt].lf && e.y2==tr[rt].rf)
     {
         tr[rt].c += e.f;
         calen(rt);
         return;
     }
     if(e.y2<=tr[*rt].rf) update(*rt, e);
     else if(e.y1>=tr[*rt+].lf) update(*rt+, e);
     else
     {
         Line tmp;
         tmp = e;
         tmp.y2 = tr[*rt].rf; update(*rt, tmp);
         tmp = e;
         tmp.y1 = tr[*rt+].lf; update(*rt+, tmp);
     }
     calen(rt);
 }
 int main()
 {
     int t, i, cnt;
     double x1, x2, y1, y2, ans;
     scanf("%d", &t);
     while(t--)
     {
         scanf("%d", &n);
         cnt = ; ans = ;
         for(i = ; i <= n; i++)
         {
             scanf("%lf%lf%lf%lf", &x1, &y1, &x2, &y2);
             line[cnt].x = x1; line[cnt].y1 = y1;
             line[cnt].y2 = y2; line[cnt].f = ;
             y[cnt++] = y1;
             line[cnt].x = x2; line[cnt].y1 = y1;
             line[cnt].y2 = y2; line[cnt].f = -;
             y[cnt++] = y2;
         }
         sort(y+, y+cnt);
         sort(line+, line+cnt, cmp);
         cnt --;
         build(, cnt, );
         update(, line[]);
         for(i = ; i <= cnt; i++)
         {
             ans += tr[].more*(line[i].x - line[i-].x);
             update(, line[i]);
         }
         printf("%.2lf\n", ans);
     }
     return ;
 }