HDU 1255 覆盖的面积 ( 扫描线 + 离散 求矩阵大于k次面积并 )
覆盖的面积
Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 4040 Accepted Submission(s): 1995
入数据的第一行是一个正整数T(1<=T<=100),代表测试数据的数量.每个测试数据的第一行是一个正整数N(1<=N&
lt;=1000),代表矩形的数量,然后是N行数据,每一行包含四个浮点数,代表平面上的一个矩形的左上角坐标和右下角坐标,矩形的上下边和X轴平行,
左右边和Y轴平行.坐标的范围从0到100000.
注意:本题的输入数据较多,推荐使用scanf读入数据.
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
#include <map>
using namespace std;
#define root 1,(n<<1)+10,1
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define lr rt<<1
#define rr rt<<1|1 const int N = ; struct Point {
double x , y1 , y2 ;
int v ;
bool operator < ( const Point &a ) const {
return x < a.x ;
}
}p[N];
int n , tot , tt ;
vector<double>e;
double c[N<<] ; void addpoint( double x , double y1 , double y2 , int v ) {
p[tot].x = x , p[tot].y1 = y1 , p[tot].y2 = y2 , p[tot].v = v , tot++ ;
}
map<double,int>mp; void lisan() {
mp.clear();
sort( e.begin() , e.end() ) ;
tt = ;
for( int i = ; i < e.size() ; ++i ) {
if( e[i] != e[i-] ) e[tt++] = e[i] ;
}
for( int i = ; i < tt ; ++i ) {
c[i+] = e[i] ;
mp[e[i]] = i+ ;
}
} int cnt[N<<] , lazy[N<<] ; double date[N<<] ; void build( int l , int r , int rt ) {
cnt[rt] = lazy[rt] = ;
if( l == r ) {
date[rt] = c[l+] - c[l];
return ;
}
int mid = (l+r)>>;
build(lson) , build(rson) ;
date[rt] = date[lr] + date[rr];
} void Down( int rt ) {
if( lazy[rt] ) {
cnt[lr] += lazy[rt] ;
cnt[rr] += lazy[rt] ;
lazy[lr] += lazy[rt] ;
lazy[rr] += lazy[rt] ;
lazy[rt] = ;
}
} void Up( int rt ) {
cnt[rt] = min( cnt[lr] , cnt[rr] );
} void update( int l , int r , int rt , int L , int R , int v ) {
if( l == L && r == R ) {
cnt[rt] += v ; lazy[rt] += v ;
return ;
}
int mid = (l+r)>> ;
Down(rt) ;
if( R <= mid ) update(lson,L,R,v) ;
else if( L > mid ) update( rson,L,R,v);
else update(lson,L,mid,v),update(rson,mid+,R,v);
Up(rt);
} double query( int l , int r , int rt , int L , int R ) {
if( cnt[rt] > ) return date[rt] ;
if( l == r ) return ;
Down(rt);
int mid = (l+r)>>;
if( R <= mid ) return query(lson,L,R);
else if( L > mid ) return query(rson,L,R);
else return query(lson,L,mid) + query(rson,mid+,R);
} int main() {
#ifdef LOCAL
freopen("in.txt","r",stdin);
#endif // LOCAL
int _ ; scanf("%d",&_);
while( _-- ) {
scanf("%d",&n);
e.clear();
tot = ;
for( int i = ; i < n ; ++i ) {
double x1 , y1 , x2 , y2 ;
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
if( y1 > y2 ) swap( y1 , y2 ) ;
addpoint( x1 , y1 , y2 , );
addpoint( x2 , y1 , y2 , - );
e.push_back(y1);
e.push_back(y2);
}
lisan();
sort( p , p + tot ) ;
build(,tt-,);
double ans = ;
update(,tt-,,mp[p[].y1],mp[p[].y2]-,p[].v);
for( int i = ; i < tot ; ++i ) {
ans += ( p[i].x - p[i-].x ) * query(,tt-,,,tt-) ;
update(,tt-,,mp[p[i].y1],mp[p[i].y2]-,p[i].v);
}
printf("%.2lf\n",ans);
}
return ;
}
HDU 1255 覆盖的面积 ( 扫描线 + 离散 求矩阵大于k次面积并 )的更多相关文章
- FZU 2187 回家种地 ( 扫描线 + 离散 求矩阵单次覆盖面积 )
2187 回家种地 Accept: 56 Submit: 230Time Limit: 1000 mSec Memory Limit : 32768 KB Problem Descript ...
- HDU 4417 Super Mario(划分树问题求不大于k的数有多少)
Super Mario Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total ...
- HDU - 1255 覆盖的面积(线段树求矩形面积交 扫描线+离散化)
链接:线段树求矩形面积并 扫描线+离散化 1.给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积. 2.看完线段树求矩形面积并 的方法后,再看这题,求的是矩形面积交,类同. 求面积时,用被覆 ...
- hdu 1255 覆盖的面积(求覆盖至少两次以上的面积)
了校赛,还有什么途径可以申请加入ACM校队? 覆盖的面积 Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/32768 K ...
- hdu 1255 覆盖的面积(线段树 面积 交) (待整理)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1255 Description 给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积. In ...
- HDU 1255 覆盖的面积 (线段树扫描线+面积交)
自己YY了一个的写法,不过时间复杂度太高了,网上的想法太6了 题意:给你一些矩阵,求出矩阵的面积并 首先按照x轴离散化线段到线段树上(因为是找连续区间,所以段建树更加好做). 然后我们可以想一下怎样 ...
- HDU 1255 覆盖的面积(线段树+扫描线)
题目地址:HDU 1255 这题跟面积并的方法非常像,仅仅只是须要再加一个变量. 刚開始我以为直接用那个变量即可,仅仅只是推断是否大于0改成推断是否大于1.可是后来发现了个问题,由于这个没有下放,没延 ...
- HDU 1255 覆盖的面积 (线段树+扫描线+离散化)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1255 题意很清楚,就是让你求矩阵之间叠加层数大于1的矩形块的面积和. 因为n只有1000,所以我离散化 ...
- hdu 1255 覆盖的面积 (线段树处理面积覆盖问题(模板))
http://acm.hdu.edu.cn/showproblem.php?pid=1255 覆盖的面积 Time Limit: 10000/5000 MS (Java/Others) Memo ...
随机推荐
- Redis的客户端Jedis
1. Redis支持消息的订阅与发布 Redis的消息订阅支持:先订阅后发布 订阅:subscribe c1 c2 发布:publish c2 hello-redis 支持通配符的订阅:psubscr ...
- html a标签链接点击闪动问题解决
<a href="#">链接点击会闪动,解决: 这三种都可以用:<a href="javascript:;"></a>< ...
- 输vim /etc/rc.d/init.d/mysqld 报错 …..localdomain.pid
在安装MySQL的过程中出现以下错误: 本人解决方法是: ps aux | grep mysqld 罗列出所有关于mysqld的进程,然后把关于mysqld的进程都kill了.
- 【改】DOS文件格式转UNIX文件格式
windows中的文本文件的换行符是"\r\n",而linux中是"\n",dos格式文件传输到unix系统时,会在每行的结尾多一个^M,当然也有可能看不到,但 ...
- Python PostgreSQL Psycopg2
[转] http://daigong.iteye.com/blog/901160 Python如果要操作Postgresql,需要一个API,这就需要Psycopg2 1. 链接PostgreSQL并 ...
- Linux 部署或升级openssh7.5p1
运维Linux系统,部署或升级openssh是经常面临的事,以下已redhat6和redhat7为例. 在redhat6中部署openssh会有什么坑,在编辑openssh源码包时会报一些类似的错误, ...
- 快速开发框架下载地址(github)
eladmin:https://github.com/elunez/eladmin bootDo:https://www.oschina.net/p/bootdo
- 在flask服务里面执行adb 安装命令
今天讲的是我在flask服务里执行adb install 命令,被坑了一天 需求是我们的安卓打包机器和我们的自动化机器不是同一台,管理人也不一样,想让每次打包后能通知到我们的服务器,然后执行后续的操作 ...
- zabbix邮件报警通过脚本来发送邮件
zabbix默认邮件报警会将各个报警接收人单独发送邮件,为了使邮件能以群发的方式统一一封邮件发送所有接收人,需要改成脚本的形式: sendemail.py: #!/usr/bin/python imp ...
- 一个能极大提高生产率的Chrome新建标签页扩展
我是一个对开发生产率有着BT需求的程序员,总是追求将自己的单位时间生产率最大化. 通过分析,我发现自己一天会反反复复使用Chrome的新建标签,然后访问常用的网站.因此,我期望新建一个默认的Chrom ...