BZOJ1007: [HNOI2008]水平可见直线

BZOJ1007: [HNOI2008]水平可见直线

https://lydsy.com/JudgeOnline/problem.php?id=1007

分析:

• 这玩意是个半平面交，但是他可以用凸包来对偶搞一搞。
• 把直线按斜率排序，然后搞一搞就做完了

代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <cmath>
using namespace std;
typedef double f2;
#define N 50050
int n,ans[N],S[N],tp;
struct A {
    f2 k,b;
    int id;
    bool operator < (const A &u) const {
        return abs(k-u.k)<1e-8 ? b<u.b : k<u.k;
    }
}a[N];
f2 X(A &x,A &y) {
    return -(x.b-y.b)/(x.k-y.k);
}
int main() {
    scanf("%d",&n);
    int i;
    for(i=1;i<=n;i++) {
        scanf("%lf%lf",&a[i].k,&a[i].b); a[i].id=i;
    }
    sort(a+1,a+n+1);

    for(i=1;i<=n;i++) {
        while(tp>1&&((abs(a[S[tp]].k-a[i].k)<1e-8)||X(a[S[tp-1]],a[S[tp]])>=X(a[S[tp]],a[i]))) tp--;
        if(tp==1&&(abs(a[S[tp]].k-a[i].k)<1e-8)) tp--;
        S[++tp]=i;
    }
    for(i=1;i<=tp;i++) ans[++ans[0]]=a[S[i]].id;
    sort(ans+1,ans+ans[0]+1);
    for(i=1;i<=ans[0];i++) printf("%d ",ans[i]);
}

转载于:https://www.cnblogs.com/suika/p/10092568.html