HDU 1542 Atlantis （线段树求矩阵覆盖面积）

题意：给你n个矩阵求覆盖面积。

思路：看了别人的结题报告

给定一个矩形的左下角坐标和右上角坐标分别为:(x1,y1)、(x2,y2)，对这样的一个矩形，我们构造两条线段，一条定位在x1，它在y坐标的区间是[y1,y2]，并且给定一个cover域值为1；另一条线段定位在x2，区间一样是[y1,y2]，给定它一个cover值为-1。根据这样的方法对每个矩形都构造两个线段，最后将所有的线段根据所定位的x从左到右进行排序

#include <iostream>#include <stdio.h>#include <string.h>#include <algorithm>using namespace std;#define M 100#define inf 0x3fffffff#define maxn 500000*2struct seg{int flag;double up,down,x;}line[M*5];int cmp(seg a,seg b){return a.x<b.x;}struct Tree{double up,down,x;int cover;bool flag;}tree[M*M];double y[M*3];//对y进行分割，建树void build(int id,int l,int r){tree[id].down=y[l];tree[id].up=y[r];tree[id].flag=false;tree[id].cover=0;tree[id].x=-1;if(l+1==r){tree[id].flag=true;return ;}int mid=(l+r)/2;build(id*2,l,mid);build(id*2+1,mid,r);}double insert(int id,double x,double l,double r,int flag)//flag表示为左边还是右边{if(tree[id].down>=r||tree[id].up<=l) return 0;if(tree[id].flag){if(tree[id].cover>0)//递归到了叶子节点{double temp_x=tree[id].x;double ans=(x-temp_x)*(tree[id].up-tree[id].down);tree[id].x=x;//定位上一次的xtree[id].cover+=flag;return ans;}else {tree[id].cover+=flag;tree[id].x=x;return 0;}}double ans1,ans2;ans1=insert(id*2,x,l,r,flag);ans2=insert(id*2+1,x,l,r,flag);return ans1+ans2;}int main(){int t,ca=1;while(scanf("%d",&t)&&t){double x1,x2,y1,y2;int k=0;for(int i=0;i<t;i++){scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);k++;y[k]=y1;line[k].down=y1;line[k].up=y2;line[k].x=x1;line[k].flag=1;//1表示左边k++;y[k]=y2;line[k].x=x2;line[k].down=y1;line[k].up=y2;line[k].flag=-1;//-1表示右边}sort(y+1,y+k+1);sort(line+1,line+1+k,cmp);build(1,1,k);double ans=0;for(int i=1;i<=k;i++)ans+=insert(1,line[i].x,line[i].down,line[i].up,line[i].flag);printf("Test case #%d\nTotal explored area: %.2f\n\n",ca++,ans);}return 0;}/*210 10 20 2015 15 25 25.50*/