BZOJ 2732 射箭

Posted 2020-08-02 GFY

http://www.lydsy.com/JudgeOnline/problem.php?id=2732

题意:给你n个靶子，让你求是否有一个经过原点的抛物线经过最多的前k个靶子，求出最大的k

思路:

就是这样的形式:y1<=ax^2+bx<=y2，这里y1,y2,x是已知的，有多组，我们发现，变量只有a和b了，而且都是一次，这样就转换成二元一次不等式组，这个问题就是高二学过的线性规划了

可以把它转换成半平面交，然后二分答案就可以了。

坑点:TM BZOJ上面要用longdouble，不然会WA一发，还有，考试的时候我居然在二分里面排序，导致复杂度退化到nlog^2 n T_T，这点要记住了。

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<algorithm>
 6 #define dou long double
 7 const dou inf=1e15;
 8 int tot,n;
 9 struct Point{
10     dou x,y;
11     Point(){}
12     Point(dou x0,dou y0):x(x0),y(y0){}
13 };
14 struct Line{
15     Point s,e;
16     dou slop;
17     int id;
18     Line(){}
19     Line(Point s0,Point e0):s(s0),e(e0){}
20 }l[200005],c[200005],L[200005];
21 int read(){
22     int t=0,f=1;char ch=getchar();
23     while (ch<‘0‘||ch>‘9‘){if (ch==‘-‘)f=-1;ch=getchar();}
24     while (‘0‘<=ch&&ch<=‘9‘){t=t*10+ch-‘0‘;ch=getchar();}
25     return t*f;
26 }
27 dou operator *(Point p1,Point p2){
28     return p1.x*p2.y-p1.y*p2.x;
29 }
30 Point operator -(Point p1,Point p2){
31     return Point(p1.x-p2.x,p1.y-p2.y);
32 }
33 bool cmp(Line p1,Line p2){
34     if (p1.slop==p2.slop) return (p2.e-p1.s)*(p1.e-p1.s)>=0;
35     else return p1.slop<p2.slop;
36 }
37 Point inter(Line p1,Line p2){
38     dou k1=(p2.e-p1.s)*(p1.e-p1.s);
39     dou k2=(p1.e-p1.s)*(p2.s-p1.s);
40     dou t=(k2)/(k1+k2);
41     dou x=p2.s.x+(p2.e.x-p2.s.x)*t;
42     dou y=p2.s.y+(p2.e.y-p2.s.y)*t;
43     return Point(x,y);
44 }
45 bool jud(Line p1,Line p2,Line p3){
46     Point p=inter(p1,p2);
47     return (p-p3.s)*(p3.e-p3.s)>0;
48 }
49 bool check(int mid){
50     int Tot=0;
51     for (int i=1;i<=tot;i++)
52      if (l[i].id<=mid) L[++Tot]=l[i];
53     int cnt=1;
54     for (int i=2;i<=Tot;i++)
55      if (L[i].slop!=L[i-1].slop) L[++cnt]=L[i]; 
56     int ll=1,rr=2;
57     c[1]=L[1];c[2]=L[2];
58     for (int i=3;i<=Tot;i++){
59         while (ll<rr&&jud(c[rr],c[rr-1],L[i])) rr--;
60         while (ll<rr&&jud(c[ll],c[ll+1],L[i])) ll++;
61         c[++rr]=L[i];
62     } 
63     while (ll<rr&&jud(c[rr],c[rr-1],c[ll])) rr--;
64     while (ll<rr&&jud(c[ll],c[ll+1],c[rr])) ll++;
65     if (rr-ll+1<3) return 0;
66     else return 1;
67 }
68 int main(){
69     n=read();
70     l[++tot].s=Point(-inf,inf);l[tot].e=Point(-inf,-inf);
71     l[++tot].s=Point(-inf,-inf);l[tot].e=Point(inf,-inf);
72     l[++tot].s=Point(inf,-inf);l[tot].e=Point(inf,inf);
73     l[++tot].s=Point(inf,inf);l[tot].e=Point(-inf,inf);
74     for (int i=1;i<=n;i++){
75         dou x=read(),y1=read(),y2=read();
76         l[++tot].s.x=-1;l[tot].s.y=y1/x-(-1)*x;
77         l[tot].e.x=1;l[tot].e.y=y1/x-x;
78         l[tot].id=i;
79         l[++tot].s.x=1;l[tot].s.y=y2/x-x;
80         l[tot].e.x=-1;l[tot].e.y=y2/x+x;
81         l[tot].id=i;
82     }
83     for (int i=1;i<=tot;i++) l[i].slop=atan2(l[i].e.y-l[i].s.y,l[i].e.x-l[i].s.x);
84     std::sort(l+1,l+1+tot,cmp);
85     int LL=1,RR=n,ans=0;
86     while (LL<=RR){
87         int mid=(LL+RR)>>1;
88         if (check(mid)) ans=mid,LL=mid+1;
89         else RR=mid-1;
90     }
91     printf("%d",ans);
92 }