【bzoj2732】[HNOI2012]射箭

题目链接

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 = ax^2+bx$
那么对于一个障碍我们可以列出这样的不等式：

$$y_1\leq ax_1^2+bx_1\leq y_2$$
变形一下就可以得到两个不等式：
$$\begin{cases} b\geq y_1/x_1-ax_1\\ b\leq y_2/x_1-ax_1 \end{cases}$$
每个障碍都可以两个不等式，这样的话二分答案mid判断前mid个障碍产生的2mid个不等式是否有解就好了。这样就是一个轻松愉快的半平面交了不是吗？
然而WA，接连的WA。
为什么会变成这样呢……第一次写了半平面交，又写对了二分答案。两件快乐事情重合在一起。而这两份快乐，又给我带来更多的快乐。得到的，本该是像梦境一般幸福的时间……但是，为什么，会变成这样呢……
首先呢，double是行不通的（我也不知道为什么），反正就用long double就对了。
然后呢，考虑到这可是射出去的箭矢，显然要满足 $\begin{cases}a<0\\b>0\end{cases}$，那么我们在半平面交里的包围框就要调整成包围了整个第二象限，而且不能包围坐标轴（这个用-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;
}