BZOJ3190 [JLOI2013]赛车（单调栈+半平面交）

【题解】

数形结合思考：
画出v-t图像，若一条直线在一、四象限有不被覆盖的部分，它代表的车就可以领跑 
将直线按他们的斜率从小到大排序后，要维护一个下凸壳，由于一条新加入的直线，能覆盖的左边的直线是从右向左单调的，所以用单调栈来维护 

需要注意的细节：
1. 要在y轴右边做半平面交，我加入了一条过原点，斜率负无穷大的直线，覆盖掉了所有直线的左边一半 
2. 平行线无法算交点，重叠的直线都要算，要特殊判断 

【代码】

#include<stdio.h>
#include<stdlib.h>
#define eps 1e-13
double b[10005],k[10005];
int num[10005],ans[10005],sta[10005];
void jh(int* a,int* b)
{
	int t=*a;
	*a=*b;
	*b=t;
}
void jh(double* a,double* b)
{
	double t=*a;
	*a=*b;
	*b=t;
}
void kp(int low,int high)
{
	int i=low,j=high;
	double mk=k[(i+j)/2],mb=b[(i+j)/2];
	while(i<j)
	{
		while( k[i]<mk || (k[i]==mk&&b[i]<mb) ) i++;
		while( k[j]>mk || (k[j]==mk&&b[j]>mb) ) j--;
		if(i<=j)
		{
			jh(&k[i],&k[j]);
			jh(&b[i],&b[j]);
			jh(&num[i],&num[j]);
			i++;
			j--;
		}
	}
	if(j>low) kp(low,j);
	if(i<high) kp(i,high);
}
void kpans(int low,int high)
{
	int i=low,j=high,mid=ans[(i+j)/2];
	while(i<j)
	{
		while(ans[i]<mid) i++;
		while(ans[j]>mid) j--;
		if(i<=j)
		{
			jh(&ans[i],&ans[j]);
			i++;
			j--;
		}
	}
	if(j>low) kpans(low,j);
	if(i<high) kpans(i,high);
}
double getx(double k1,double b1,double k2,double b2)
{
	return (b2-b1)/(k1-k2);
}
int main()
{
	int n,i,j,top=0;
	scanf("%d",&n);
	for(i=1;i<=n;i++)
		scanf("%lf",&b[i]);
	for(i=1;i<=n;i++)
		scanf("%lf",&k[i]);
	for(i=1;i<=n;i++)
		num[i]=i;
	n++;//加入直线y=0
	k[n]=-1/eps;
	kp(1,n);
	for(i=1;i<=n;i++)
	{
		if(i>1&&k[i]==k[i-1])
		{
			if(b[i]==b[i-1])
			{
				sta[++top]=i;
				continue;
			}
			else
			{
				j=i-1;
				while(j>=1&&sta[top]==j) top--;
			}
		}
		while( top>1 && getx(k[i],b[i],k[sta[top]],b[sta[top]])-getx(k[sta[top-1]],b[sta[top-1]],k[sta[top]],b[sta[top]])<=-eps ) top--;
		sta[++top]=i;
	}
	for(i=1;i<=top;i++)
		ans[i]=num[sta[i]];
	kpans(1,top);
	printf("%d\n",top-1);
	for(i=2;i<=top;i++)
	{
		if(ans[i]!=0) printf("%d",ans[i]);
		if(i<top) printf(" ");
	}
	return 0;
}