Atlantis

原创 2013年12月05日 11:31:44

题目大意:

给出N个矩阵,矩阵会重叠,求矩阵围成的面积

SUCH AS

黑色部分就代表所求的面积。

题目思路:

线段树扫面线,画个图,我们可以把矩阵变成这样一个图形


按照
X轴进行排序,从左到右扫描,将各个不同颜色的长方形进行想家就可以得到答案

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 300
#define L(a) a<<1
#define R(a) a<<1|1
#define B(a) printf("***********************%d************************\n",a);

int n;

struct node
{
	int l,r;
	double y1,y2;
	double len;
	int zhuangtai;
}tree[N*4];

struct node2
{
	double x;
	double y1,y2;
	int flag;
	void in(double _x,double _y1,double _y2,int _flag)
	{
		x = _x;
		y1 = _y1;
		y2 = _y2;
		flag = _flag;
	}
}note[N];

double y[N*2];

int cmp(double a,double b)
{
	return a < b; 
}

bool cmp2(node2 a,node2 b)
{
	return a.x < b.x;
} 
void build_tree(int l,int r,int pos)
{
	tree[pos].l = l;
	tree[pos].r = r;
	tree[pos].y1 = y[l];
	tree[pos].y2 = y[r];
	tree[pos].zhuangtai = 0;
	tree[pos].len = 0;
	if( l + 1 == r )
		return ;
	int mid = ( l + r ) >> 1;
	build_tree(l,mid,L(pos));
	build_tree(mid,r,R(pos));
}

void get_len(int pos)
{
	if( tree[pos].zhuangtai > 0 )
		tree[pos].len = tree[pos].y2 - tree[pos].y1;
	else if ( tree[pos].r - tree[pos].l == 1 )
		tree[pos].len = 0;
	else
		tree[pos].len = tree[L(pos)].len + tree[R(pos)].len;
}

void updata(int pos,node2 a)
{
	//cout<<pos<<endl;
	//printf("a.y1 == %.2f a.y2 == %.2f\n",a.y1,a.y2);
	//printf("tree[%d].y1 == %.2f tree[%d].ye == %.2f",pos,tree[pos].y1,pos,tree[pos].y2);
	if( tree[pos].y1 == a.y1 && tree[pos].y2 == a.y2 )
	{
		tree[pos].zhuangtai += a.flag;
		get_len(pos);
		return ;
	}
	
	if( tree[L(pos)].y2 >= a.y2 )
		updata(L(pos),a);
	else if ( tree[R(pos)].y1 <= a.y1 )
		updata(R(pos),a);
	else
	{
		node2 tmp = a;
		tmp.y2 = tree[L(pos)].y2;
		updata(L(pos),tmp);
		tmp = a;
		tmp.y1 = tree[R(pos)].y1;
		updata(R(pos),tmp);
	}
	get_len(pos);
}

int main()
{
	int kas = 1;
	while(scanf("%d",&n)&&n)
	{
		double x1,x2,y1,y2;
		int t = 1;
		for(int i = 0 ; i < n ; i ++ )
		{
			scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
			note[t].in(x1,y1,y2,1);
			y[t ++] = y1;
			note[t].in(x2,y1,y2,-1);
			y[t ++] = y2;
		}
	//	B(1);
		sort(y+1,y+t,cmp);
		sort(note + 1,note+t,cmp2);
		build_tree(1,t-1,1);
	//	B(2);
	//	for(int i = 1 ; i < t ; i ++ )
	//		printf("y[%d] = %lf\n",i,y[i]);
	//	for(int i = 1 ; i < t ; i ++ )
	//		printf("note[%d].y1 = %lf note[%d].y2 = %lf note[%d].x =%lf\n",i,note[i].y1,i,note[i].y2,i,note[i].x);
		updata(1,note[1]);
		//B(3);
		double sum = 0.0;
		for(int i = 2 ; i < t ; i ++ )
		{
			sum += tree[1].len * ( note[i].x - note[i - 1].x);
			updata(1,note[i]);
		}
	//	B(4);
		printf("Test case #%d\n",kas++);  
        printf("Total explored area: %.2lf\n\n",sum); 
	}
	return 0;
}

相关文章推荐

Atlantis解题报告

  • 2012年10月08日 14:45
  • 53KB
  • 下载

The Rise Of Atlantis 含破解檔

  • 2010年11月09日 13:42
  • 19.6MB
  • 下载

Atlantis hdu 1542 线段树 扫面线 区间合并

http://acm.hdu.edu.cn/showproblem.php?pid=1542

Atlantis(P1151)

先用好像不是线段树的方法模仿做出,,然后再做一下线段做的方法吧 #include #include #include #include #include #include #include ...

hdu 1542 Atlantis(线段树&扫描线&面积并)

Atlantis Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total...

poj 1151 Atlantis(线段树 扫描线)

Atlantis Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 8759   Accep...
  • fp_hzq
  • fp_hzq
  • 2011年04月11日 18:22
  • 744

hdu 1542 Atlantis(矩形面积并)

Atlantis Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total...

hdu 1542 Atlantis 矩形面积并

// hdu 1542 Atlantis 矩形面积并 // // 解题思路: // // 将矩形沿着x轴或者y轴切割.然后求和 // 具体做法: // 将每个矩形看成是两条线段.将x方向 ...

hdu-1542-Atlantis 线段树

http://xuyemin520.is-programmer.com/posts/25988 http://hi.baidu.com/nicker2010/item/d88fe842c22163a1...

[POJ1151] Atlantis - 矩形切割

题目描述给出n个矩形,求这些矩形面积的并。输入格式输入包含多组测试数据; 每组数据的第一行包含一个整数n(1...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Atlantis
举报原因:
原因补充:

(最多只允许输入30个字)