Lightoj 1285

题目链接: 点这里
题意: n个点连成一个多边形,每个点只能连两条边。输出一种解。
思路:一开始我请教别人的思路是每次求凸包,在求剩下点的凸包,一直把所有的点求完,但是我每次连接外面的凸包和里面相邻的凸包的时候有错误。上网搜了一下题解,看到邝斌大神有极角排序做的,我有做了好长时间代码一直有错,后来找到一组反例,才发现邝斌大神的代码真是厉害。
邝斌大神的代码: 点这里

AC代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <stack>
#include <queue>
#include <vector>
#include <cmath>
#include <map>
#include <set>
#define ll long long
#define llu unsigned long long
using namespace std;
const int maxn = 100010;
const double eps = 1e-10;
const double PI  = acos(-1.0);
struct Point {
    double x,y;
    int w;
    Point (double x = 0,double y = 0) : x(x),y(y) {}
};
Point p[3050];
bool dcmp(double x) {
    if(fabs(x)-0.0 <= eps) return true;
    return false;
}
//
double dis(Point a,Point b) {
    return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y));
}

double cross(Point a,Point b,Point c) {
//    if(dcmp((b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y))) return 0;
    return ((b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y));
}
bool cmp1(Point a,Point b) {
    if(a.x != b.x) return a.x < b.x;
    else return a.y<b.y;
}
bool cmp2(Point a,Point b) {
    double m = cross(p[0],a,b);
    if(dcmp(m)) {
        return dis(p[0],a) < dis(p[0],b);
//        return dis(p[1],b) > dis(p[1],a) ? true : false;
    }else {
        if(m > 0) return true;
        else return false;
    }
}

int n;

int main(){
    int T; scanf("%d",&T);
    for(int cas=1; cas<=T; cas++) {
        scanf("%d",&n);
        for(int i=0;i<n;i++) {
            scanf("%lf%lf",&p[i].x,&p[i].y);
            p[i].w = i;
        }
        sort(p,p+n,cmp1);
        sort(p+1,p+n,cmp2);
        int f = 0;
        for(int i=n-2;i>=1;i--) {
            if(!dcmp(cross(p[0],p[n-1],p[i]))) {
                f=i; break;     
            }
        }
        printf("Case %d:\n",cas);
        if(!f) {
            printf("Impossible\n"); continue;
        }
        printf("%d",p[0].w);
        reverse(p+f+1,p+n);
        for(int i=1;i<n;i++) {
            printf(" %d",p[i].w);
        }
        printf("\n");
    }
    return 0;
}

/*
4
3
0 0 3 0 5 1
3
0 0 3 0 5 0
8
0 0 10 0 3 3 8 2 4 6 7 5 2 10 10 10
6
-1 -1 0 0 5 5 8 8 6 -3 9 9
9
0 0 1 1 2 2 3 3 4 4 1 0 2 0 3 0 4 0
*/
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值