[bzoj2732][半平面交]射箭

2732: [HNOI2012]射箭

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 3709 Solved: 1146
[Submit][Status][Discuss]
Description

沫沫最近在玩一个二维的射箭游戏，如下图 1 所示，这个游戏中的 x 轴在地面，第一象限中有一些竖直线段作为靶子，任意两个靶子都没有公共部分，也不会接触坐标轴。沫沫控制一个位于(0,0)的弓箭手，可以朝 0 至 90?中的任意角度（不包括 0度和 90度），以任意大小的力量射出带有穿透能力的光之箭。由于游戏中没有空气阻力，并且光之箭没有箭身，箭的轨迹会是一条标准的抛物线，被轨迹穿过的所有靶子都认为被沫沫射中了，包括那些 只有端点被射中的靶子。这个游戏有多种模式，其中沫沫最喜欢的是闯关模式。在闯关模式中，第一关只有一个靶 子，射中这个靶子即可进入第二关，这时在第一关的基础上会出现另外一个靶子，若能够一箭 双雕射中这两个靶子便可进入第三关，这时会出现第三个靶子。依此类推，每过一关都会新出 现一个靶子，在第 K 关必须一箭射中前 K 关出现的所有 K 个靶子才能进入第 K+1 关，否则游戏 结束。沫沫花了很多时间在这个游戏上，却最多只能玩到第七关“七星连珠”，这让她非常困惑。 于是她设法获得了每一关出现的靶子的位置，想让你告诉她，最多能通过多少关

Input

输入文件第一行是一个正整数N，表示一共有N关。接下来有N行，第i+1行是用空格隔开的三个正整数xi，yi1，yi2(yi1

sol:

显然是二分题，然后判断是否有a,b使得限制成立。我们考虑题目其实是给定了若干个限制y1<=ax^2+bx<=y2
有二次的，有点麻烦。但我们发现x是已知的。我们移项一下
x^2a+xb-y1>=0
x^2a+xb-y2<=0
把x^2看成整体这就是一个半平面的形式了。求半平面交即可。
但是注意到a<0和b>0，在外面加个大框即可。而且记得不要把这个框放在x=0,y=0的地方，要稍微偏一点，我是- 1e-10,我开-1e-7的时候wa掉了，大家注意一下精度。

#include<iostream>
#include<algorithm>
#include<string>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<cmath>
using namespace std;
typedef long long ll;
typedef double s64;
int m;
inline int read()
{
    char c;
    int res,flag=0;
    while((c=getchar())>'9'||c<'0') if(c=='-')flag=1;
    res=c-'0';
    while((c=getchar())>='0'&&c<='9') res=(res<<3)+(res<<1)+c-'0';
    return flag?-res:res;
}
const int N=210000;
s64 a[N][3];
const s64 eps=1e-10;
struct P
{
    s64 x,y;
    friend P operator +(P a,P b) {return (P){a.x+b.x,a.y+b.y};}
    friend P operator -(P a,P b) {return (P){a.x-b.x,a.y-b.y};}
    friend P operator *(P a,s64 b) {return (P){a.x*b,a.y*b};}
    friend s64 operator *(P a,P b) {return a.x*b.y-a.y*b.x;}
    friend s64 operator /(P a,P b) {return a.x*b.x+a.y*b.y;}
};
struct L
{
    P a,b,v;
    s64 ang;
    friend bool operator <(L a,L b) {return a.ang<b.ang||a.ang==b.ang&&a.v*(b.b-a.a)>0;}
    friend inline P inter(L a,L b)
    {
        P nw=b.a-a.a;
        s64 tt=(nw*a.v)/(a.v*b.v);
        return b.a+b.v*tt;
    }
    friend inline bool jud(P a,L b) {return b.v*(a-b.a)<0;}
}q[N],l[N],p[N];
inline s64 sqr(s64 x) {return x*x;}
const s64 inf=1e9+7;
int num;
inline bool check(int n)
{
    n=4+n*2;
    int tot=0;
    for(int i=1;i<=n;++i) l[i]=p[i];
    for(int i=1;i<=n;++i) l[i].ang=atan2(l[i].v.y,l[i].v.x);
    sort(l+1,l+1+n);
    tot=0;
    for(int i=1;i<=n;++i)
    if(l[i].ang!=l[i-1].ang) l[++tot]=l[i];
    n=tot;
    q[1]=l[1];
    q[2]=l[2];
    int R=2;
    int L=1;
    for(int i=3;i<=n;++i)
    {
        while(L<R&&jud(inter(q[R-1],q[R]),l[i])) --R;
        while(L<R&&jud(inter(q[L],q[L+1]),l[i])) ++L;
        q[++R]=l[i];
    }
    while(L<R&&jud(inter(q[R],q[R-1]),q[L])) --R;
    while(L<R&&jud(inter(q[L+1],q[L]),q[R])) ++L;
    if(R-L+1>=3) return 1;
    return 0;
}
int main()
{
//  freopen("2732.in","r",stdin);
//  freopen(".out","w",stdout);
    m=read();
    p[1].a.x=0;p[1].a.y=eps;
    p[1].v.x=1;p[1].v.y=0;
    p[1].b=p[1].a+p[1].v;
    p[2].a.x=-eps;p[2].a.y=0;
    p[2].v.x=0;p[2].v.y=1;
    p[2].b=p[2].a+p[2].v;
    p[3].a.x=0;p[3].a.y=inf;
    p[3].b.x=-inf;p[3].b.y=inf;
    p[3].v=p[3].b-p[3].a;
    p[4].a.x=-inf;p[4].a.y=inf;
    p[4].b.x=-inf;p[4].b.y=0;
    p[4].v=p[4].b-p[4].a;
    num=4;
    for(int i=1;i<=m;++i)
    {
        scanf("%lf%lf%lf",&a[i][0],&a[i][1],&a[i][2]);
        a[i][1]-=eps;
        a[i][2]+=eps;
        ++num;
        p[num].a=(P){0,a[i][1]/a[i][0]};
        p[num].v=(P){a[i][0],-sqr(a[i][0])};
        p[num].b=p[num].a+p[num].v;
        ++num;
        p[num].a=(P){0,a[i][2]/a[i][0]};
        p[num].v=(P){-a[i][0],sqr(a[i][0])};
        p[num].b=p[num].a+p[num].v;
    }
    int l,r,mid;
    l=1;r=m;
    while(l+1<r)
    {
        mid=l+r>>1;
        if(check(mid)) l=mid;
        else r=mid;
    }
    if(check(r)) printf("%d\n",r);
    else printf("%d\n",l);
}