扫描线 求多个矩形合并后的周长 hdu1828 Picture

题目链接   hdu1828 Picture 

题目大意：给你多个矩形，求他们合并后的周长，被覆盖的边不能算进周长之内。

解题思路： 其实周长并和面积并没什么很大的区别，只不过周长并增加了判断左右端点是否被覆盖的标记 lbd 和rbd 数组， 以及numseg 数组 记录连续区间段数。

numseg : 一根扫描线扫描过去，会记录有多少个连续的区间段，每个连续的区间段都有两条相等的竖边，而每次扫描过去竖边长度都相等。

#include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 using namespace std;
 
 #define lz 2*u,l,mid
 #define rz 2*u+1,mid+1,r
 const int maxn=22222;
 int sum[4*maxn];              ///记录被覆盖区间的长度
 int lbd[4*maxn], rbd[4*maxn]; ///记录左右端点是否被覆盖
 int numseg[4*maxn];            ///记录该区间连续的段数
 int flag[4*maxn];              /// 记录该区间是否被覆盖
 
 struct Node
 {
     int lx, rx, y, s;
     int lbd, rbd;
     Node() {};
     Node(int lx_, int rx_, int y_, int s_)
     {
         lx=lx_, rx=rx_, y=y_, s=s_;
     }
     bool operator <(const Node &S) const
     {
         if(y==S.y) return s>S.s;
         return y<S.y;
     }
 } line[maxn];
 
 void push_up(int u, int l, int r)
 {
     if(flag[u])
     {
         lbd[u]=1;
         rbd[u]=1;
         sum[u]=r-l+1;
         numseg[u]=2;
     }
     else if(l==r)
     {
         sum[u]=numseg[u]=lbd[u]=rbd[u]=0;
     }
     else
     {
         lbd[u]=lbd[2*u];
         rbd[u]=rbd[2*u+1];
         sum[u]=sum[2*u]+sum[2*u+1];
         numseg[u]=numseg[2*u]+numseg[2*u+1];
         if(rbd[2*u]&&lbd[2*u+1]) numseg[u]-=2;
     }
 }
 
 void Update(int u, int l, int r, int tl, int tr, int c)
 {
     if(tl<=l&&r<=tr)
     {
         flag[u]+=c;
         push_up(u,l,r);
         return ;
     }
     int mid=(l+r)>>1;
     if(tr<=mid) Update(lz,tl,tr,c);
     else if(tl>mid) Update(rz,tl,tr,c);
     else
     {
         Update(lz,tl,mid,c);
         Update(rz,mid+1,tr,c);
     }
     push_up(u,l,r);
 }
 
 int main()
 {
     int n;
     while(cin >> n)
     {
         int x1, x2, y1, y2;
         int num=0, lbd=10000, rbd=-10000;
         for(int i=0; i<n; i++)
         {
             scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
             line[++num]=Node(x1,x2,y1,1);
             line[++num]=Node(x1,x2,y2,-1);
             lbd=min(lbd,x1);
             rbd=max(rbd,x2);
         }
         sort(line+1,line+num+1);
         int ans=0, last=0;
         for(int i=1; i<=num; i++)
         {
             Update(1,lbd,rbd-1,line[i].lx,line[i].rx-1,line[i].s);
             ans+=numseg[1]*(line[i+1].y-line[i].y);
             ans+=abs(sum[1]-last);
             last=sum[1];
         }
         cout << ans <<endl;
     }
     return 0;
 }