题目链接:点击打开链接

题目描写叙述:给定一些矩形,求这些矩形的总面积。假设有重叠。仅仅算一次

解题思路:扫描线+线段树+离散(代码从上往下扫描)

代码:

#include<cstdio>
#include <algorithm>
#define MAXN 110
#define LL ((rt<<1)+1)
#define RR ((rt<<1)+2)
using namespace std;
int n;
struct segment{
double l,r,h;
int f;
bool operator<(const segment& b)const{
return h>b.h;
}
}sg[2*MAXN];
double pos[2*MAXN];
int id;
void addSegment(double x1,double y1,double x2,double y2){
sg[id].l=x1;sg[id].r=x2;
sg[id].h=y1;sg[id].f=1;
pos[id++]=x1;
sg[id].l=x1;sg[id].r=x2;
sg[id].h=y2;sg[id].f=-1;
pos[id++]=x2;
}
int binary(double key,int low,int high){
while(low<=high){
int mid=(low+high)/2;
if(pos[mid]==key)
return mid;
else if(key<pos[mid])
high=mid-1;
else
low=mid+1;
}
return -1;
}
struct Tree{
int l,r;
int cover;
double len;
}tree[8*MAXN];
void build(int rt,int l,int r){
tree[rt].l=l;
tree[rt].r=r;
tree[rt].cover=0;
tree[rt].len=0;
if(l==r-1)
return;
int mid=(l+r)>>1;
build(LL,l,mid);
build(RR,mid,r);
}
void pushup(int rt){
if(tree[rt].cover)
tree[rt].len=pos[tree[rt].r]-pos[tree[rt].l];
else if(tree[rt].l==tree[rt].r-1)
tree[rt].len=0;
else
tree[rt].len=tree[LL].len+tree[RR].len;
}
void update(int rt,int l,int r,int f){
if(tree[rt].l==l&&tree[rt].r==r){
tree[rt].cover+=f;
pushup(rt);
return;
}
int mid=(tree[rt].l+tree[rt].r)>>1;
if(r<=mid)
update(LL,l,r,f);
else if(l>=mid)
update(RR,l,r,f);
else{
update(LL,l,mid,f);
update(RR,mid,r,f);
}
pushup(rt);
}
int main(){
int Case=0;
while(scanf("%d",&n)!=EOF&&n!=0){
id=0;
double x1,y1,x2,y2;
for(int i=0;i<n;++i){
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
addSegment(x1,y1,x2,y2);
}
n=(n<<1);
sort(sg,sg+n);
sort(pos,pos+n);
int m=1;
for(int i=1;i<n;++i)
if(pos[i]!=pos[i-1])
pos[m++]=pos[i];
build(0,0,m-1);
double ans=0;
int l=binary(sg[0].l,0,m-1);
int r=binary(sg[0].r,0,m-1);
update(0,l,r,sg[0].f);
for(int i=1;i<n;i++){
ans+=(sg[i-1].h-sg[i].h)*tree[0].len;
l=binary(sg[i].l,0,m-1);
r=binary(sg[i].r,0,m-1);
update(0,l,r,sg[i].f);
}
printf("Test case #%d\n",++Case);
printf("Total explored area: %.2f\n\n",ans);
}
return 0;
}

hdu1542 Atlantis(扫描线+线段树+离散)矩形相交面积的更多相关文章

  1. [HDU1542]Atlantis(扫描线+线段树)

    Atlantis Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Su ...

  2. hdu1542 Atlantis (线段树+扫描线+离散化)

    Atlantis Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total S ...

  3. [POJ1151][HDU1542]Atlantis(线段树,扫描线)

    英文题面,我就只放个传送门了. Solution  题意是算矩形面积并,这是扫描线算法能解决的经典问题. 算法的大致思想是,把每一个矩形拆成上边和下边(以下称作扫描线),每条扫描线有四个参数l,r,h ...

  4. hdu-1542 Atlantis(离散化+线段树+扫描线算法)

    题目链接: Atlantis Time Limit: 2000/1000 MS (Java/Others)     Memory Limit: 65536/32768 K (Java/Others) ...

  5. poj1151 Atlantis——扫描线+线段树

    题目:http://poj.org/problem?id=1151 经典的扫描线问题: 可以用线段树的每个点代表横向被矩形上下边分割开的每一格,这样将一个矩形的出现或消失化为线段树上的单点修改: 每个 ...

  6. POJ 1151 Atlantis (扫描线+线段树)

    题目链接:http://poj.org/problem?id=1151 题意是平面上给你n个矩形,让你求矩形的面积并. 首先学一下什么是扫描线:http://www.cnblogs.com/scau2 ...

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

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

  8. POJ 1151Atlantis 扫描线+线段树求矩形面积并

    题目链接 #include <iostream> #include <vector> #include <cstdio> #include <cstring& ...

  9. 【POJ1151】Atlantis(线段树,扫描线)

    [POJ1151]Atlantis(线段树,扫描线) 题面 Vjudge 题解 学一学扫描线 其实很简单啦 这道题目要求的就是若干矩形的面积和 把扫描线平行于某个轴扫过去(我选的平行\(y\)轴扫) ...

随机推荐

  1. codeforces——思路与规律

    codeforces 804B     http://codeforces.com/problemset/problem/804/B /* 题意:给定一个只含ab的序列,每次操作可将ab变为bba 问 ...

  2. HangFire的定时任务和Quartz.NET总结(三)Quartz 配置

    在一个Net Core需求中,需要在每天的凌晨三点去抓取两个电商仓库的剩余的每个料号的数量来写会自己的表中, 用到了HangFire的定时任务 这篇文章讲的很详细记录下   文章2  这篇更简单 Qu ...

  3. servlet-后台获取form表单传的参数

    前台代码: <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> <html> & ...

  4. Table is specified twice, both as a target for 'UPDATE' and as a separate source

    UPDATE Bins b SET b.ShopSn =’111201611111168706’ WHERE b.Id IN (SELECT b.Id FROM Bins b JOIN BinInve ...

  5. c# ado.net eftity framework 返回多表查询结果

    public static IQueryable GetWeiXinTuWenList() { using (var Model = new Model.WeiXinEntities()) { var ...

  6. VHDL之concurrent之when

    WHEN (simple and selected) It is one of the fundamental concurrent statements (along with operators ...

  7. C#的WebBrowser操作frame

    刚学c#不久,也不太懂什么IHTMLDocument.IHTMLDocument2.IWebBrowser2等等.自己琢磨了好久,终于知道了怎么用WebBrowser操作frame和iframe. 1 ...

  8. 【JSP】上传图片到数据库中

    第一步:建立数据库 create table test_img(id number(4),name varchar(20),img long raw); 第二步:(NewImg.html) <h ...

  9. extension Array where Element 代码学习

    var fieldNames: [String] { let p = UnsafePointer<Int32>(self.pointer) return Array(utf8Strings ...

  10. BZOJ 4817: [Sdoi2017]树点涂色 LCT+Access的性质+DFS序+线段树

    Code: #include<bits/stdc++.h> #define maxn 200003 #define inf -1000000 using namespace std; vo ...