poj 1389 线段树+扫描线 求矩形交的面积

【题意】

求矩形交的面积

【题解】

相较之下，比求矩形交的周长简单不少。

维护这样几个量：

cov：线段被覆盖几次

sum：线段被覆盖的总长

【教训】

一开始想整段更新一定要用延迟标记，后来发现根本不用。

【代码】

#include <iostream>
#include <algorithm>
using namespace std;
const int maxl=50000,maxn=1000;

struct node
{
       int x,y1,y2,kk;
}a[maxn*2];
int s[maxl*4],cov[maxl*4];
int tot,ans,maxy,x1,x2,y1,y2;
void ins(int i,int l,int r,int x,int y,int k)
{
     if (x<=l && y>=r) cov[i]+=k;
     else 
     {
         int mid=(l+r)/2;
         if (x<mid) ins(i*2,l,mid,x,y,k);
         if (y>mid) ins(i*2+1,mid,r,x,y,k);
     }
     if (cov[i]>0) s[i]=r-l;
     else if (l+1==r) s[i]=0;
     else s[i]=s[i*2]+s[i*2+1];
}

bool cmp(node a,node b)
{
     return a.x<b.x;
}

int main()
{
    freopen("pin.txt","r",stdin);
    freopen("pou.txt","w",stdout);
    int i;
    while (scanf("%d%d%d%d",&x1,&y1,&x2,&y2)!=EOF)
    {
          if (x1==-1 && x2==-1 && y1==-1 && y2==-1) break;
          memset(cov,0,sizeof(cov));
          memset(s,0,sizeof(s));
          maxy=tot=0;
          while (x1!=-1 && x2!=-1 && y1!=-1 && y2!=-1)
          {
                maxy=max(maxy,y2);
                a[++tot].x=x1;a[tot].y1=y1;a[tot].y2=y2;a[tot].kk=1;
                a[++tot].x=x2;a[tot].y1=y1;a[tot].y2=y2;a[tot].kk=-1;
                scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
          }
          sort(a+1,a+tot+1,cmp);
          ans=0;//maxy+=25;
          ins(1,0,maxy,a[1].y1,a[1].y2,a[1].kk);
          for (i=2;i<=tot;i++)
          {
              ans+=s[1]*(a[i].x-a[i-1].x);
              ins(1,0,maxy,a[i].y1,a[i].y2,a[i].kk);
          }
          cout << ans << endl;
    }
    return 0;
}