1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66
| #include<bits/stdc++.h> using namespace std; int n; struct node { double x, y; int bian; }a[1010]; double f[1010][1010][2]; int pre[1010][1010][2]; int st; inline double dis(int x, int y) { return sqrt((a[x].x - a[y].x) * (a[x].x - a[y].x) + (a[x].y - a[y].y) * (a[x].y - a[y].y)); }
#define po(x) ((st + n - (x)) % n) #define pr(x) ((st + x) % n) int minn, minnc; void dfs(int x, int y, int c) { if(x || y) { if(!c) dfs(x - 1, y, pre[x][y][c]); else dfs(x, y - 1, pre[x][y][c]); } printf("%d ", (c) ? po(y) + 1 : pr(x) + 1); } int main() { scanf("%d", &n); for (int i = 0; i < n; ++i) { scanf("%lf %lf", &a[i].x, &a[i].y); a[i].bian = i; if (a[i].y > a[st].y) st = i; } for (int i = 0; i < n; ++i) { for (int j = 0; j < n; ++j) { f[i][j][0] = f[i][j][1] = 1e18; } } f[0][0][0] = f[0][0][1] = 0.0; for (int i = 1; i < n; ++i) { f[i][0][0] = f[i - 1][0][0] + dis(pr(i - 1), pr(i)); pre[i][0][0] = 0; } for (int i = 1; i < n; ++i) { f[0][i][1] = f[0][i - 1][1] + dis(po(i - 1), po(i)); pre[0][i][1] = 1; } for (int i = 1; i < n; ++i) { for (int j = 1; j < n - i; ++j) { double tmp1 = f[i - 1][j][0] + dis(pr(i - 1), pr(i)), tmp2 = f[i - 1][j][1] + dis(po(j), pr(i)); f[i][j][0] = min(tmp1, tmp2); if (f[i][j][0] == tmp1) pre[i][j][0] = 0; else pre[i][j][0] = 1; tmp1 = f[i][j - 1][1] + dis(po(j - 1), po(j)), tmp2 = f[i][j - 1][0] + dis(pr(i), po(j)); f[i][j][1] = min(tmp1, tmp2); if (f[i][j][1] == tmp1) pre[i][j][1] = 1; else pre[i][j][1] = 0; } } for (int i = 0; i < n; ++i) { if (f[i][n - i - 1][0] < f[minn][n - minn - 1][minnc]) minn = i, minnc = 0; if (f[i][n - i - 1][1] < f[minn][n - minn - 1][minnc]) minn = i, minnc = 1; } dfs(minn, n - minn - 1, minnc); return 0; }
|