JLOI2013 赛车

没想到吧这是一道半平面交（其实黄学长说过）

将问题转化到S-t图像上会发现问题其实是从y为正无穷处观察最顶端直线数及原始编号。

所以是可以半平面交的。

另外要注意的是，批注处作用为当tme结果小于零时能弹出该直线

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<cctype>
#include<iomanip>
#define eps 1e-6
using namespace std;
inline int read(){
int i=0,f=1;
char ch;
for(ch=getchar();ch<'0'||ch>'9';ch=getchar())
if(ch=='-') f=-1;
for(;ch>='0'&&ch<='9';ch=getchar())
i=(i<<3)+(i<<1)+(ch^48);
return i*f;
}
int n,topos,nop;
int ans[30300];
struct autoracer{
int pos,velo,id;
friend bool operator <(const autoracer &a,const autoracer &b){
if(a.velo!=b.velo) return a.velo<b.velo;
else return a.pos<b.pos;
}
friend double tme(const autoracer &a,const autoracer &b){
return (double) (a.pos-b.pos)/(b.velo-a.velo);
}
}c[30300],q[30300];
bool check(autoracer a,autoracer b,autoracer c){
return tme(a,b)>tme(a,c);
}
int main(){
n=read();
for(int i=1;i<=n;++i){
c[i].pos=read();
c[i].id=i;
}
for(int i=1;i<=n;++i)
c[i].velo=read();
sort(c+1,c+n+1);
topos=1;
for(int i=2;i<=n;++i){
if(c[i].velo!=c[i-1].velo||c[i].velo==c[i-1].velo&&c[i].pos==c[i-1].pos)
++topos;
c[topos]=c[i];
}
n=topos,nop=1;
q[nop]=c[1];ans[nop]=c[1].id;
for(int i=2;i<=n;++i){
while(nop>=1&&tme(q[nop],c[i])<-eps) --nop；//
while(nop>=2&&check(q[nop-1],q[nop],c[i])) --nop;
q[++nop]=c[i];
ans[nop]=c[i].id;
}
printf("%d\n",nop);
sort(ans+1,ans+nop+1);
for(int i=1;i<nop;++i)
printf("%d ",ans[i]);
printf("%d",ans[nop]);
return 0;
}