题目链接: 点这里
题意: 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
*/