Codevs P3044矩形面积求并

Codevs P3044矩形面积求并

---

题目描述

Description

输入n个矩形，求他们总共占地面积（也就是求一下面积的并）

输入描述 Input Description
可能有多组数据，读到n=0为止(不超过15组)

每组数据第一行一个数n，表示矩形个数(n<=100)

接下来n行每行4个实数x1,y1,x2,y1(0 <= x1 < x2 <= 100000;0 <= y1 < y2 <= 100000)，表示矩形的左下角坐标和右上角坐标

输出描述 Output Description

每组数据输出一行表示答案

样例输入 Sample Input

2
10 10 20 20
15 15 25 25.5
0

样例输出 Sample Output

180.00

数据范围及提示 Data Size & Hint

无

---

题解

（其实，这题目是我在某本将结构题的书中看到的……，数据比这个大了100倍，QAQ）

线段树

把矩形看做两条平行的线段（可以与x轴平行也可以与y轴平行，这里我的代码是与x轴平行的），然后，我们只需要通过离散将线段的纵坐标排成有序排列后，按纵坐标从小到大做线段树的线段覆盖（一个矩形所分离出的两条线段中，上面的线段的赋值为-1，下面的线段的赋值为1）得出覆盖的总长度为$l$，并求出这条线段到下一条线段的距离$d$，然后只要把$l$*$d$累计入$ans$即可

（$ans$忘清0，WA了几发……尴尬……）

---

代码

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;

int n;
int c[805];
double yy[805],sum[805];
struct data
{
    double l,r,y;
    int flag;
} tr[805];

bool cmp(data a,data b) {return a.y<b.y;}

int find(double x)
{
    int l=1,r=2*n,mid;
    while (l<=r)
    {
        mid=(l+r)>>1;
        if (yy[mid]<x) l=mid+1;
        else if (yy[mid]==x) return mid;
        else r=mid-1;
    }
}

void pushup(int l,int r,int rt)
{
    if (c[rt]) sum[rt]=yy[r+1]-yy[l];
    else if (l==r) sum[rt]=0;
    else sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}

void updata(int ll,int rr,int l,int r,int add,int rt)
{
    if (ll==l&&rr==r) {c[rt]+=add;pushup(ll,rr,rt);return;}
    int mid=(ll+rr)>>1;
    if (r<=mid) updata(ll,mid,l,r,add,rt<<1);
    else if (mid<l) updata(mid+1,rr,l,r,add,rt<<1|1);
    else {updata(ll,mid,l,mid,add,rt<<1);updata(mid+1,rr,mid+1,r,add,rt<<1|1);}
    pushup(ll,rr,rt);
}

int main()
{
    while (true)
    {
        scanf("%d",&n);
        if (n==0) break;
        memset(c,0,sizeof(c));
        memset(sum,0,sizeof(sum));
        double x1,x2,y1,y2;
        for (int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            tr[2*i-1].l=x1;tr[2*i-1].r=x2;tr[2*i-1].y=y1;
            tr[2*i].l=x1;tr[2*i].r=x2;tr[2*i].y=y2;
            tr[2*i-1].flag=1;tr[2*i].flag=-1;
            yy[2*i-1]=x1;yy[2*i]=x2;
        }
        sort(yy+1,yy+2*n+1);
        sort(tr+1,tr+2*n+1,cmp);
        double ans=0;
        for (int i=1;i<2*n;i++)
        {
            int l=find(tr[i].l),r=find(tr[i].r)-1;
            if (l<=r) updata(1,2*n,l,r,tr[i].flag,1);
            ans+=sum[1]*(tr[i+1].y-tr[i].y);
        }
        printf("%.2lf\n",ans);
    }
    return 0;
}