链接:线段树求矩形面积并 扫描线+离散化

1、给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积.

2、看完线段树求矩形面积并 的方法后,再看这题,求的是矩形面积交,类同。

求面积时,用被覆盖2次以上的那一段乘以扫描线的距离即可,具体实现见代码。

3、

/*
HDU 1255 覆盖的面积
求矩形面积交(离散化+线段树)
给定一些矩形
求被这些矩形覆盖过至少两次的区域的面积这里的方法是:线段树求矩形面积交 扫描线+离散化
左右扫描(x轴扫描),把y轴上的线段离散化
*/
#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;const int MAXN=2005;struct Node
{int l,r;//离散化后的y值double lf,rf;//原来的y值,lf代表l原来的y值,rf代表r原来的y值int cnt;//当前段被覆盖过几次double lenOnce;//当前段中覆盖一次以上的长度double lenTwice;//当前段中覆盖两次以上的长度
} segTree[MAXN*4];
struct Line
{double x;//也就是对应的扫描线x坐标double y1,y2;//线段的端点坐标double f;//1表示一个矩形左边的边,-1表示右边的边
} line[MAXN];double y[MAXN];bool cmp(Line a,Line b)
{return a.x<b.x;
}void Build(int i,int l,int r)
{segTree[i].l=l;segTree[i].r=r;segTree[i].lenOnce=0;segTree[i].lenTwice=0;segTree[i].lf=y[l];segTree[i].rf=y[r];if(l+1==r)return;int mid=(l+r)>>1;Build(i<<1,l,mid);Build((i<<1)|1,mid,r);
}
void calen(int i)
{if(segTree[i].cnt>=2)//被线段覆盖次数>=2
    {segTree[i].lenOnce=segTree[i].rf-segTree[i].lf;segTree[i].lenTwice=segTree[i].rf-segTree[i].lf;return;}else if(segTree[i].cnt==1)//被线段覆盖次数==1
    {segTree[i].lenOnce=segTree[i].rf-segTree[i].lf;if(segTree[i].l+1==segTree[i].r)segTree[i].lenTwice=0;//当是叶子节点时else segTree[i].lenTwice=segTree[i<<1].lenOnce+segTree[(i<<1)|1].lenOnce;//因为当前线段被覆盖过1次,所以覆盖2次以上的长度是其左右孩子的覆盖1次以上的长度的和
    }else//被线段覆盖次数==0
    {if(segTree[i].l+1==segTree[i].r)//当是叶子节点时
        {segTree[i].lenOnce=segTree[i].lenTwice=0;}else{//因为当前节点被覆盖次数为0,所以lenOnce和lenTwice的值由左右孩子决定segTree[i].lenOnce=segTree[i<<1].lenOnce+segTree[(i<<1)|1].lenOnce;segTree[i].lenTwice=segTree[i<<1].lenTwice+segTree[(i<<1)|1].lenTwice;}}
}
void update(int i,Line e)
{if(e.y1==segTree[i].lf&&segTree[i].rf==e.y2){segTree[i].cnt+=e.f;calen(i);//更新当前节点信息return;}if(e.y2<=segTree[i<<1].rf) update(i<<1,e);else if(e.y1>=segTree[(i<<1)|1].lf) update((i<<1)|1,e);else{Line temp=e;temp.y2=segTree[i<<1].rf;update(i<<1,temp);temp=e;temp.y1=segTree[(i<<1)|1].lf;update((i<<1)|1,temp);}calen(i);//更新当前节点信息
}
int main()
{//freopen("in.txt","r",stdin);//freopen("out.txt","w",stdout);int T;int n;double x1,y1,x2,y2;scanf("%d",&T);while(T--){scanf("%d",&n);int t=1;for(int i=1; i<=n; i++){scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);//这里题目描述有问题?左下角和右上角line[t].y1=y1;line[t].y2=y2;line[t].x=x1;line[t].f=1;//左边的边y[t]=y1;//y坐标离散化(下端点)t++;line[t].y1=y1;line[t].y2=y2;line[t].x=x2;line[t].f=-1;//右边的边y[t]=y2;//y坐标离散化(上端点)t++;}sort(line+1,line+t,cmp);//线段按x的值升序排序sort(y+1,y+t);//y坐标从小下到大排序Build(1,1,t-1);//用离散化后的y值(也就是t)建立线段树
update(1,line[1]);//从第一条扫描线向右扫描double ans=0;for(int i=2; i<t; i++){ans+=segTree[1].lenTwice*(line[i].x-line[i-1].x);update(1,line[i]);}printf("%.2lf\n",ans);}return 0;
}

View Code

转载于:https://www.cnblogs.com/bofengyu/p/4964641.html

HDU - 1255 覆盖的面积(线段树求矩形面积交 扫描线+离散化)相关推荐

  1. 线段树求矩形面积并 扫描线+离散化

    顾名思义,扫描法就是用一根想象中的线扫过所有矩形,在写代码的过程中,这根线很重要.方向的话,可以左右扫,也可以上下扫.方法是一样的,这里我用的是由下向上的扫描法. 如上图所示,坐标系内有两个矩形.位置 ...

  2. hdu 1255(线段树求重叠面积)

    扫描线求矩形重叠面积:http://www.cnblogs.com/scau20110726/archive/2013/04/12/3016765.html http://www.tuicool.co ...

  3. 线段树 求矩形并 面积

    http://www.cnblogs.com/ka200812/archive/2011/11/13/2247064.html 给定一个矩形的左下角坐标和右上角坐标分别为:(x1,y1).(x2,y2 ...

  4. hdu 1255 覆盖的面积 (扫描线求矩形交)

    覆盖的面积 Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Sub ...

  5. hdu 2461(线段树求面积并)

    题意:给出N个矩形,M次询问 每次询问给出R个,问这R个矩形围成的面积 解题思路:对于每次询问,做一次线段树求面积的并操作. 每个节点保存的信息有l,r,cover,len分别表示该节点表示的区间[l ...

  6. poj 1151(线段树求面积并)

    解题思路:线段树求面积并,水题 #include<iostream> #include<cstdio> #include<cstring> #include< ...

  7. hdu 1255 覆盖的面积

    覆盖的面积 Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total Sub ...

  8. 2019CCPC网络赛 1002 HDU 6703(权值线段树)

    2019CCPC网络赛 1002 HDU 6703(权值线段树) 思路:用权值线段树存题目给的数据后,2操作就是求权值线段树中大于等于k的部分中,靠近左端点的第一个大于r的值(这个求出来的只是原序列中 ...

  9. 线段树求逆序对(hdu1394Minimum Inversion Number)

    说实话,线段树求逆序对我理解了半天诶,不知是否有人像我一样. 对于每个数来说,只有和已经出现过的.比它大的数才能形成逆序对,那么在给定的数列中,每给一个数就向前找比它大的数. 样例:10 1 3 6 ...

最新文章

  1. vector机器人 WHERE TO USE VECTOR 在哪里使用 VECTOR
  2. 模糊数学笔记:六、模糊模型识别-II(择近原则)
  3. oracle 导入流程,Oracle数据库(旧)空间数据表创建与数据导入流程说明V2.0
  4. postgresql 创建用户_国内源安装postgresql
  5. 【转】最牛B的编码套路
  6. 八款Js框架介绍及比较~转载
  7. 问题解决: 此文件来自其他计算机,可能被阻止以帮助保护该计算机
  8. git 怎么查看合并过来哪些代码_git整理纷乱的历史合并记录
  9. mysql命令大全_Keeplived+mysql双master高可用如何实现?
  10. java---servlet与filter的联系与区别
  11. BT4中文版(集成spoonwep2/spoonwpa) 破解无线网卡非常容易
  12. 论文参考文献格式写法
  13. symfony配置文件之三:routing和路由配置
  14. CAD中插入外部参照字体会变繁体_CAD中形形色色的“替换”技巧你确定都知道了?...
  15. 徒手写代码之《机器学习实战》-----决策树算法(2)(使用决策树预测隐形眼镜类型)
  16. FL STUDIO 21中文版已经发布啦,有什么新功能?
  17. mybatis 自动填充无效_mybatisPlus踩坑之--自动填充
  18. MySQL 5.6.21下载安装之安装篇(二)
  19. 用伪造的TCP协议头花式欺骗核心转发设备?
  20. Linux系列:二进制 Byte为什么是-127—127范围

热门文章

  1. 第十一周项目实践4 BFS(广度优先搜索)基本模板
  2. [YTU]_2641 9 填空题:静态成员---计算学生个数)
  3. 1.2.2一个数可以有多少种用连续素数之和表示POJ 2739
  4. 判断二进制数1的个数
  5. apache mysql php实现最大负载的方法_如何架设高性能nginx+php+mysql搭配的服务器,解决高并发问题...
  6. 五种排序方式gif展示【python】
  7. STL中的栈结构和队列结构
  8. ElasticSearch启动报错,bootstrap checks failed
  9. 写爬虫,不会正则怎么行?
  10. Freemarker 内置函数 数字、字符串、日期格式化用法介绍