POJ 1151 Atlantis 求矩阵面积并 扫描线 具体解释
题意:
给定n个矩阵的左下角和右上角坐标,求矩阵面积并(矩阵总是正放的,即与x轴y轴都平行)
思路:
扫描线裸题
http://www.cnblogs.com/fenshen371/p/3214092.html
对于一个矩阵,我们仅仅关心他的上下底边。线段树维护的是当前有哪些区间是被线段覆盖的(即把线段都移动到x轴上,x轴被覆盖的长度是多少)
我们从上往下扫,则对于随意一条(如15,10) 我们会计算(20,13)-(15,10)之间的矩形面积
而对于(11,8),我们会计算(15,10)-(11,8)之间的矩形面积。
矩形面积= 底边长*高
1、显然这里的高就是(15,10)-(11,8) ; 更通俗的就是:line[i-1].y - line[i].y;
2、而底边长则是取决于从(11,8) 向上(不包含(11,8))有多长的区间被线段覆盖了,这个长度则是用线段树来维护。
3、而(11,8)这条线段对于线段树或者说 对区间是覆盖还是删除则取决于这条线段是上底边还是下底边(上底边则是覆盖,下底边则是删除)
对于线段树的建树也要注意一个要点:
平时建树的习惯是:
[1,4]
/ \
[1,2] [3,4]
而我们用线段树维护的是x轴坐标,是连续的区间,假设这样建树则忽略掉了区间[2,3]
所以我们建成如此:
[1,4]
/ \
[1,2] [2,4]
即: build(l,mid,L(id)); build(mid,r,R(id)); mid不+1了
推断某区间是否为最底层区间时则写成 if(tree[id].l +1 == tree[id].r) ;
由于[1,1]这个区间对我们是没有意义的,即 l < r才是有长度的区间,才是有意义的。而最小的l < r 就是 l +1 = r 的情况。
代码:
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
#define N 205
#define L(x) (x<<1)
#define R(x) (x<<1|1)
int Mid(int a,int b){return (a+b)>>1;}
//这份代码是对y轴建树,从左到右扫
struct node{ //这个区间被c条线段覆盖了
int l, r;
int c; //用来记录重叠情况
double cnt, lf, rf;//离散点l相应的实际点lf cnt表示这个区间内存在的线段长度
}tree[N*4];
struct Line{
double x,y1,y2;
int f;
Line(double a=0.0,double b=0.0,double c=0.0, int d=0):x(a),y1(b),y2(c),f(d){}
}line[N*2];
bool cmp(Line a, Line b){return a.x<b.x;}
double y[N];
void build(int l, int r, int id){
tree[id].l = l, tree[id].r = r;
tree[id].cnt = tree[id].c = 0;
tree[id].lf = y[l]; tree[id].rf = y[r];
if(l+1==r)return ;
int mid = Mid(l,r);
build(l, mid, L(id));
build(mid,r,R(id));
}
void push_up(int id){
if(tree[id].c){
tree[id].cnt = tree[id].rf - tree[id].lf;
return ;
}
if(tree[id].l+1 == tree[id].r) tree[id].cnt = 0;
else tree[id].cnt = tree[L(id)].cnt + tree[R(id)].cnt;
}
void update(int id, Line e){
if(e.y1 == tree[id].lf && tree[id].rf == e.y2){
tree[id].c += e.f;
push_up(id); return ;
}
if(e.y2 <= tree[L(id)].rf)
update(L(id),e);
else if(tree[R(id)].lf <= e.y1)
update(R(id),e);
else {
Line tmp = e; tmp.y2 = tree[L(id)].rf; update(L(id),tmp);
tmp = e; tmp.y1 = tree[R(id)].lf; update(R(id),tmp);
}
push_up(id);
}
set<double>myset;
set<double>::iterator p;
int main(){
int i, j, n, Cas = 1;
double x1,x2,y1,y2;
while(scanf("%d",&n),n){
myset.clear();
int top = 1;
for(i = 1;i <= n; i++){
scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2);
myset.insert(y1); myset.insert(y2);
Line E = Line(x1,y1,y2,1);
line[top++] = E;
Line E2= Line(x2,y1,y2,-1);
line[top++] = E2;
}
sort(line+1, line+top, cmp);
int siz = 1;
for(p=myset.begin(); p!=myset.end(); p++) y[siz++] = *p;
build(1, siz-1, 1);
update(1, line[1]);
double ans = 0;
for(i = 2; i < top; i++){
ans += tree[1].cnt * (line[i].x - line[i-1].x);
update(1, line[i]);
}
printf("Test case #%d\nTotal explored area: %.2lf\n\n",Cas++,ans);
}
return 0;
}
POJ 1151 Atlantis 求矩阵面积并 扫描线 具体解释的更多相关文章
- POJ 1151 Atlantis(扫描线)
题目原链接:http://poj.org/problem?id=1151 题目中文翻译: POJ 1151 Atlantis Time Limit: 1000MS Memory Limit: 10 ...
- HDU - 1255 覆盖的面积(线段树求矩形面积交 扫描线+离散化)
链接:线段树求矩形面积并 扫描线+离散化 1.给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积. 2.看完线段树求矩形面积并 的方法后,再看这题,求的是矩形面积交,类同. 求面积时,用被覆 ...
- POJ 1151 Atlantis(经典的线段树扫描线,求矩阵面积并)
求矩阵的面积并 采用的是区间更新 #include <iostream> #include <stdio.h> #include <string.h> #inclu ...
- POJ 1151 Atlantis 线段树+离散化+扫描线
这次是求矩形面积并 /* Problem: 1151 User: 96655 Memory: 716K Time: 0MS Language: G++ Result: Accepted */ #inc ...
- POJ 1151 Atlantis(离散化)
点我看题目 题意 : 就是给你n个矩形的最左下角和最右上角的点的坐标,然后将这n个矩形的面积求出来. 思路 : 离散化求矩形面积并.离散化具体介绍.将横纵坐标离散开来分别存,然后排序,也可以按照黑书上 ...
- Area - POJ 1654(求多边形面积)
题目大意:从原点开始,1-4分别代表,向右下走,向右走,向右上走,向下走,5代表回到原点,6-9代表,向上走,向左下走,向左走,向左上走.求出最后的多边形面积. 分析:这个多边形面积很明显是不规则的, ...
- POJ 1151 Atlantis 矩形面积求交/线段树扫描线
Atlantis 题目连接 http://poj.org/problem?id=1151 Description here are several ancient Greek texts that c ...
- hdu 1542&&poj 1151 Atlantis[线段树+扫描线求矩形面积的并]
Atlantis Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Total S ...
- POJ 1151 Atlantis(线段树-扫描线,矩形面积并)
题目链接:http://poj.org/problem?id=1151 题目大意:坐标轴上给你n个矩形, 问这n个矩形覆盖的面积 题目思路:矩形面积并. 代码如下: #include<stdio ...
随机推荐
- H5新特性:video与audio的使用
HTML5 DOM 为 <audio> 和 <video> 元素提供了方法.属性和事件. 这些方法.属性和事件允许您使用 JavaScript 来操作 <audio> ...
- Ubuntu 登陆循环启动 无法进入桌面 libGL error: failed to load driver: swrast
导致无法进入Ubuntu图形界面的原因很多,关键是要找到问题原因所在. 最佳的方法是查看.xsession-errors这个日志文件,依据问题解决~ 这里先汇总一下问题: 1 环境变量导致的,进入tt ...
- OBJECT_ID()的使用方法
数据库中每个对像都有一个唯一的ID值,用Object_name(id)可以根据ID值得到对像的名称,object_id(name)可以根据对像名称得到对象的ID object_id()只能返回用户创建 ...
- Python学习笔记:个税起征点上调至5000,算一算少交多少税?
一.旧税率表与新税率表比较 以前起征点是3500,2018年10月1日起起征点正式修改为5000,下面我们用Python来分别计算新旧个人所得税分别为多少? 二.旧的个人所得税 import sys ...
- 编程六月定律 | 外刊IT评论网
编程六月定律 上周,我被迫对一个很老的项目做一些修改.麻烦是,当开始着手时,我真的记不清这个项目究竟有多老了. 这实际上是我使用Codeigniter实现的第一个MVC项目.打开项目文件后,很多东西都 ...
- 课堂实验-String类和Arrays类
课堂实验 在IDEA中以TDD的方式对String类和Arrays类进行学习 测试相关方法的正常,错误和边界情况 String类 charAt split Arrays类 sort binarySea ...
- 判断一个字符是否为数字的两种方法(C/C++)
在平时,我们经常遇见判断字符是否为数字这种题目,虽然感觉还是很简单,不过我是个更喜欢用函数的人,因为我觉得这样更便捷,所以我更推荐第二种方式. 1.直接判断 #include <stdio.h& ...
- JS黑魔法之this, setTimeout/setInterval, arguments
最近发现了JavaScript Garden这个JS黑魔法收集处,不过里面有一些东西并没有说得很透彻,于是边看边查文档or做实验,写了一些笔记,顺手放在博客.等看完了You don't know JS ...
- jquery全面判断是否IE6浏览器
今天在写一个登录回车提交表单的操作时,出现keydown在IE6下不能提交的兼容问题,随之无奈,找到可以使用keyup或者keypress事件时间来兼容所有浏览器,但是呢体验效果很不友好,所以只能委屈 ...
- HTTP Status 500 - Request processing failed; nested exception is org.apache.ibatis.binding.BindingException
在使用Maven工程管理工具整合SSM框架时,Mybatis使用逆向工程生成的pojo,mapper接口及映射文件,把mapper接口和映射文件放在DAO工程的同一级src/main/java目录下. ...