bzoj 3190: [JLOI2013]赛车 （半平面交）

3190: [JLOI2013]赛车

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 1721   Solved: 529
[ Submit][ Status][ Discuss]

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

Source

[ Submit][ Status][ Discuss]

题解：半平面交

其实这道题跟bzoj1007水平可见直线是基本上相同的。

有些细节需要注意，在计算的过程中可能会炸int。

如果好几条直线相交于一点，那么这些直线都要输出。

如果存在重合的直线，那么这些直线都要输出。

#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<algorithm>
#define N 10003
#define inf 1000000000
#define LL long long 
using namespace std;
struct data{
	double k,b; int id;
}a[N];
bool operator <(data a,data b) {
	return a.k>b.k||a.k==b.k&&a.b<b.b;
}
int n,cnt,mark[N],q[N],ans[N];
bool pd(data x1,data x2,data x3)
{
	LL r1=(LL)(x3.k-x2.k)*(LL)(x2.b-x1.b);
	LL r2=(LL)(x1.k-x2.k)*(LL)(x2.b-x3.b);
	return r1>r2;
}
int main()
{
	freopen("race12.in","r",stdin);
	freopen("my.out","w",stdout);
	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].id=i;
	a[0].k=-inf; a[0].b=0;
	sort(a,a+n+1);
	//for (int i=0;i<=n;i++) cout<<a[i].id<<" ";
	//cout<<endl;
	int head=0,tail=0;
	for (int i=0;i<=n;i++) {
		while (tail>0&&(pd(a[q[tail-1]],a[q[tail]],a[i])||a[i].k==a[q[tail]].k&&a[i].b>a[q[tail]].b)) 
		 mark[a[q[tail]].id]=0,tail--;
		if (a[i].k<a[q[tail]].k&&a[i].b<a[q[tail]].b) continue;
		q[++tail]=i; mark[a[q[tail]].id]=1;
	}
	//for (int i=1;i<=n;i++) cout<<mark[i]<<" ";
	//cout<<endl;
	cnt=0;
	for (int i=1;i<=n;i++) 
	 if (mark[i]) ans[++cnt]=i;
	printf("%d\n",cnt);
	for (int i=1;i<=cnt-1;i++) printf("%d ",ans[i]);
	printf("%d\n",ans[cnt]);
}