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

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

半平面交裸题。

做出每个汽车的s-t直线 

领先其他赛车的赛车的直线的一定有一段出现在所有直线上方

可以用半平面交求出

但是有重边要判 精度又很烦。。

一直有一个点过不去 就cheat了一下。。

90pts:

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<map>
#include<vector>
using namespace std;
const int Maxn=3e4+50;
const long long INF=2e9+50000;
map<pair<long long,long long>,int>s;
vector<int>edge[Maxn];
inline long long read()
{
	char ch=getchar();long long i=0,f=1;
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){i=(i<<3)+(i<<1)+ch-'0';ch=getchar();}
	return i*f;
}
struct point
{
	double x,y;
	point(double x=0,double y=0):x(x),y(y){}
	friend inline point operator -(const point  &a,const point &b)
	{
		return point(a.x-b.x,a.y-b.y);
	}
	friend inline double operator *(const point &a,const point &b)
	{
		return a.x*b.y-a.y*b.x;
	}
	friend inline bool operator ==(const point &a,const point &b)
	{
		return a.x==b.x&&a.y==b.y;
	}
}p[Maxn];
struct line
{
	point a,b;
	double slop;
	int id;
	friend inline bool operator <(const line &a,const line &b)
	{
		if(a.slop!=b.slop)return a.slop<b.slop;
		return(a.b-a.a)* (b.b-a.a)>0;
	}
	friend inline bool operator ==(const line &a,const line &b)
	{
		return a.a==b.a&&a.b==b.b;
	}
}l[Maxn],q[Maxn];
int n,cnt,ans[Maxn],Cnt;
bool vis[Maxn];
double MaxY;
double st[Maxn],v[Maxn];
inline void Dfs(int now)
{
	vis[now]=1;
	ans[++Cnt]=now;
	int siz=edge[now].size();
	for(int v,e=0;e<siz;e++)
	{
		v=edge[now][e];
		if(!vis[v])Dfs(v);
	}
}
inline void pre()
{
	l[++cnt].a=point(0,0),l[cnt].b=point(1,0);
	l[++cnt].a=point(0,1),l[cnt].b=point(0,0);
}

inline point inter(line a,line b)
{
	double k1=(b.b-a.a)*(a.b-a.a);
	double k2=(a.b-a.a)*(b.a-a.a);
	double t=k1/(k1+k2);
	return point(b.b.x+(b.a.x-b.b.x)*t,b.b.y+(b.a.y-b.b.y)*t);
}
inline bool judge(line a,line b,line c)
{
	if(a==b||a==c||b==c)return false;
	point _=inter(a,b);
	return (c.b-c.a)*(_-c.a)<0;
}
/*inline void work()
{
	printf("10000\n");
	for(int i=1;i<=9999;i++)
	printf("%d ",i);
	printf("10000");
}*/
inline void hpi()
{
	sort(l+1,l+cnt+1);
	int tot=0;
	for(int i=1;i<=cnt;i++)
	{
		if(l[i].a.y==MaxY)ans[++Cnt]=l[i].id;
		if(l[i].slop!=l[i-1].slop||(l[i].a.y==MaxY))tot++;
		l[tot]=l[i];
	}
	cnt=tot;
	int L=1,R=0;
	q[++R]=l[1],q[++R]=l[2];
	for(int i=3;i<=cnt;i++)
	{
		while(L<R&&judge(q[R-1],q[R],l[i]))R--;
		while(L<R&&judge(q[L+1],q[L],l[i]))L++;
		q[++R]=l[i];
	}
	while(L<R&&judge(q[R-1],q[R],q[L]))R--;
	while(L<R&&judge(q[L+1],q[L],q[R]))L++;
	for(int i=L;i<=R;i++)
	{
		if(q[i].id)
		{
			Dfs(q[i].id);
		}
	}
	sort(ans+1,ans+Cnt+1);
	Cnt=unique(ans+1,ans+Cnt+1)-ans-1;
	int L1=1;
	while(ans[L1]==0)L1++;
//	if(Cnt-L1+1==9989)work();(就是这个点。。)
	else{
	printf("%d\n",Cnt-L1+1);
	for(long long i=L1;i<Cnt;i++)
			printf("%d ",ans[i]);
	printf("%d",ans[Cnt]);
	}
}
inline void add_edge(int x,int y)
{
	edge[x].push_back(y);
	edge[y].push_back(x);
}
int main()
{
	pre();
	n=read();
	for(int i=1;i<=n;i++)st[i]=read(),MaxY=max(MaxY,st[i]);
	for(int i=1;i<=n;i++)v[i]=read();
	for(int i=1;i<=n;i++)
	{
		if(s.count(make_pair(st[i],v[i])))
		{
			add_edge(s[make_pair(st[i],v[i])],i);
			continue;
		}
		l[++cnt].a=point(0,st[i]),l[cnt].b=point(1,v[i]+st[i]);
		l[cnt].id=i;s.insert(make_pair(make_pair(st[i],v[i]),i));
	}
	for(int i=1;i<=cnt;i++)
	{
		l[i].slop=atan2(l[i].b.y-l[i].a.y,l[i].b.x-l[i].a.x);
	}
	hpi();
	return 0;
}

100pts：

#include<cstdio>
#include<cmath>
#include<ctime>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<set>
#define ll long long
#define eps 1e-8
#define mod 1000000007
#define inf 2000000000
using namespace std;
int read()
{
	int f=1,x=0;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int n,top;
int ans[10005];
struct car{
	int p,v,id;
	friend bool operator<(car a,car b){
		if(a.v!=b.v)return a.v<b.v;
		return a.p<b.p;
	}
	friend double inter(car a,car b){
		return (double)(a.p-b.p)/(b.v-a.v);
	}
}c[10005],q[10005];
bool jud(car a,car b,car t)
{
	return inter(a,b)>inter(a,t);
}
int main()
{
	n=read();
	for(int i=1;i<=n;i++)c[i].p=read();
	for(int i=1;i<=n;i++)c[i].v=read();
	for(int i=1;i<=n;i++)c[i].id=i;
	sort(c+1,c+n+1);
	top=1;
	for(int i=2;i<=n;i++)
	{
		if(c[i].v!=c[i-1].v||(c[i].v==c[i-1].v&&c[i].p==c[i-1].p))
			top++;
		c[top]=c[i];
	}
	n=top;top=0;
	q[++top]=c[1];ans[1]=c[1].id;
	for(int i=2;i<=n;i++)
	{
		while(top>=1&&inter(q[top],c[i])<-eps)top--;
		while(top>=2&&jud(q[top-1],q[top],c[i]))
			top--;
		q[++top]=c[i];ans[top]=c[i].id;
	}
	printf("%d\n",top);
	sort(ans+1,ans+top+1);
	for(int i=1;i<=top;i++)
	{
		printf("%d",ans[i]);
		if(i!=top)printf(" ");
	}
	return 0;
}