[bzoj3190][单调栈]赛车

Description

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

Input

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

Output

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

Sample Input

4

1 1 0 0

15 16 10 20

Sample Output

3

1 2 4

HINT

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

2016.1.17新加数据一组 By Nano_ape

题解

你先把这堆赛车的s-t图像搞出来，然后不就和水平可见直线差不多的吗？？
一开始我以为这堆射线是散落在第一象限里的，然后搞了搞发现原来是全部从y轴射出的。。
把射线按斜率排序，维护一个单调栈，下凸包
踢栈顶的方法有两个
1：斜率比栈顶不小，且出发点比栈顶高
2：这条直线与栈中第二个顶这条直线能完全遮蔽该直线，注意不能交于一点！
然后注意一下细节就过了？？？
那为什么网上都是半平面交

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
struct node
{
    double k,b;//y=kx+b
    int op;
}a[11000];int n,len;
bool cmp(node n1,node n2){return n1.k!=n2.k?n1.k<n2.k:n1.b<n2.b;}
double getpt(int u,int v){return (a[u].b-a[v].b)/(a[v].k-a[u].k);}
int sta[11000],tp;
int ans[11000];
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%lf",&a[i].b);
    for(int i=1;i<=n;i++)scanf("%lf",&a[i].k),a[i].op=i;
    sort(a+1,a+1+n,cmp);
    sta[++tp]=1;
    for(int i=2;i<=n;i++)
    {
        while(tp)
        {
            if(a[i].k==a[sta[tp]].k && a[i].b==a[sta[tp]].b)break;
            double fx=getpt(i,sta[tp]),gx=getpt(sta[tp],sta[tp-1]);
            if((a[i].k>=a[sta[tp]].k && a[i].b>a[sta[tp]].b) || fx<gx)tp--;
            else break;
        }
        sta[++tp]=i;
    }
    for(int i=1;i<=tp;i++)ans[i]=a[sta[i]].op;
    sort(ans+1,ans+1+tp);
    printf("%d\n",tp);
    for(int i=1;i<tp;i++)printf("%d ",ans[i]);
    printf("%d\n",ans[tp]);
    return 0;
}