【题目描述】
“你们搞的这个导弹啊,excited!”
Asm.Def通过数据链发送了算出的疑似目标位置,几分钟后,成群结队的巡航导弹从“无蛤”号头顶掠过,布满了天空。
“一共发射了多少导弹?”
“十亿美元。”斯科特·华莱士回答,“单价100万,现在天上有1000多枚。这玩意能自动搜索10个可疑点,找到目标就发动攻击。”
“什么?10个?我给了它10万个点!”
“这会让它的程序崩溃的。好在你还有时间手动输入路径。”
“多久?”
“零……还有一秒,他们又给续上了一秒。”
“我想静静,别问我静静是谁。”
Asm.Def在第一象限内找到了n个可疑点。他需要为导弹规划路径。
如图所示,导弹一开始在(0,0)。它只能朝着一定的方向——即严格夹在图中两条射线间的方向(白色部分)前进。注意,它不能沿着这两条射线前进,当然也不能停在原地。
当导弹到达某个可疑点后,它仍然只能朝着该范围内的方向前进,如图。
Asm.Def想要让导弹经过尽可能多的可疑点。他需要在一秒钟内知道,最多能经过多少个可疑点。
【输入格式】
第1行1个整数n。
第2行4个整数a b c d:代表两条射线的斜率分别是a/b和c/d。保证0<=a,b,c,d<=10^5,a/b<c/d(即a/b是靠下的那条射线),a/b≠0/0,c/d≠0/0.
接下来n行,每行2个整数xi,yi(1<=xi,yi<=10^5),代表i号可疑点的坐标。
【输出格式】
一行一个整数,即最多能经过几个可疑点。
【样例输入】
15 1 3 2 1 3 1 6 2 4 2 2 5 4 5 6 6 3 4 1 6 2 1 7 4 9 3 5 3 1 3 15 5 12 4
【样例输出】
4
【提示】
这是最佳路径。注意,导弹不能前往位于射线上的点。
对于30%的数据,n<=1000,a=0,b=1,c=1,d=0。
对于60%的数据,n<=1000。
对于100%的数据,n<=10^5。
【来源】
在此键入。
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~计算几何+树状数组~
(事实证明COGS确实是要写文件操作的,不写会T掉啊……)
假设现在停留在点(x1,y1),如果可以向(x2,y2)转移,那么需要满足:
c/d<(y1-y2)/(x1-x2)<a/b,
化简(注意y1-y2<0,乘的时候要变号)得:x(点1)=y1*b-x1*a>y2*b-x2*a=x(点2),cd同理。
所以我们把每个点记为两个状态x,y,那么每个点可以向x、y均小于它的点转移。
所以我们先用unique和lower_bound求出每个点y性质的大小序号,然后用树状数组维护小于x性质的点最多走多少步,最后转移即可。
(要添加一个原点,容易计算~)
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long
int n,x1,y1,x2,y2,c[100001],tot,f[100001];
ll q[100001];
struct node{
ll x,y;
int id;
}a[100001];
ll read()
{
ll x=0,f=1;char ch=getchar();
while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}
while(ch>='0' && ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
void add(int u,int v)
{
for(;u<=n+1;u+=u&(-u)) c[u]=max(c[u],v);
}
int cal(int u)
{
int tot=0;
for(;u;u-=u&(-u)) tot=max(tot,c[u]);
return tot;
}
bool operator < (node u,node v)
{
return u.x==v.x ? u.y<v.y:u.x<v.x;
}
int main()
{
freopen("asm_second.in","r",stdin);
freopen("asm_second.out","w",stdout);
n=read();
y1=read();x1=read();y2=read();x2=read();
for(int i=1;i<=n;i++) a[i].x=read(),a[i].y=read(),a[i].id=i;
for(int i=0;i<=n;i++)
{
ll x=x2*a[i].y-y2*a[i].x,y=y1*a[i].x-x1*a[i].y;
a[i].x=x;a[i].y=y;q[++tot]=y;
}
sort(q+1,q+tot+1);
tot=unique(q+1,q+tot+1)-q;
for(int i=0;i<=n;i++) a[i].y=lower_bound(q+1,q+tot+1,a[i].y)-q;
sort(a,a+n+1);
for(int i=0,j=0;i<=n;i++)
{
while(a[j].x<a[i].x)
{
add(a[j].y,f[a[j].id]);j++;
}
f[a[i].id]=cal(a[i].y-1)+1;
}
printf("%d\n",f[0]-1);
return 0;
}