1007: [HNOI2008]水平可见直线[维护下凸壳]

#include<cstdio>
#include<algorithm>
using namespace std;
inline int read(){
    int x=0,f=1;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;
}
const int N=5e5+5;
struct node{int a,b,id;}c[N];
int n,cnt,pos[N];double slop[N];bool vis[N];
bool operator <(const node &x,const node &y){
    return x.a<y.a;
}
int main(){
    n=read();
    for(int i=1;i<=n;i++) c[i].a=read(),c[i].b=read(),c[i].id=i;
    sort(c+1,c+n+1);
    int now=1;
    for(int i=2;i<=n;i++){
        for(;c[i].a==c[now].a;i++){
            if(c[i].b>c[now].b){
                c[now]=c[i];
            }
        }
        if(i<=n) c[++now]=c[i];
    }
    pos[cnt=1]=1;slop[1]=-0x7fffffffffffffffLL;
    for(int i=2;i<=now;i++){
        double low;
        for(;;){//相对斜率关系(避免求交点) 
            low=(double)(c[pos[cnt]].b-c[i].b)/(double)(c[i].a-c[pos[cnt]].a);
            if(low<=slop[cnt])
                cnt--;
            else 
                break;
        }
        pos[++cnt]=i;
        slop[cnt]=low;
    }
    for(int i=1;i<=cnt;i++) vis[c[pos[i]].id]=1;
    for(int i=1;i<=n;i++) if(vis[i]) printf("%d ",i);
    return 0;
}