poj_1389 Area of Simple Polygons 线段树+扫描线
来源:互联网 发布:淘宝网站客服电话 编辑:程序博客网 时间:2024/06/06 09:54
Area of Simple Polygons
Time Limit: 1000MS Memory Limit: 65536KTotal Submissions: 3139 Accepted: 1601
Description
There are N, 1 <= N <= 1,000 rectangles in the 2-D xy-plane. The four sides of a rectangle are horizontal or vertical line segments. Rectangles are defined by their lower-left and upper-right corner points. Each corner point is a pair of two nonnegative integers in the range of 0 through 50,000 indicating its x and y coordinates.
Assume that the contour of their union is defi ned by a set S of segments. We can use a subset of S to construct simple polygon(s). Please report the total area of the polygon(s) constructed by the subset of S. The area should be as large as possible. In a 2-D xy-plane, a polygon is defined by a finite set of segments such that every segment extreme (or endpoint) is shared by exactly two edges and no subsets of edges has the same property. The segments are edges and their extremes are the vertices of the polygon. A polygon is simple if there is no pair of nonconsecutive edges sharing a point.
Example: Consider the following three rectangles:
rectangle 1: < (0, 0) (4, 4) >,
rectangle 2: < (1, 1) (5, 2) >,
rectangle 3: < (1, 1) (2, 5) >.
The total area of all simple polygons constructed by these rectangles is 18.
Assume that the contour of their union is defi ned by a set S of segments. We can use a subset of S to construct simple polygon(s). Please report the total area of the polygon(s) constructed by the subset of S. The area should be as large as possible. In a 2-D xy-plane, a polygon is defined by a finite set of segments such that every segment extreme (or endpoint) is shared by exactly two edges and no subsets of edges has the same property. The segments are edges and their extremes are the vertices of the polygon. A polygon is simple if there is no pair of nonconsecutive edges sharing a point.
Example: Consider the following three rectangles:
rectangle 1: < (0, 0) (4, 4) >,
rectangle 2: < (1, 1) (5, 2) >,
rectangle 3: < (1, 1) (2, 5) >.
The total area of all simple polygons constructed by these rectangles is 18.
Input
The input consists of multiple test cases. A line of 4 -1's separates each test case. An extra line of 4 -1's marks the end of the input. In each test case, the rectangles are given one by one in a line. In each line for a rectangle, 4 non-negative integers are given. The first two are the x and y coordinates of the lower-left corner. The next two are the x and y coordinates of the upper-right corner.
Output
For each test case, output the total area of all simple polygons in a line.
Sample Input
0 0 4 41 1 5 21 1 2 5-1 -1 -1 -10 0 2 21 1 3 32 2 4 4-1 -1 -1 -1-1 -1 -1 -1
Sample Output
1810
#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>using namespace std;#define INF 0x3f3f3f3f#define inf -0x3f3f3f3f#define FOR(a,b) for(int i = a ; i < b ; i++)#define mem0(a) memset(a,0,sizeof(a))const int N = 50000+10;struct Node{ int x,y1,y2; int flag ;} node[N];//把一段段平行于y轴的线段表示成数组 ,//x是线段的x坐标,y1,y2线段对应的下端点和上端点的坐标//一个矩形 ,左边的那条边f为1,右边的为-1,//用来记录重叠情况,可以根据这个来计算,nod节点中的cbool cmp(Node a,Node b){ return a.x - b.x <0.0000001;}struct node{ int l,r;//线段树左右的整点 int ml, mr;//分别对应左右真实的浮点数端点 int s, len ;//s用来记录重叠情况,len计算实际长度} a[N*3];int y[N];//记录y坐标的数组void build(int i,int l,int r){ a[i].l = l ; a[i].r = r; a[i].ml = y[l]; a[i].mr = y[r]; a[i].s = 0 ; a[i].len = 0; if(a[i].l + 1 == a[i].r) { return ; } int mid =(l+r)>>1; build(i<<1,l,mid); build(i<<1|1,mid,r);}void callen(int i ){ if(a[i].s > 0 ) { a[i].len = a[i].mr -a[i].ml; } else if(a[i].r - a[i].l == 1) { a[i].len = 0 ; } else { a[i].len = a[i<<1].len+a[i<<1|1].len; } return ;}void update(int i,Node b) //加入线段b,从第一条开始,更新线段树{ if(a[i].ml == b.y1 && a[i].mr == b.y2) { a[i].s += b.flag ; callen(i); return ; } if(b.y2 <= a[i<<1].mr ) update(i<<1,b); else if(b.y1 >= a[i<<1|1].ml) update(i<<1|1,b); else { Node temp = b; temp.y2 = a[i<<1].mr; update(i<<1,temp); temp =b ; temp.y1 = a[i<<1|1].ml; update(i<<1|1,temp); } callen(i); return ;}int main(){ int x1,y1,x2,y2; int t= 1,flag=1; while(scanf("%d%d%d%d",&x1,&y1,&x2,&y2)!=EOF) { if(x1==-1&&x2==-1&&y1==-1&&y2==-1) { flag = 0; sort(node+1,node+t,cmp); sort(y+1,y+t); build(1,1,t-1); update(1,node[1]);//先加入第一条线段 int sum = 0 ; for(int i = 2 ; i < t ; i++) { sum += a[1].len * (node[i].x - node[i-1].x); update(1,node[i]);//从第二条开始加入 } printf("%d\n",sum); t =1 ; scanf("%d%d%d%d",&x1,&y1,&x2,&y2); if(x1==-1&&x2==-1&&y1==-1 &&y2==-1) break; } node[t].x= x1; node[t].y1 = y1; node[t].y2 = y2; node[t].flag = 1;//入边 y[t++] = y1; node[t].x = x2; node[t].y1 = y1; node[t].y2 = y2; node[t].flag = -1; y[t++] = y2; } return 0;}
0 0
- poj_1389 Area of Simple Polygons 线段树+扫描线
- [POJ1389]Area of Simple Polygons(扫描线+线段树)
- POJ 1389 Area of Simple Polygons 线段树 扫描线
- poj - 1389 - Area of Simple Polygons(离线化+扫描线+线段树)
- poj 1389 Area of Simple Polygons(线段树+扫描线)
- 【POJ 1389】Area of Simple Polygons(线段树+扫描线)
- POJ 1389 Area of Simple Polygons(线段树+扫描线+离散化)
- POJ1389Area of Simple Polygons【离散化+线段树+扫描线】
- hdoj 1542 && poj 1151 Atlantis && poj 1389 Area of Simple Polygons 【线段树 + 离散化 + 扫描线】
- POJ 1389 Area of Simple Polygons(线段树+扫描面积)
- POJ 1389 Area of Simple Polygons(扫描线)
- poj 1389 Area of Simple Polygons(线段树做法)
- poj1389Area of Simple Polygons(线段树+线扫描求矩形面积并)
- POJ1389-Area of Simple Polygons
- POJ1389[Area of Simple Polygons]
- poj1389 Area of Simple Polygons
- POJ 1389 Area of Simple Polygons 线段树求矩形面积
- (概念题)线段树 HOJ 1495 Area of Simple Polygons
- iOS常用控件、图片尺寸大集合
- xml字符串转为XmlDocument
- spring bean abstract="true"
- Protobuf开发环境搭建
- c#实现随鼠标移动窗体
- poj_1389 Area of Simple Polygons 线段树+扫描线
- java 调用存储过程 实例
- Xmanager Enterprise 4 图形工具
- XmlDocument 转换为xml字符串
- 初次项目,页面布局
- php 下载页面为doc文件
- android 生成pdf文件
- 声音的录制(VC)(保存为WAV文件) (转)
- 海南出差报告总结(案件录入与案件追踪系统)