poj1177
题意:在平面直角坐标系内给出一些与坐标轴平行的矩形,将这些矩形覆盖的区域求并集,然后问这个区域的周长是多少。(边与边重合的地方不计入周长)
分析:线段树。曾经做过类似的求矩形覆盖的总面积的题。这道题同样要使用扫描线算法。属于线保留型线段树。
我们先领扫描线与y轴平行。
线段树内每个节点除了要记录该区间被覆盖了几层之外,还要记录当前状态下扫描线在该区间(开区间)内与多少条与x轴平行的边相交。
节点上还有两个bool型变量,记录该区间内(包括子树)线段覆盖是否接触到该区间的起始和结束点。
在父节点如果没被整个覆盖,则需要从子区间的起始和结束点来更新父节点两端点的覆盖情况。
更新过程在返回时,父节点的交点数量应等于两子节点交点数量的和,另外特判一下两子区间的公共点(父节点的中点),判断这里是不是覆盖与未覆盖的分界点,如果是则还需要在父节点上增加这个交点。
这样就得知了每段与x轴平行的距离内有多少条线段需要计算,距离乘以数量即可。这样就计算出了所有与x轴平行的周长上的边的总长度。
之后让扫描线与x轴平行即可计算出与y轴平行的所有周长上的边的总长度。
#include <cstdio>
#include <algorithm>
using namespace std;
//scanning from left to right
//discretionize Ys
#define MAX_REC_NUM 5005
#define MAX_INTERVAL MAX_REC_NUM * 2 struct Node
{
int l, r;
Node *pleft, *pright;
int num;
bool to_left, to_right;
int edge_num;
}; int node_cnt;
Node tree[MAX_INTERVAL * ]; struct Interval
{
int start, end;
int pos;
int value;
Interval()
{}
Interval(int start, int end, int pos, int value):start(start), end(end), pos(pos), value(value)
{}
bool operator < (const Interval &a)const
{
if (pos != a.pos)
return pos < a.pos;
return value > a.value;
}
}interval[MAX_REC_NUM * ]; struct Rectangle
{
int l, d, u, r;
}rec[MAX_REC_NUM]; int discrete[MAX_REC_NUM * ];
int discrete_num;
int rec_num;
int interval_num; int get_index(int a)
{
return lower_bound(discrete, discrete + discrete_num, a) - discrete;
} void discretization(int discrete[], int &discrete_num)
{
sort(discrete, discrete + discrete_num);
discrete_num = unique(discrete, discrete + discrete_num) - discrete;
} void input()
{
scanf("%d", &rec_num);
for (int i = ; i < rec_num; i++)
{
int l, d, r, u;
scanf("%d%d%d%d", &l, &d, &r, &u);
rec[i].l = l;
rec[i].r = r;
rec[i].u = u;
rec[i].d = d;
}
} void make_xscan()
{
discrete_num = ;
interval_num = ;
for (int i = ; i < rec_num; i++)
{
int l, d, r, u;
l = rec[i].l;
r = rec[i].r;
u = rec[i].u;
d = rec[i].d;
interval[interval_num++] = Interval(d, u, l, );
interval[interval_num++] = Interval(d, u, r, -);
discrete[discrete_num++] = u;
discrete[discrete_num++] = d;
}
} void make_yscan()
{
discrete_num = ;
interval_num = ;
for (int i = ; i < rec_num; i++)
{
int l, d, r, u;
l = rec[i].l;
r = rec[i].r;
u = rec[i].u;
d = rec[i].d;
interval[interval_num++] = Interval(l, r, d, );
interval[interval_num++] = Interval(l, r, u, -);
discrete[discrete_num++] = l;
discrete[discrete_num++] = r;
}
} void buildtree(Node *proot, int s, int e)
{
proot->l = s;
proot->r = e;
proot->to_left = false;
proot->to_right = false;
proot->num = ;
proot->edge_num = ;
if (s == e)
{
proot->pleft = proot->pright = NULL;
return;
}
node_cnt++;
proot->pleft = tree + node_cnt;
node_cnt++;
proot->pright = tree + node_cnt;
buildtree(proot->pleft, s, (s + e) / );
buildtree(proot->pright, (s + e) / + , e);
} void recount(Node *p)
{
if (p->num > )
{
p->edge_num = ;
p->to_right = p->to_left = true;
return;
}
if (p->pleft == NULL || p->pright == NULL)
{
p->edge_num = ;
p->to_right = p->to_left = false;
return;
}
p->to_left = p->pleft->to_left;
p->to_right = p->pright->to_right;
p->edge_num = p->pleft->edge_num + p->pright->edge_num;
if (p->pleft->to_right != p->pright->to_left)
p->edge_num++;
} void insert(Node *proot, int s, int e, int value)
{
if (s > proot->r || e < proot->l)
return;
s = max(s, proot->l);
e = min(e, proot->r);
if (s == proot->l && e == proot->r)
{
proot->num += value;
recount(proot);
return;
}
insert(proot->pleft, s, e, value);
insert(proot->pright, s, e, value);
recount(proot);
} long long work()
{
long long ans = ;
for (int i = ; i < interval_num; i++)
{
int s = get_index(interval[i].start);
int e = get_index(interval[i].end) - ;
insert(tree, s, e, interval[i].value);
long long line_num = tree->edge_num;
if (tree->to_left)
line_num++;
if (tree->to_right)
line_num++;
if (i != interval_num - )
ans += (interval[i + ].pos - interval[i].pos) * line_num;
}
return ans;
} int main()
{
input();
long long ans = ;
make_xscan();
sort(interval, interval + interval_num);
discretization(discrete, discrete_num);
buildtree(tree, , discrete_num);
ans += work(); make_yscan();
sort(interval, interval + interval_num);
discretization(discrete, discrete_num);
buildtree(tree, , discrete_num);
ans += work(); printf("%lld\n", ans);
return ;
}
poj1177的更多相关文章
- 【poj1177】 Picture
http://poj.org/problem?id=1177 (题目链接) 题意 求矩形周长并. Solution 转自:http://www.cnblogs.com/Booble/archive/2 ...
- [POJ1177]Picture
[POJ1177]Picture 试题描述 A number of rectangular posters, photographs and other pictures of the same sh ...
- 【HDOJ1828&&POJ1177】Picture(线段树,扫描线)
题意:给定n个矩形,求他们的并的周长 n<=5e3,abs(x[i])<=1e4 思路:From https://www.cnblogs.com/kuangbin/archive/2013 ...
- POJ1177 Picture —— 求矩形并的周长 线段树 + 扫描线 + 离散化
题目链接:https://vjudge.net/problem/POJ-1177 A number of rectangular posters, photographs and other pict ...
- 扫描线矩形周长的并 POJ1177
//扫描线矩形周长的并 POJ1177 // 我是按x轴 #include <iostream> #include <cstdio> #include <cstdlib& ...
- POJ1177 Picture 线段树+离散化+扫描线
求最终的覆盖图形周长,写这种代码应该短而精确,差的比较远 /* Problem: 1177 User: 96655 Memory: 348K Time: 32MS Language: C++ Resu ...
- POJ1177+线段树+扫描线
思路: 以y的值进行离散化 根据x的值 对每一条y轴边进行处理,如果是"左边"则插入,是"右边"则删除. /* 扫描线+线段树+离散化 求多个矩形的周长 */ ...
- POJ-1177 Picture 矩形覆盖周长并
题目链接:http://poj.org/problem?id=1177 比矩形面积并麻烦点,需要更新竖边的条数(平行于x轴扫描)..求横边的时候,保存上一个结果,加上当前长度与上一个结果差的绝对值就行 ...
- IOI1998 hdu1828 poj1177 Picture
写了一发扫描线竟然狂WA不止,hdu死活过不了,poj和当时IOI的数据(还花了我1dsdn积分..)都过了. 然后看到谋篇blog里有评论,把数据拿下来发现WA了. 数据是 20 0 1 11 0 ...
随机推荐
- Web启动服务器上的某一个服务
情景是这样的.. 网页打开一个数据列表..数据要求实时从其他多个平台上获取.. 所以就有了一个Web页面..还有个WinService的服务来定时获取这些数据... 问题来了.. 发现Service有 ...
- oracle vm virtualbox右ctrl切换显示模式
转自: http://blog.csdn.net/lyc_daniel/article/details/44195515 virtualbox里面有个HOME键,注意这个HOME键不一定是键盘上的HO ...
- glusterFS的常用命令 (转)
1. 启动/关闭/查看glusterd服务 # /etc/init.d/glusterd start # /etc/init.d/glusterd stop # /etc/init.d/g ...
- python简明手册学习
1.行末单独一个反斜杠表示字符串在下一行继续,而不是开始一个新的行. >>> "This is the first sentence.\ ... This is the s ...
- 【解决】同一url的http请求所获取的结果总是相同
曾经在WP7写过一个通过HTTP获取网页内容的小程序,当时一直没能够解决: 有一个网址,在每次点击刷新之后页面所呈现的内容都是不同的.但是进行HTTP请求时,结果将会一直重复. 从网上查资料得知,在请 ...
- Lexicography(数学推论>>求按字典序排第k个排列)
Lexicography Time Limit:1000MS Memory Limit:131072KB 64bit IO Format:%lld & %llu Submit ...
- Redis命令
redis的常用命令主要分为两个方面.一个是键值相关命令.一个是服务器相关命令(redis-cli进入终端) 1.键值相关命令 keys * 取出当前所有的key exists name 查看n是否有 ...
- lianjie
数值策划入门:如何确定游戏中的资源价值和定价http://bbs.gameres.com/thread_494366.html 一张常规的RPG游戏地图的制作流程 http://bbs.gameres ...
- Mac Pro 编译安装 PHP 5.6.21 及 问题汇总
[系统环境] 操作系统:OS X 10.11.5 Xcode:7.3.1 [注意] 编译之前,需要安装 xcode.Homebrew 套件! Mac Pro 安装 Homebrew 软件包管理工具 1 ...
- JS中的prototype(原文地址:http://www.cnblogs.com/yjf512/archive/2011/06/03/2071914.html)
JS中的phototype是JS中比较难理解的一个部分 本文基于下面几个知识点: 1 原型法设计模式 在.Net中可以使用clone()来实现原型法 原型法的主要思想是,现在有1个类A,我想要创建一个 ...