数据分析
这两个OJ网站各有特点,其中BZOJ对一些特殊情况有预判(最重要的是抛物线的位置有限制),洛谷则对时间卡的很严,所以两个都有必要AC一下(我是洛谷第4个AC这题的,在这撒个花吧,heheda)
题目里说的很清楚了。
解题分析
这道题给的初始条件就是一条条线段,那么先试试转化一下
已知x0,y1,y2是常量,a和b是自变量,有没有发现,右边的式子可以看成平面上的一条直线呢?那再结合左边,就是判断是直线的上边还是下边,还有,直线把一个平面分成了两部分,每一部分就是一个——半平面,好吗,那就想到半平面交了。把每一关转换成两半平面,如果能穿过这K个靶子,那就说明这2K个半平面的交非空。然后核心思想就完了。
下面是温馨提示
1
有没有发现,如果加入的半平面越多,那么半平面交只会变小不会变大,那么这样一个递减的状态,可以考虑二分优化,而不是每次加入求一次,然后再加入再求一次。
2
求半平面交的时候要先对极角排序,但是每次二分拍一次很烦,所以可以一开始先对所有半平面标号,然后对极角排序,求半平面交的时候先扫一趟把符合范围内的加进去,然后再求半平面交。
3
还有一个莫名其妙的优化,根据前面写的推导公式(以上面的为例),可以得到对应有向直线点为(0,y1/x0),方向向量为(1,-x0),然而把方向向量改成(1/x0,-1),然后时间就优下来了。
4
洛谷的时间问题解决了,然后就是BZOJ,首先求半平面交的时候会给四边加个框,但是注意抛物线y=a* x^2+b*x中根据题意必须a<0,b>0,所以框要定在第二象限,还有,由于不带等号,所以不能为0,可以是一个很小的常数1e-15
5
所求的半平面交有可能是一个点或一条线段,这种时候队列中也til-hed<=1,这不是加个等号就能解决的,需要将靶子两端分别延长。这样点和线段就成了多边形,然后就解决了。
时间:O(n*logn); 空间:O(n);
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
struct fdata{
double x,y;
fdata (double x=0,double y=0):x(x),y(y){}
}p[200010];
struct hdata{
fdata p,v; double ang; int id;
hdata () {}
hdata (fdata p,fdata v,int id):p(p),v(v),id(id){ang=atan2(v.y,v.x);}
bool operator < (const hdata &b)const{
return ang<b.ang;
}
}a[200010],que[200010],lne[200010];
int n,tot,L,R,ans;
fdata operator + (const fdata a,const fdata b) {return fdata(a.x+b.x,a.y+b.y);}
fdata operator - (const fdata a,const fdata b) {return fdata(a.x-b.x,a.y-b.y);}
fdata operator * (const fdata a,const double x) {return fdata(a.x*x,a.y*x);}
double cross(const fdata a,const fdata b) {return a.x*b.y-a.y*b.x;}
inline void readi(int &x){
x=0; int f=1; char ch=getchar();
while ('0'>ch||ch>'9') {if (ch=='-') f=-f; ch=getchar();}
while ('0'<=ch&&ch<='9') {x=x*10+ch-'0'; ch=getchar();}
x*=f;
}
void _init()
{
freopen("archery.in","r",stdin);
freopen("archery.out","w",stdout);
readi(n); tot=4;
a[1]=hdata(fdata(0,1e-15),fdata(1,0),0);
a[2]=hdata(fdata(-1e-15,0),fdata(0,1),0);
a[3]=hdata(fdata(-1e-15,1e15),fdata(-1,0),0);
a[4]=hdata(fdata(-1e15,1e15),fdata(0,-1),0);
for (int i=1;i<=n;i++)
{
int x,y,z; readi(x); readi(y); readi(z);
double xx=x,yx=y-1e-15,yy=z+1e-15;
a[++tot]=hdata(fdata(0,yx/xx),fdata(1.0/xx,-1),i);
a[++tot]=hdata(fdata(0,yy/xx),fdata(-1.0/xx,1),i);
}
sort(a+1,a+tot+1);
}
bool pd_lft(fdata p,hdata a){return cross(a.v,p-a.p)>=0;}
fdata getJ(hdata a,hdata b){
fdata u=a.p-b.p;
double t=cross(b.v,u)/cross(a.v,b.v);
return a.v*t+a.p;
}
bool _check(int tem)
{
int k=0;
for (int i=1;i<=tot;i++)
if (a[i].id<=tem) lne[++k]=a[i];
int hed=1,til=1; que[1]=lne[1];
for (int i=2;i<=k;i++){
while (hed<til&&!pd_lft(p[til-1],lne[i])) til--;
while (hed<til&&!pd_lft(p[hed],lne[i])) hed++;
que[++til]=lne[i];
if (fabs(cross(que[til].v,que[til-1].v))<1e-15){
til--;
if (pd_lft(lne[i].p,que[til])) que[til]=lne[i];
}
if (hed<til) p[til-1]=getJ(que[til-1],que[til]);
}
while (hed<til&&!pd_lft(p[til-1],que[hed])) til--;
return til-hed>1;
}
void _solve()
{
L=1; R=n;
while (L<=R)
{
int mid=(R-L)/2+L;
if (_check(mid)) {ans=mid; L=mid+1;}
else R=mid-1;
}
printf("%d",ans);
}
int main()
{
_init();
_solve();
return 0;
}
PS:求管理员解释为什么上一篇博客莫名失踪了。555555……我还得重新写一次。