BZOJ1007: [HNOI2008]水平可见直线

计算几何相关

题目传送门

画画图就可以发现其实是维护一个下凸壳，像这样（图中绿色部分）：

其它的都不可能在答案里。

那么我们以斜率为第一关键字，截距为第二关键字排个序，斜率相同的保留截距最大的那一个。然后用一个单调栈维护。

当当前直线与栈顶直线交点在栈顶直线与栈顶下面那个直线交点的左边（就是比栈顶先交到，相当于把栈顶盖住）时弹出栈顶，如下图所示，绿色为当前直线，红色为栈顶直线。

那么最后栈中留下的就是所有可见的直线了。

因为不卡精度我就没写eps。

代码：

#include<cmath>
#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 50005
#define F inline
#define eps 1e-6
using namespace std;
struct L { int x,y,id; }a[N];
int n,tp,stk[N],ans[N];
F char readc(){
    static char buf[100000],*l=buf,*r=buf;
    if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
    if (l==r) return EOF; return *l++;
}
F int _read(){
    int x=0,f=1; char ch=readc();
    while (!isdigit(ch)) { if (ch=='-') f=-1; ch=readc(); }
    while (isdigit(ch)) x=(x<<3)+(x<<1)+(ch^48),ch=readc();
    return x*f;
}
F void writec(int x){
    if (x<0) putchar('-'),x=-x;
    if (x>9) writec(x/10); putchar(x%10+48);
}
F void _write(int x){ writec(x),putchar(' '); }
#define abs(x) ((x)>0?(x):-(x))
F bool cmp(L a,L b){ return a.x==b.x?a.y>b.y:a.x>b.x; }
F double calc(int x){//求交点
    return 1.0*(a[stk[tp]].y-a[x].y)/(a[x].x-a[stk[tp]].x);
}
int main(){
    n=_read();
    for (int i=1;i<=n;i++)
        a[i].x=_read(),a[i].y=_read(),a[i].id=i;
    sort(a+1,a+n+1,cmp);
    for (int i=1;i<=n;i++){
        if (i!=1&&a[i].x==a[i-1].x) continue;
        while (tp>1&&calc(i)>=calc(stk[tp-1])) tp--;
        stk[++tp]=i;
    }
    for (int i=1;i<=tp;i++) ans[i]=a[stk[i]].id;
    sort(ans+1,ans+tp+1);
    for (int i=1;i<=tp;i++) _write(ans[i]);
    return 0;
}