http://acm.hdu.edu.cn/showproblem.php?pid=1542
给出n个矩形的左下角和右上角的坐标,求矩形面积的并
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <queue>
#include <map>
#include <algorithm>
using namespace std;
const int maxn=1005;
int n;
double Y[maxn*2];
struct Line{
double x; //横坐标的值
double y1,y2;//纵坐标区间
int flag; //flag==1时表示他是长方形左边;flag==-1时,表示他是长方形右边
}line[maxn*2];
struct Node{
double x; //最新到过该节点的线对应的x
double left,right; //左边界、右边界
bool sign; //标记该节点是否是叶子节点
int cover; //覆盖值
}node[maxn<<3];
///这里不能使MAXN左移两位,因为有N个矩形,对应2*N个线段,则应该是2*N左移2位。
bool cmp(struct Line a,struct Line b){
return a.x<b.x;
}
void Build(int root,int l,int r){
node[root].left=Y[l];
node[root].right=Y[r];
node[root].x=-1;
node[root].sign=false; //说明该节点不是叶子节点
node[root].cover=0;
if (r-l==1){
node[root].sign=true; //是叶子节点
return;
}
int mid=(l+r)>>1;
Build(root<<1,l,mid);
Build(root<<1|1,mid,r);
}
///参数含义,root为节点,line_x代表当前线对应的横坐标,l下界,r上界
///flag代表当前的线是长方形的左侧的边还是右侧的边
double Calculate(int root,double line_x,double l,double r,int flag){
if (r<=node[root].left||l>=node[root].right){
return 0;
}
if (node[root].sign){ //如果是叶子节点
if (node[root].cover>0){
double pre=node[root].x;
double ans=(line_x-pre)*(node[root].right-node[root].left);
node[root].x=line_x;
node[root].cover+=flag;
return ans;
}
else{
node[root].x=line_x;
node[root].cover+=flag;
return 0;
}
}
double ans1=Calculate(root<<1,line_x,l,r,flag);
double ans2=Calculate(root<<1|1,line_x,l,r,flag);
return ans1+ans2;
}
int main(){
long long a,b,c,d;
while (cin >> n && n){
int num=0;
for (int i=0;i<n;i++){
cin >> a >> b >> c >> d;
Y[num]=b;
line[num].x=a;
line[num].y1=b;
line[num].y2=d;
line[num].flag=1;
num++;
Y[num]=d;
line[num].x=c;
line[num].y1=b;
line[num].y2=d;
line[num].flag=-1;
num++;
}
sort(Y,Y+num);
sort(line,line+num,cmp);
Build(1,0,num-1);
long long union_area=0;
for (int i=0;i<num;i++){
union_area+=Calculate(1,line[i].x,line[i].y1,line[i].y2,line[i].flag);
}
cout << union_area << endl;
}
cout << "*" << endl;
}