HDU 1828 Picture (线段树:扫描线周长)

依然是扫描线，只不过是求所有矩形覆盖之后形成的图形的周长。

容易发现，扫描线中的某一条横边对答案的贡献。

其实就是 加上/去掉这条边之前的答案 和 加上/去掉这条边之后的答案 之差的绝对值

然后横着竖着都做一遍就行了

 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define N 10010
 #define ll long long
 using namespace std;

 int n,sz;
 int a[N],cnt[N<<],sum[N<<];
 struct SQU{
     int a1,b1,a2,b2;
 }q[N];
 struct node{
     int l,r,la,ra,h,f;
 }sc[N];
 int cmp1(node s1,node s2){
     if(s1.h!=s2.h) return s1.h<s2.h;
     else return s1.f>s2.f;
 }
 void pushup(int l,int r,int rt)
 {
     if(cnt[rt]>) sum[rt]=a[r+]-a[l];
     else if(l==r) sum[rt]=;
     else sum[rt]=sum[rt<<]+sum[rt<<|];
 }
 void update(int L,int R,int l,int r,int rt,int w)
 {
     if(L<=l&&r<=R)
     {
         cnt[rt]+=w;
         pushup(l,r,rt);
         return;
     }
     int mid=(l+r)>>;
     if(L<=mid) update(L,R,l,mid,rt<<,w);
     if(R>mid) update(L,R,mid+,r,rt<<|,w);
     pushup(l,r,rt);
 }
 void clr()
 {
     memset(a,,sizeof(a));
     memset(sc,,sizeof(sc));
     memset(cnt,,sizeof(cnt));
     memset(sum,,sizeof(sum));
 }
 int solvex()
 {
     int ans=,lst=;
     for(int i=;i<=n;i++)
     {
         a[*i-]=q[i].a1,a[*i]=q[i].a2;
         sc[*i-].l=q[i].a1,sc[*i].l=q[i].a1;
         sc[*i-].r=q[i].a2,sc[*i].r=q[i].a2;
         sc[*i-].h=q[i].b1,sc[*i].h=q[i].b2;
         sc[*i-].f=,sc[*i].f=-;
     }
     sort(a+,a+*n+);
     sz=unique(a+,a+*n+)-(a+);
     sort(sc+,sc+*n+,cmp1);
     for(int i=;i<=*n;i++)
     {
         int la=lower_bound(a+,a+sz+,sc[i].l)-a;
         int ra=lower_bound(a+,a+sz+,sc[i].r)-a;
         lst=sum[];
         update(la,ra-,,sz,,sc[i].f);
         ans+=abs(sum[]-lst);
     }
     return ans;
 }
 int solvey()
 {
     int ans=,lst=;
     for(int i=;i<=n;i++)
     {
         a[*i-]=q[i].b1,a[*i]=q[i].b2;
         sc[*i-].l=q[i].b1,sc[*i].l=q[i].b1;
         sc[*i-].r=q[i].b2,sc[*i].r=q[i].b2;
         sc[*i-].h=q[i].a1,sc[*i].h=q[i].a2;
         sc[*i-].f=,sc[*i].f=-;
     }
     sort(a+,a+*n+);
     sz=unique(a+,a+*n+)-(a+);
     sort(sc+,sc+*n+,cmp1);
     for(int i=;i<=*n;i++)
     {
         int la=lower_bound(a+,a+sz+,sc[i].l)-a;
         int ra=lower_bound(a+,a+sz+,sc[i].r)-a;
         lst=sum[];
         update(la,ra-,,sz,,sc[i].f);
         ans+=abs(sum[]-lst);
     }
     return ans;
 }

 int main()
 {
     //freopen("data.in","r",stdin);
     scanf("%d",&n);
     for(int i=;i<=n;i++)
         scanf("%d%d%d%d",&q[i].a1,&q[i].b1,&q[i].a2,&q[i].b2);
     int ret=;
     ret+=solvex();
     clr();
     ret+=solvey();
     printf("%d\n",ret);
     return ;
 }