cogs 2085 Asm.Def的一秒

“你们搞的这个导弹啊，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号可疑点的坐标。

【输出格式】

一行一个整数，即最多能经过几个可疑点。

如果夹角是$90^\circ$，直接按照x坐标排序，问题就变成了y坐标上的LIS。
本来坐标系也不非要是直角，按照他给的角度重新建系就行了。

#include<cstdio>
#include<cmath>
#include<vector>
#include<algorithm>
using namespace std;
const double eps=1e-8;
int dcmp(double x)
{
    if (x>eps) return 1;
    if (fabs(x)<eps) return 0;
    return -1;
}
int cmp1(double x,double y)
{
    return dcmp(x-y)==-1;
}
int cmp3(double x,double y)
{
    return dcmp(x-y)==0;
}
int cmp2(int x,int y)
{
    return x>y;
}
int rd()
{
    int x=0;
    char c=getchar();
    while (c<'0'||c>'9') c=getchar();
    while (c>='0'&&c<='9')
    {
        x=x*10+c-'0';
        c=getchar();
    }
    return x;
}
struct Vector
{
    double x,y;
    void read2()
    {
        y=rd();
        x=rd();
    }
    void read()
    {
        x=rd();
        y=rd();
    }
    Vector operator + (const Vector &v) const
    {
        return (Vector){x+v.x,y+v.y};
    }
    Vector operator - (const Vector &v) const
    {
        return (Vector){x-v.x,y-v.y};
    }
    Vector operator * (const int &k) const
    {
        return (Vector){x*k,y*k};
    }
    Vector operator / (const int &k) const
    {
        return (Vector){x/k,y/k};
    }
}a[100010],u,v,p;
double dot(Vector v1,Vector v2)
{
    return v1.x*v2.x+v1.y*v2.y;
}
double cross(Vector v1,Vector v2)
{
    return v1.x*v2.y-v1.y*v2.x;
}
int n,m,tot,mx[200010];
double ord[200010];
vector<int> f[200010];
vector<int>::iterator it;
int qry(int p)
{
    int ret=0;
    for (;p;p-=p&-p) ret=max(ret,mx[p]);
    return ret;
}
void modi(int p,int x)
{
    for (;p<=m;p+=p&-p) mx[p]=max(mx[p],x);
}
int main()
{
    freopen("asm_second.in","r",stdin);
    freopen("asm_second.out","w",stdout);
    double x1,y1;
    int x2,y2,ans=0;
    n=rd();
    u.read2();
    v.read2();
    for (int i=1;i<=n;i++)
    {
        p.read();
        x1=-cross(u,p)/cross(v,u);
        y1=-cross(v,p)/cross(u,v);
        if (dcmp(x1)==1&&dcmp(y1)==1) a[++tot]=(Vector){x1,y1};
    }
    for (int i=1;i<=tot;i++)
    {
        ord[++m]=a[i].x;
        ord[++m]=a[i].y;
    }
    sort(ord+1,ord+m+1,cmp1);
    m=unique(ord+1,ord+m+1,cmp3)-ord-1;
    for (int i=1;i<=tot;i++)
    {
        x2=lower_bound(ord+1,ord+m+1,a[i].x,cmp1)-ord;
        y2=lower_bound(ord+1,ord+m+1,a[i].y,cmp1)-ord;
        f[x2].push_back(y2);
    }
    for (int i=1;i<=m;i++)
    {
        sort(f[i].begin(),f[i].end(),cmp2);
        for (it=f[i].begin();it!=f[i].end();it++)
        {
            x2=qry(*it);
            ans=max(x2+1,ans);
            modi(*it+1,x2+1);
        }
    }
    printf("%d\n",ans);
}