洛谷 P3256 [JLOI2013]赛车【半平面交？？一次函数互相覆盖的判断】

【表示感觉跟那个水平可见直线一样。。。】
P3256 [JLOI2013]赛车
题目描述
这里有一辆赛车比赛正在进行，赛场上一共有N辆车，分别称为个g1，g2……gn。赛道是一条无限长的直线。最初，gi位于距离起跑线前进ki的位置。比赛开始后，车辆gi将会以vi单位每秒的恒定速度行驶。在这个比赛过程中，如果一辆赛车曾经处于领跑位置的话（即没有其他的赛车跑在他的前面），这辆赛车最后就可以得奖，而且比赛过程中不用担心相撞的问题。现在给出所有赛车的起始位置和速度，你的任务就是算出那些赛车将会得奖。

输入格式：
第一行有一个正整数N表示赛车的个数。接下来一行给出N个整数，按顺序给出N辆赛车的起始位置。再接下来一行给出N个整数，按顺序给出N辆赛车的恒定速度。

输出格式：
输出包括两行，第一行为获奖的赛车个数。第二行按从小到大的顺序输出获奖赛车的编号，编号之间用空格隔开，注意最后一个编号后面不要加空格。

输入样例#1：
4
1 1 0 0
15 16 10 20
输出样例#1：
3
1 2 4

对于100%的数据N<=10000, 0<=ki<=10^9, 0<=vi<=10^9

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=100005;
int n,d,z[N];
struct ss{
    int x,y,z;
}a[N];bool flag[N];
//Writer : SM
bool cmp(ss x,ss y){return x.x<y.x||x.x==y.x&&x.y<y.y;}
bool cmp1(int x,int y){return a[x].z<a[y].z;}
double X(int x,int y){
    double xx=double(a[x].y-a[y].y)/(a[y].x-a[x].x);
    return xx;
}bool pd(int x,int y,int z){
    double xx=X(x,y);
    double yy=X(x,z);
    return xx>yy;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        a[i].z=i;
        scanf("%d",&a[i].y);
    }for(int i=1;i<=n;i++)
        scanf("%d",&a[i].x);
    sort(a+1,a+n+1,cmp);
    for(int i=1;i<=n;i++){
        while(d>=1&&X(z[d],i)<0)d--;
        while(d>=2&&pd(z[d-1],z[d],i))d--;
        z[++d]=i;
    }sort(z+1,z+d+1,cmp1);
    printf("%d\n",d);
    for(int i=1;i<=d;i++)
        printf("%d ",a[z[i]].z);
}