Description
沫沫最近在玩一个二维的射箭游戏,如下图 1 所示,这个游戏中的 x 轴在地面,第一象限中有一些竖直线段作为靶子,任意两个靶子都没有公共部分,也不会接触坐标轴。沫沫控制一个位于(0,0)的弓箭手,可以朝 0 至 90?中的任意角度(不包括 0度和 90度),以任意大小的力量射出带有穿透能力的光之箭。由于游戏中没有空气阻力,并且光之箭没有箭身,箭的轨迹会是一条标准的抛物线,被轨迹穿过的所有靶子都认为被沫沫射中了,包括那些 只有端点被射中的靶子。这个游戏有多种模式,其中沫沫最喜欢的是闯关模式。在闯关模式中,第一关只有一个靶 子,射中这个靶子即可进入第二关,这时在第一关的基础上会出现另外一个靶子,若能够一箭 双雕射中这两个靶子便可进入第三关,这时会出现第三个靶子。依此类推,每过一关都会新出 现一个靶子,在第 K 关必须一箭射中前 K 关出现的所有 K 个靶子才能进入第 K+1 关,否则游戏 结束。沫沫花了很多时间在这个游戏上,却最多只能玩到第七关“七星连珠”,这让她非常困惑。 于是她设法获得了每一关出现的靶子的位置,想让你告诉她,最多能通过多少关
Input
输入文件第一行是一个正整数N,表示一共有N关。接下来有N行,第i+1行是用空格隔开的三个正整数xi,yi1,yi2(yi1< yi2 ),表示第i关出现的靶子的横坐标是xi,纵坐标的范围是从yi1到yi2 。
输入保证30%的数据满足N≤100,50%的数据满足N≤5000,100%的数据满足N≤100000且给 出的所有坐标不超过109 。
Output
仅包含一个整数,表示最多的通关数。
Sample Input
5
2 8 12
5 4 5
3 8 10
6 2 3
1 3 7
Sample Output
3
HINT
题解
先感谢这位dalao的debug和数据:
http://www.cnblogs.com/liu-runda/p/6439792.html
不放先设我们的抛物线为
y=ax2+bx
那么对于一个障碍我们可以列出这样的不等式:
变形一下就可以得到两个不等式:
每个障碍都可以两个不等式,这样的话二分答案mid判断前mid个障碍产生的2mid个不等式是否有解就好了。这样就是一个轻松愉快的半平面交了不是吗?
然而WA,接连的WA。
首先呢,double是行不通的(我也不知道为什么),反正就用long double就对了。
然后呢,考虑到这可是射出去的箭矢,显然要满足 {a<0b>0 ,那么我们在半平面交里的包围框就要调整成包围了整个第二象限,而且不能包围坐标轴(这个用-eps代替0就好了)。
嗯,反正我是把这两个处理一下以后就A了。
#include<bits/stdc++.h>
using namespace std;
inline int read(){
int x = 0, f = 1; char c = getchar();
while(!isdigit(c)) { if(c == '-') f = -1; c = getchar(); }
while(isdigit(c)) { x = x * 10 + c - '0'; c = getchar(); }
return x * f;
}
#define db long double
#define inf 1e15
#define eps 1e-18
const int N = 200000 + 10;
int n, m, hd, tl;
struct pt{
db x, y;
pt(){}
pt(db a, db b){x = a, y = b;}
pt operator - (const pt &a) const{
return pt(x - a.x, y - a.y);
}
pt operator + (const pt &a) const{
return pt(x + a.x, y + a.y);
}
pt operator * (const db k) const{
return pt(x * k, y * k);
}
db operator * (const pt &a) const{
return x * a.y - y * a.x;
}
};
struct ln{
pt a, b;
db ang;
int id;
bool operator < (const ln &rhs) const{
return ang < rhs.ang || ang == rhs.ang && (b - a) * (rhs.a - a) > 0;
}
}l[N], a[N], q[N];
pt inter(ln x, ln y){
pt u = x.a - y.a, v1 = x.b - x.a, v2 = y.b - y.a;
db t = (v2 * u) / (v1 * v2);
return x.a + v1 * t;
}
bool jud(ln a, ln b, ln t){
pt p = inter(a, b);
return (p - t.a) * (t.b - t.a) > 0;
}
db cal(db a, db b, db x){return b / a - a * x;}
void init(){
scanf("%d", &n);
l[++m].a = pt(-1e-13, 1e-13), l[m].b = pt(-1e-13, inf);
l[++m].a = pt(-1e-13, inf), l[m].b = pt(-inf, inf);
l[++m].a = pt(-inf, inf), l[m].b = pt(-inf, 1e-13);
l[++m].a = pt(-inf, 1e-13), l[m].b = pt(-1e-13, 1e-13);
db x, ya, yb;
for(int i = 1; i <= n; i++){
x = read(), ya = read(), yb = read();
l[++m].a.x = -1; l[m].a.y = cal(x, ya, -1);
l[m].b.x = 1; l[m].b.y = cal(x, ya, 1);
l[++m].a.x = 1; l[m].a.y = cal(x, yb, 1);
l[m].b.x = -1; l[m].b.y = cal(x, yb, -1);
l[m].id = l[m-1].id = i;
}
for(int i = 1; i <= m; i++)
l[i].ang = atan2(l[i].b.y - l[i].a.y, l[i].b.x - l[i].a.x);
sort(l + 1, l + m + 1);
}
bool hpi(int k){
int tot = 0;
for(int i = 1; i <= m; i++)
if(l[i].id <= k){
if(l[i].ang != a[tot].ang) tot++;
a[tot] = l[i];
}
hd = 1; tl = 0;
q[++tl] = a[1]; q[++tl] = a[2];
for(int i = 3; i <= tot; i++){
while(hd < tl && jud(q[tl-1], q[tl], a[i])) tl--;
while(hd < tl && jud(q[hd+1], q[hd], a[i])) hd++;
q[++tl] = a[i];
}
while(hd < tl && jud(q[tl-1], q[tl], q[hd])) tl--;
while(hd < tl && jud(q[hd+1], q[hd], q[tl])) hd++;
}
void work(){
int l = 1, r = n, ans;
while(l <= r){
int mid = (l + r) >> 1;
hpi(mid);
if(tl - hd >= 2) ans = mid, l = mid + 1;
else r = mid - 1;
}
printf("%d\n", ans);
}
int main(){
init();
work();
return 0;
}