BZOJ3076 & 洛谷3081:[USACO2013 MAR]Hill Walk 山走——题解
http://www.lydsy.com/JudgeOnline/problem.php?id=3076
https://www.luogu.org/problemnew/show/P3081#sub
有N(1 <= N <= 100,000)座小山,每座山所占的区域用直线(x1, y1) 到 (x2, y2)来表示(x1 < x2 并且 y1 < y2)。也就是说这些山用笛卡尔坐标系里的线段来表示,这些用于表示小山的线段都没有任何交点,第一座山的一端位于(x1, y1) = (0,0)
贝西从(0,0)开始在第一座山上漫步,一旦贝西到了一座山,她会一直走到该山的终点,这时,她会从边缘处起跳,如果她降落到另一座山上,她会继续在新的山上漫步。贝西起跳后沿y轴方向下落,如果贝西不能降落到一座山上,她会一直下落,直到到达y轴的负无穷大位置(y = -infinity)。
每座用线段表示的山 (x1, y1) -> (x2, y2)包含(x1, y1)这个点,但不包含(x2, y2) ,请计算出贝西总共在多少座山上漫步了。
……这道题网上几乎没找到题解……所以是参考了标程的(我也看不懂英文啊……)
首先说一下解题思路:我们显然要走一遍线段,并且在掉落的时候以O(logn)的速度查找到你应当掉落在那个线段上。那么复杂度就是O(nlogn)了。
显然的思路,我们记录当前走在了哪个线段,用set维护接下来**可能**走的线段,取出我们所要的线段即可。
(PS:“可能”指你所在的坐标的x在某个线段的x区间内)
走在了哪个线段没有难度,本代码的cur就是做这个的。
找接下来可能走的线段(并且及时剔除掉走过/不能走到的线段),扫描线即可解决。
那么这题的难点就在于我们如何排序(也就是set的比较方法)来取出我们所要的线段。
先画个图感性理解一下,黑点是我们所在的位置,编号就是我们排好序之后的编号。
(我们所在的位置是3,显然只需要跳到3-1=2的编号的那条边,我们就实现了走的动作。)
设比较的两条线段a,b,则分为两种情况:
1.ax2<bx2(eg:a=5,b=4或a=2,b=4):
我们把a的右端点和b的左端点连线,显然发现这条线段将这两条线划分开来,上面的线编号大,下面的线编号小。
2.ax2>bx2(eg:a=1,b=3或a=4,b=3):
同1的思路,其实只是把a和b颠倒而已。
我们可以通过斜率来比较两条线的位置。
(PPS:貌似可以两条线分别取一个不同方向的点连线貌似就可以,不必非得按照题解的连线方法)
(PPPS:可能有这种情况(举一个1的例子):
这时候我们只需要平移一下就好了。
所以用斜率判断是最准确的。)
————————————————————
一些小细节,比如说你已经在最下层线段没法在往下走时及时跳出
#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=1e5+;
inline int read(){
int X=,w=;char ch=;
while(ch<''||ch>''){if(ch=='-')w=-;ch=getchar();}
while(ch>=''&&ch<='')X=(X<<)+(X<<)+ch-'',ch=getchar();
return X*w;
}
struct edge{
int x1,x2,y1,y2,id;
bool operator <(edge const& a)const{
if(x2<a.x2){
return (ll)(y2-a.y1)*(ll)(a.x2-a.x1)<(ll)(a.y2-a.y1)*(ll)(x2-a.x1);
}else{
return (ll)(a.y2-y1)*(ll)(x2-x1)>(ll)(y2-y1)*(ll)(a.x2-x1);
}
}
}e[N];
struct event{
int x,y,id;
bool operator <(event const& a)const{
return x<a.x||(x==a.x&&y<a.y);
}
}f[*N];
set<edge>s;
set<edge>::iterator it1;
set<edge>::iterator it2;
int main(){
int n=read();
for(int i=;i<=n;i++){
e[i].x1=read();e[i].y1=read();
e[i].x2=read();e[i].y2=read();
e[i].id=i;
f[i*-].x=e[i].x1;f[i*].x=e[i].x2;
f[i*-].y=e[i].y1;f[i*].y=e[i].y2;
f[i*-].id=f[i*].id=i;
}
sort(f+,f+*n+);
s.insert(e[]);
int cur=,tot=;
for(int i=;i<=*n;i++){
event ev=f[i];
edge ed=e[ev.id];
if(ev.x==ed.x1)s.insert(ed);
else if(ev.id==cur){
it1=s.find(ed);
if(it1==s.begin())break;
it2=it1;--it2;
cur=it2->id;
s.erase(it1);
tot++;
}else s.erase(ed);
}
printf("%d\n",tot);
return ;
}
+++++++++++++++++++++++++++++++++++++++++++
+本文作者:luyouqi233。 +
+欢迎访问我的博客:http://www.cnblogs.com/luyouqi233/+
+++++++++++++++++++++++++++++++++++++++++++
BZOJ3076 & 洛谷3081:[USACO2013 MAR]Hill Walk 山走——题解的更多相关文章
- 【洛谷P3369】【模板】普通平衡树题解
[洛谷P3369][模板]普通平衡树题解 题目链接 题意: 您需要写一种数据结构(可参考题目标题),来维护一些数,其中需要提供以下操作:1. 插入x数2. 删除x数(若有多个相同的数,因只删除一个)3 ...
- 洛谷 P1992 不想兜圈的老爷爷 题解
洛谷 P1992 不想兜圈的老爷爷 题解 题目描述 一位年过古稀的老爷爷在乡间行走 而他不想兜圈子 因为那会使他昏沉 偶然路过小A发扬助人为乐优良传统 带上地图 想知道路况是否一定使他清醒 usqwe ...
- BZOJ5285 & 洛谷4424 & UOJ384:[HNOI/AHOI2018]寻宝游戏——题解
https://www.lydsy.com/JudgeOnline/problem.php?id=5285 https://www.luogu.org/problemnew/show/P4424 ht ...
- 洛谷p3384【模板】树链剖分题解
洛谷p3384 [模板]树链剖分错误记录 首先感谢\(lfd\)在课上调了出来\(Orz\) \(1\).以后少写全局变量 \(2\).线段树递归的时候最好把左右区间一起传 \(3\).写\(dfs\ ...
- 洛谷 P1789 【Mc生存】插火把 题解
P1789 [Mc生存]插火把 题目背景 初一党应该都知道...... 题目描述 话说有一天 linyorson 在"我的世界"开了一个 \(n\times n(n\le 100) ...
- 洛谷P1066 2^k进制数(题解)(递推版)
https://www.luogu.org/problemnew/show/P1066(题目传送) (题解)https://www.luogu.org/problemnew/solution/P106 ...
- 洛谷P1373小a和uim大逃离题解
题目 这个题好坑啊,首先是他会卡空间,然后我们就只能把一种比较好理解的状态给舍弃,因为空间开不下,然而采用一种难理解的状态就是\(dp[i][j][l][0/1]\)表示\(i\),\(j\)位置,两 ...
- 洛谷 P1202 [USACO1.1]黑色星期五Friday the Thirteenth 题解
题目传送门 这道题暴力就能解决. #include<bits/stdc++.h> using namespace std; int xi; ,ans[]; int main() { int ...
- 洛谷 P2241统计方形(数据加强版) 题解
题目传送门 说是加强版,其实可以把棋盘那道题的代码粘过来(注意要开long long): #include<bits/stdc++.h> using namespace std; ,c; ...
随机推荐
- SpringBoot学习:整合shiro(验证码功能和登录次数限制功能)
项目下载地址:http://download.csdn.NET/detail/aqsunkai/9805821 (一)验证码 首先login.jsp里增加了获取验证码图片的标签: <body s ...
- 使用PowerDesign15反向生成数据库
在Pd15中建立物理模型后,可以通过反向工程直接生成数据库的表结构.主要有以下几个步骤: 1. 首先设置一下数据库配置,选择对应要使用的数据库(此处选择Sql Server 2008 R ...
- jenkins--Jenkins+Git+coding+maven 实现自动化测试持续集成
1.打开Jenkins官网,下载jenkins.war https://jenkins.io/download/ 2.将该war包直接放置到Tomcat的webapp下. 3.查看自己Tomcat的端 ...
- python基础之变量和简单数据类型
1.1 变量的命名和使用规范 变量名可以包含数字.字母.下划线,但是不能以数字开头. 变量名不能包含空格,可使用下划线来分割其中的单词. 不要将Python关键字和函数名用作变量名. 变量名应既简短又 ...
- unity发布自定义分辨率
如果你需要发布unity时想要使用自己设置的分辨率仅需要一下几个步骤: 打开Build Setting->PlayerSetting->Resolution and Presentatio ...
- Anyproxy抓包工具
1.安装Nodejs和AnyProxy以及安装模拟器(移动端抓包工具) 第一步:nodejs下载地址:http://nodejs.cn/download/ 下载Windows版本,直接运行安装即可, ...
- mapReduce入门教程
什么是MapReduce MapReduce是Google提出的一个软件架构,用于大规模数据集(大于1TB)的并行运算.概念"Map(映射)"和"Reduce(归纳)&q ...
- opencv-学习笔记(4)-模糊
opencv-学习笔记(4)-模糊 本章要点: 4种模糊方式 2d卷积 Cv2.filter2D(‘图像对象’,‘目标图像这里直接设为-1即可’,kernal,anchor(-1,-1)) 一般后一个 ...
- facebook演讲
任何为了更大愿景工作的人,可能会被称为疯子,即使你最终获得成功. 任何为了复杂问题工作的人,都会因为不能全面了解挑战而被指责,即使你不可能事先了解一切. 任何抓住主动权先行一步的人,都会因为步子太快而 ...
- ZOJ 2532 Internship(最大流找关键割边)
Description CIA headquarter collects data from across the country through its classified network. Th ...