【 HDU 1255】 覆盖的面积(矩阵面积交,线段树,扫描法)
【题目】
覆盖的面积
Problem Description给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积.
Input输入数据的第一行是一个正整数T(1<=T<=100),代表测试数据的数量.每个测试数据的第一行是一个正整数N(1<=N<=1000),代表矩形的数量,然后是N行数据,每一行包含四个浮点数,代表平面上的一个矩形的左上角坐标和右下角坐标,矩形的上下边和X轴平行,左右边和Y轴平行.坐标的范围从0到100000.注意:本题的输入数据较多,推荐使用scanf读入数据.
Output对于每组测试数据,请计算出被这些矩形覆盖过至少两次的区域的面积.结果保留两位小数.Sample Input2
5
1 1 4 2
1 3 3 7
2 1.5 5 4.5
3.5 1.25 7.5 4
6 3 10 7
3
0 0 1 1
1 0 2 1
2 0 3 1Sample Output7.63
0.00AuthorIgnatius.L & weigang Lee
- #include<cstdio>
- #include<cstdlib>
- #include<cstring>
- #include<iostream>
- #include<algorithm>
- #include<cmath>
- using namespace std;
- #define Maxn 100010
- struct hp
- {
- int x1,x2,flag;
- double y;
- }t[Maxn*];int tl;
- struct lsh
- {
- double x;
- int id;
- }q[Maxn*];int ql;
- bool cmp(lsh x,lsh y) {return x.x<y.x;}
- bool cmp2(hp x,hp y) {return x.y<y.y;}
- int mymin(int x,int y) {return x<y?x:y;}
- double dis[*Maxn];
- struct node
- {
- int l,r,lc,rc,cnt;
- double sm,fg;
- }tr[Maxn*];int len;
- int build(int l,int r)
- {
- int x=++len;
- tr[x].l=l;tr[x].r=r;
- tr[x].sm=tr[x].cnt=;tr[x].fg=;
- if(l<r-)
- {
- int mid=(tr[x].l+tr[x].r)>>;
- tr[x].lc=build(l,mid);
- tr[x].rc=build(mid,r);
- }
- else tr[x].lc=tr[x].rc=;
- return x;
- }
- void upd(int x)
- {
- int l=tr[x].l,r=tr[x].r,lc=tr[x].lc,rc=tr[x].rc;
- if(tr[x].cnt!=) tr[x].fg=dis[r]-dis[l];
- else tr[x].fg=tr[lc].fg+tr[rc].fg;
- if(tr[x].cnt>=) tr[x].sm=dis[r]-dis[l];
- else if(tr[x].cnt==) tr[x].sm=tr[lc].fg+tr[rc].fg;
- else tr[x].sm=tr[lc].sm+tr[rc].sm;
- }
- void change(int x,int l,int r,int c)
- {
- if(tr[x].l==l&&tr[x].r==r)
- {
- tr[x].cnt+=c;
- upd(x);
- return;
- }
- int mid=(tr[x].l+tr[x].r)>>;
- if(r<=mid) change(tr[x].lc,l,r,c);
- else if(l>=mid) change(tr[x].rc,l,r,c);
- else
- {
- change(tr[x].lc,l,mid,c);
- change(tr[x].rc,mid,r,c);
- }
- upd(x);
- }
- int main()
- {
- int T;
- int n;
- scanf("%d",&T);
- while(T--)
- {
- scanf("%d",&n);
- if(n==) break;
- tl=;ql=;
- for(int i=;i<=n;i++)
- {
- double x1,y1,x2,y2;
- scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
- //t[++tl].x1,
- t[++tl].flag=;t[tl].y=y1;
- t[++tl].flag=-;t[tl].y=y2;
- q[++ql].x=x1,q[ql].id=tl;
- q[++ql].x=x2,q[ql].id=tl+*n;
- }
- sort(q+,q++ql,cmp);
- int p=;
- for(int i=;i<=ql;i++)
- {
- if(q[i].x!=q[i-].x||p==) p++,dis[p]=dis[p-]+q[i].x-q[i-].x;
- if(q[i].id<=*n) t[q[i].id].x1=t[q[i].id-].x1=p;
- else t[q[i].id-*n].x2=t[q[i].id--*n].x2=p;
- }
- sort(t+,t++tl,cmp2);
- len=;tr[].sm=;
- build(,p);
- double ans=;
- for(int i=;i<tl;i++)
- {
- // printf("%d %d\n",t[i].x1,t[i].x2);
- change(,t[i].x1,t[i].x2,t[i].flag);
- ans+=tr[].sm*(t[i+].y-t[i].y);
- }
- printf("%.2lf\n",ans);
- // printf("Test case #%d\nTotal explored area: %.2lf\n\n",++kase,ans);
- }
- return ;
- }
2016-11-10 14:37:49
【 HDU 1255】 覆盖的面积(矩阵面积交,线段树,扫描法)的更多相关文章
- hdu1255 矩阵的交 线段树+扫描线
/* 不是叶子节点 ,且cnt=1.注意这里,cnt=1确切的意义是什么, 应该是,可以确定,这个区间被完全覆盖了1次, 而有没有被完全覆盖两次或以上则不知道无法确定,那么怎么怎么办了, 只要加上t[ ...
- hdu 1255 覆盖的面积(线段树 面积 交) (待整理)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1255 Description 给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积. In ...
- HDU - 1255 覆盖的面积(线段树求矩形面积交 扫描线+离散化)
链接:线段树求矩形面积并 扫描线+离散化 1.给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积. 2.看完线段树求矩形面积并 的方法后,再看这题,求的是矩形面积交,类同. 求面积时,用被覆 ...
- hdu 1255 覆盖的面积(求覆盖至少两次以上的面积)
了校赛,还有什么途径可以申请加入ACM校队? 覆盖的面积 Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/32768 K ...
- hdu 5700区间交(线段树)
区间交 Time Limit: 8000/4000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submiss ...
- hdu 5266 pog loves szh III(lca + 线段树)
I - pog loves szh III Time Limit:6000MS Memory Limit:131072KB 64bit IO Format:%I64d & %I ...
- HDU 2795 Billboard(宣传栏贴公告,线段树应用)
HDU 2795 Billboard(宣传栏贴公告,线段树应用) ACM 题目地址:HDU 2795 Billboard 题意: 要在h*w宣传栏上贴公告,每条公告的高度都是为1的,并且每条公告都要 ...
- HDU 1255 覆盖的面积 (线段树扫描线+面积交)
自己YY了一个的写法,不过时间复杂度太高了,网上的想法太6了 题意:给你一些矩阵,求出矩阵的面积并 首先按照x轴离散化线段到线段树上(因为是找连续区间,所以段建树更加好做). 然后我们可以想一下怎样 ...
- hdu 1255 覆盖的面积 (扫描线求矩形交)
覆盖的面积 Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Subm ...
- hdu 1255 覆盖的面积 (线段树处理面积覆盖问题(模板))
http://acm.hdu.edu.cn/showproblem.php?pid=1255 覆盖的面积 Time Limit: 10000/5000 MS (Java/Others) Memo ...
随机推荐
- Python中如何写控制台进度条的整理
本文实例讲述了Python显示进度条的方法,是Python程序设计中非常实用的技巧.分享给大家供大家参考.具体方法如下: 首先,进度条和一般的print区别在哪里呢? 答案就是print会输出一个\n ...
- FWFT FIFO读操作注意
FWFT:First Word Fall Through的缩写,好像是Xilinx的说法,Altera对应的概念是Show-ahead synchronous(SASO).即数据在rdreq有效之前就 ...
- .Net Framework 各个版本新特性总结 (一)
.Net Framework 4.5 新特性 最近面试时又看到有问.Net Framework 新特性的问题,一时被问到了.平时也是拿起来就用,新版本出来了,新特性也就是瞄一眼,也没去仔细查看.这次干 ...
- 介绍SmartUpload很好的网站
附带链接:http://www.cnblogs.com/elleniou/archive/2012/09/24/2700583.html
- 首页的sitecontent地址
当无法出现工具栏时,利用url地址,可以方便直接进入sitecontent http://sp2013/sites/bentest/_layouts/15/viewlsts.aspx
- mysql数据库文件默认保存目录(windows)
如果没有自己去设置安装路径,MYSQL默认安装在C:/Program Files/MySQL/MySQL Server 5.1,新建的数据库文件在C:/Documents and Settings/A ...
- oracle数据库TNS
TNS是Oracle Net的一部分,专门用来管理和配置Oracle数据库和客户端连接的一个工具,在大多数情况下客户端和数据库要通讯,必须配置TNS,下面看一如何配置它吧: TNS简要介绍与应用 :O ...
- mysql学习笔记(1)
参考教材<MySQL入门经典> 王雨竹 高飞 机械工业出版社 软件下载:http://www.mysql.com 安装好后打开命令提示符 (黑窗口) net start mys ...
- [翻译][MVC 5 + EF 6] 3:排序、过滤、分页
原文:Sorting, Filtering, and Paging with the Entity Framework in an ASP.NET MVC Application 1.添加排序: 1. ...
- 一次plsql 问题记录
环境 : window 7 x64 oracle 10.2g plsql 10.0.5 问题是 新装的 oracle10.2 plsql 一直连接不上 ,oracle_home 配置都对 .sql ...