[HNOI2012]射箭

题目描述：

有n(<=10w)个靶子，射一条过原点的抛物线，问从靶子1起最多能穿过多少连续标号的靶子。

输入：各种靶子。

输出：最多靶子数。

非常裸的一道半平面交的题目，二分答案判定即可。

不过半平面交的算法倒是YY了好久，一直纠结于各种细节……不过写过一遍之后感觉还是比较简单的。

值得注意的有两点：

1.初始的矩形框有方向。（这题虽然不用）  

2.直线加入完毕后还要判断L1切队尾的情况。

另外，这题精度很坑很坑，所以不要乱用eps。

计算几何的题目最好套个模板，写起来顺手多了。

Code：

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int SN=200100;
const double oo=1e20, pi=acos(-1), eps=1e-12;

struct dbl { double x,y; };
struct line { dbl P,v; double jj; int od; };
int n,m,i,l,r,md,ql,qr;
line R[SN],que[SN];
dbl p[SN];

dbl operator + (dbl a, dbl b) { return (dbl){a.x+b.x, a.y+b.y}; }
dbl operator - (dbl a, dbl b) { return (dbl){a.x-b.x, a.y-b.y}; }
dbl operator * (dbl a, double b) { return (dbl){a.x*b, a.y*b}; }
dbl operator / (dbl a, double b) { return (dbl){a.x/b, a.y/b}; }

int dcmp(double d) { if(d<-eps) return -1; else return d>eps; }
double cj(dbl a, dbl b) { return a.x*b.y-a.y*b.x; }
double dj(dbl a, dbl b) { return a.x*b.x+a.y*b.y; }
double dis(dbl a, dbl b) { return sqrt( (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y) ); }
double fjj(dbl v) { double rn=atan2(v.y, v.x); rn=(rn<0)?rn+2*pi:rn; return rn; }
void rotate(dbl& v, double g) { double C=cos(g),S=sin(g); v=(dbl){C*v.x-S*v.y, S*v.x+C*v.y}; }

bool Cmp(line a, line b) { return a.jj<b.jj; }
bool OnRight(dbl p, line L) { return dcmp(cj(L.v, p-L.P))<0; }
dbl LLjd(line a, line b) { dbl u=b.P-a.P; return a.P+a.v*(cj(u, b.v)/cj(a.v, b.v)); }

void AddL(double k, double b, bool ck)
{
  m++, R[m].od=i, R[m].P=(dbl){0, b}, R[m].v=(dbl){1, k};
  if(ck) R[m].v=(dbl){-1, -k};
}
bool Chk(int n)
{
  ql=1,qr=0;
  for(i=1; i<=m; i++) if(R[i].od<=n) 
  {
    while(ql<qr && OnRight(p[qr-1], R[i])) qr--;
    while(ql<qr && OnRight(p[ ql ], R[i])) ql++;
    que[++qr]=R[i]; 
    if(ql<qr && dcmp(que[qr].jj-que[qr-1].jj)==0) qr--;
    if(ql<qr) p[qr-1]=LLjd(que[qr], que[qr-1]);
  }
  while(ql<qr && OnRight(p[qr-1], que[ql])) qr--;
  return qr-ql>1;
}

int main()
{
  double x,y1,y2;
  freopen("oj.in", "r", stdin);
  freopen("oj.out", "w", stdout);
  for(i=1,scanf("%d", &n); i<=n; i++) 
  {
    scanf("%lf%lf%lf", &x, &y1, &y2);
    AddL(-x, y1/x, 0); AddL(-x, y2/x, 1);
  }
  m++, R[m].P=(dbl){-1, -1}, R[m].v=(dbl){1, 0};  m++, R[m].P=(dbl){oo, -1}, R[m].v=(dbl){0, 1};
  m++, R[m].P=(dbl){oo, oo}, R[m].v=(dbl){-1, 0}; m++, R[m].P=(dbl){-1, oo}, R[m].v=(dbl){0, -1};
  for(i=1; i<=m; i++) R[i].jj=fjj(R[i].v);
  sort(R+1, R+m+1, Cmp);
  for(l=1,r=n,md=(l+r+1)/2; l<r; md=(l+r+1)/2)
    if(Chk(md)) l=md;
    else r=md-1;
  printf("%d\n", l);
  return 0;
}