欢迎访问~原文出处——博客园-zhouzhendong

去博客园看该题解


题目传送门 - POJ1151


题意概括

  给出n个矩形,求他们的面积并。

  n<=100


题解

  数据范围极小。

  我们分3种算法逐步优化。

  算法1: O(n3)

  如果这n个矩形的坐标都是整数,而且比较小,那么我们显然可以用最暴力的方法:一个一个打标记。

  但是不是这样的。

  坐标大小很大,而且是实数。

  然而我们发现差不多,只要先离散化一下,然后再打标记即可。

  算法2:O(n2)

  实际上,上面的方法十分慢。如果n的范围到了1000,上面的就无济于事了。

  而实际上,基于上面的打标记的算法,我们可以通过差分的方法n2解决。

  我们通过差分,可以用n2的时间标记,n2的时间判断每一个区域是否被覆盖。

  空间复杂度O(n2)

  算法3:O(n logn) 扫描线

  实际上,这类问题的数据范围可以到100000这个级别。

  矩形面积并可以用扫描线算法来解决。先看原理,后面讲具体实现。

  比如下图:

  

  当前我们的扫描线到达了淡黄色部分。

  由于之前没有记录,所以答案不增加。

  然而我们记下当前横向覆盖的长度。

  然后我们到了第二条扫描线,加上原来记录的横向覆盖长度乘以增加的高度就是当前增加的答案。

  然后,我们更新了横向覆盖的长度。

  继续。

  

  然后第三条。现在的横向覆盖长度是两边加起来,所以增加的面积是两块了。

  然后更新横向覆盖的长度,加上了中间的那一条。

  然后继续。

  

  现在有这么长的一条都是被横向覆盖的了。

  所以新增的面积是浅蓝色部分。

  然后我们发现左上那条线是出边,所以要删除这一条线。

  所以横向覆盖的长度为如下:

  

  同理,接下来是:

  

  然后就OK了。

  

  那么具体怎么实现呢?

  我们开一棵线段树来维护!

  在读入之后,我们把所有的横线都拆开,分成下边和上边两类。某一区间在进入下边的时候+1,离开上边的时候-1,所以我们分别给上下边标记+1和-1。

  对于Y,我们离散化一下。

  对于X,我们按照边的X排一个序。

  然后按照刚才那样的处理。

  具体如何维护详见代码。


代码

#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
using namespace std;
const int N=+,M=N*;
const double Eps=1e-;
int T=,n,m,tot_Y,tot_s;
double Y[M];
struct Segment{
double x,L,R;
int v;
void set(double x_,double L_,double R_,int v_){
x=x_,L=L_,R=R_,v=v_;
}
}s[M];
struct SegTree{
int cnt;
double sum;
}t[M*];
bool cmp_s(Segment a,Segment b){
return a.x<b.x;
}
void build(int rt,int le,int ri){
t[rt].cnt=;
t[rt].sum=;
if (le==ri)
return;
int mid=(le+ri)>>,ls=rt<<,rs=ls|;
build(ls,le,mid);
build(rs,mid+,ri);
}
void pushup(int rt,int le,int ri){
int ls=rt<<,rs=ls|;
if (t[rt].cnt)
t[rt].sum=Y[ri+]-Y[le];
else if (le==ri)
t[rt].sum=;
else
t[rt].sum=t[ls].sum+t[rs].sum;
}
void update(int rt,int le,int ri,int xle,int xri,int d){
if (le>xri||ri<xle)
return;
if (xle<=le&&ri<=xri){
t[rt].cnt+=d;
pushup(rt,le,ri);
return;
}
int mid=(le+ri)>>,ls=rt<<,rs=ls|;
update(ls,le,mid,xle,xri,d);
update(rs,mid+,ri,xle,xri,d);
pushup(rt,le,ri);
}
int find_double(double x){
int le=,ri=m,mid;
while (le<=ri){
mid=(le+ri)>>;
if (abs(x-Y[mid])<Eps)
return mid;
if (Y[mid]<x)
le=mid+;
else
ri=mid-;
}
}
int main(){
while (scanf("%d",&n)&&n){
tot_Y=tot_s=;
for (int i=;i<=n;i++){
double xA,yA,xB,yB;
scanf("%lf%lf%lf%lf",&xA,&yA,&xB,&yB);
if (yB-yA<Eps||xB-xA<Eps)
continue;
Y[++tot_Y]=yA,Y[++tot_Y]=yB;
s[++tot_s].set(xA,yA,yB,);
s[++tot_s].set(xB,yA,yB,-);
}
sort(Y+,Y+tot_Y+);
sort(s+,s+tot_s+,cmp_s);
m=;
for (int i=;i<=tot_Y;i++)
if (Y[i]-Y[i-]>Eps)
Y[++m]=Y[i];
build(,,m);
double ans=;
for (int i=;i<=tot_s;i++){
ans=ans+(s[i].x-s[i-].x)*t[].sum;
int L=find_double(s[i].L);
int R=find_double(s[i].R);
update(,,m,L,R-,s[i].v);
}
printf("Test case #%d\nTotal explored area: %.2lf\n\n",++T,ans);
}
return ;
}

POJ1151Atlantis 矩形面积并 扫描线 线段树的更多相关文章

  1. 矩形面积并-扫描线 线段树 离散化 模板-poj1151 hdu1542

    今天刚看到这个模板我是懵逼的,这个线段树既没有建树,也没有查询,只有一个update,而且区间成段更新也没有lazy标记....研究了一下午,我突然我发现我以前根本不懂扫描线,之所以没有lazy标记, ...

  2. POJ 1151 Atlantis 矩形面积求交/线段树扫描线

    Atlantis 题目连接 http://poj.org/problem?id=1151 Description here are several ancient Greek texts that c ...

  3. hdu1542 矩形面积并(线段树+离散化+扫描线)

    题意: 给你n个矩形,输入每个矩形的左上角坐标和右下角坐标. 然后求矩形的总面积.(矩形可能相交). 题解: 前言: 先说说做这道题的感受: 刚看到这道题顿时就懵逼了,几何 烂的渣渣.后来从网上搜题解 ...

  4. 【HDU 1542】Atlantis 矩形面积并(线段树,扫描法)

    [题目] Atlantis Problem Description There are several ancient Greek texts that contain descriptions of ...

  5. luogu P1856 [USACO5.5]矩形周长Picture 扫描线 + 线段树

    Code: #include<bits/stdc++.h> #define maxn 200007 #define inf 100005 using namespace std; void ...

  6. HDU 1255 覆盖的面积 (扫描线 线段树 离散化 矩形面积并)

    题目链接 题意:中文题意. 分析:纯手敲,与上一道题目很相似,但是刚开始我以为只是把cnt>=0改成cnt>=2就行了,. 但是后来发现当当前加入的线段的范围之前 还有线段的时候就不行了, ...

  7. 【BZOJ4418】[Shoi2013]扇形面积并 扫描线+线段树

    [BZOJ4418][Shoi2013]扇形面积并 Description 给定N个同心的扇形,求有多少面积,被至少K个扇形所覆盖. Input 第一行是三个整数n,m,k.n代表同心扇形的个数,m用 ...

  8. (HDU 1542) Atlantis 矩形面积并——扫描线

    n个矩形,可以重叠,求面积并. n<=100: 暴力模拟扫描线.模拟赛大水题.(n^2) 甚至网上一种“分块”:分成n^2块,每一块看是否属于一个矩形. 甚至这个题就可以这么做. n<=1 ...

  9. HDU - 1255 覆盖的面积(线段树求矩形面积交 扫描线+离散化)

    链接:线段树求矩形面积并 扫描线+离散化 1.给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积. 2.看完线段树求矩形面积并 的方法后,再看这题,求的是矩形面积交,类同. 求面积时,用被覆 ...

随机推荐

  1. Linux 4.10.8 根文件系统制作(三)---制作yaffs文件系统

    这里直接用的是韦东山提供的工具. yaffs文件系统是专门用于nandflash的文件系统 我们前面已经下载了yaffs 的源码,也做了文件系统目录的构建. 在yaffs2的源码目录中有一个utils ...

  2. Linux - sort 排序

    -t # 指定排序时所用的栏位分隔字符 -n # 依照数值的大小排序 -r # 以相反的顺序来排序 -f # 排序时,将小写字母视为大写字母 -d # 排序时,处理英文字母.数字及空格字符外,忽略其他 ...

  3. Java 学习札记(三)免安装版TomCat中tomcat6w.exe的运行

    1.使用环境 很多时候我们用的是官网的解压免安装版的Tomcat,相比安装Tomcat除了少了安装步骤以外还少了tomcat6w.exe运行所需要的环境变量,所以一般Java开发免安装版的已经足够使用 ...

  4. python3中的比较函数

    在py2中,比较函数是cmp,而在py3,cmp已经不存在了,Py3启用了新的比较方法 原来在py2中,a>b就会调用a对象中的__cmp__函数,而现在a>b会调用a对象中的__lt__ ...

  5. I - Older Brother Gym - 101490I

    题目链接:https://cn.vjudge.net/problem/Gym-101490I 题目大意:给你一个整数,问你这个整数能不能表示成一个素数的k次方? 具体思路:对于每一个数,我们先判断他是 ...

  6. jquery 学习(四) - 标签 添加/删除/修改

    HTML代码 <div class="a1"> <div> <span id="a2">aaa</span> & ...

  7. 关于形如--error LNK2005: xxx 已经在 msvcrtd.lib ( MSVCR90D.dll ) 中定义--的问题分析解决

    转自:http://hi.baidu.com/qinfengxiaoyue/item/ff262ccfb53b4c2ba0b50a89 引自:http://blog.csdn.net/sptoor/a ...

  8. 内核IS_ERR宏解析 【转】

    转自:http://blog.chinaunix.net/uid-20196318-id-28769.html 最近在使用filp_open打开文件时遇到到一个问题,当打开一个并不存在的文件时,fil ...

  9. 在SecureCRT中做make menuconfig乱码

      不能在SecureCRT中做(显示为乱码),从高手那里学来一招,解决了这个问题: options--terminal--emulation-- xterm  ansi color1.先设置终端为x ...

  10. windows下解压zip包,包含中文解析

    #coding=utf8 import os import zipfile import sys,locale # 本来以为需要,结果不需要 # def p(f): # #print '%s.%s() ...