题意：给你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*2
struct 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;//定位上一次的x
tree[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;
}
/*
2
10 10 20 20
15 15 25 25.5
0
*/
HDU 1542 Atlantis （线段树求矩阵覆盖面积）,布布扣,bubuko.com
HDU 1542 Atlantis （线段树求矩阵覆盖面积）
原文地址：http://blog.csdn.net/u012861385/article/details/35569373
