【集训Day3 离散化】矩形覆盖

矩形覆盖（planting）
【问题描述】
给定在一个平面坐标系上的N(1 <= N <= 100)个矩形区域，这N个矩形可能有相互覆盖的部分。求平面上被所有矩形覆盖的总面积，重复部分只算一次。
【输入格式】planting.in
第一行一个整数N。
接下来N行，每行包含4个整数x1 y1 x2 y2，表示一个矩形区域。矩形的左上角为(x1,y1)，右下角为(x2,y2)。坐标范围为–10^8到10^8之间的整数。
【输出格式】planting.out
输出被矩形覆盖的总面积。
【输入样例】
2
0 5 4 1
2 4 6 2
【输出样例】
20

【解题思路】
将输入的行和列用两个数组映射到1···n中去，将其离散化，利用另一个数组存下每一对离散的变量之间的差距
【解题反思】

• 解题时要注意细节

【参考程序】

#include<iostream>#include<cstdio>#include<algorithm>using namespace std;struct data{    int num;    int id;}x[205],y[205];bool comp(data a,data b){    return a.num<b.num;}int n,bx[205],by[205],c[205][205],ans,pic[205][205];int main(){    freopen("planting.in","r",stdin);    freopen("planting.out","w",stdout);    cin>>n;    for (int i=1;i<=n;i++)     {        cin>>x[2*i-1].num>>y[2*i].num>>x[2*i].num>>y[2*i-1].num;//由于是平面直角坐标系，列轴是向上的。         x[2*i-1].id=2*i-1;        y[2*i-1].id=2*i-1;        x[2*i].id=2*i;        y[2*i].id=2*i;    }    sort(x+1,x+1+2*n,comp);    sort(y+1,y+1+2*n,comp);    for (int i=1;i<=2*n;i++)     {        bx[x[i].id]=i;        by[y[i].id]=i;    }//离散化     for (int i=1;i<2*n;i++)        for (int j=1;j<2*n;j++)            c[i][j]=(x[i+1].num-x[i].num)*(y[j+1].num-y[j].num);//每一个离散的变量间的差距     for (int i=1;i<=n;i++)        for (int j=bx[2*i-1];j<bx[2*i];j++)            for (int k=by[2*i-1];k<by[2*i];k++) pic[j][k]=1;//传统的标记方法，简单暴力     for (int i=1;i<2*n;i++)        for (int j=1;j<2*n;j++)            if (pic[i][j]==1) ans+=c[i][j];//统计     cout<<ans;    return 0;}