先判断是否有解,从上到下dfs判断连通性,如果有从顶部到底部连通图,则无解。再判断最北的进出位置,从上边界开始遍历,沿途检查与边界相交的圆。这些圆的左边界的交点中最靠南边的一个就是所有的最北进入位置,和右边的最南交点就是所求的最北离开位置。
AC代码:
#include <iostream> #include <cstdio> #include <cstdlib> #include <cctype> #include <cstring> #include <string> #include <sstream> #include <vector> #include <set> #include <map> #include <algorithm> #include <stack> #include <queue> #include <bitset> #include <cassert> #include <cmath> using namespace std; const int maxn = 1005; const double W = 1000.0; int n, vis[maxn]; double x[maxn], y[maxn], r[maxn], leftP, rightP; bool flag; // 判断c1和c2是否相交 bool intersect(int c1, int c2) { return sqrt((x[c1] - x[c2]) * (x[c1] - x[c2]) + (y[c1] - y[c2]) * (y[c1] - y[c2])) < r[c1] + r[c2]; } void checkCircle(int u) { if (x[u] - r[u] < 0) { leftP = min(leftP, y[u] - sqrt(r[u] * r[u] - x[u] * x[u])); } if (x[u] + r[u] > W) { rightP = min(rightP, y[u] - sqrt(r[u] * r[u] - (W - x[u]) * (W - x[u]))); } } // 能不能到达底部 bool dfs(int u) { if (vis[u]) { return false; } vis[u] = 1; if (y[u] - r[u] < 0) { return true; } for (int v = 0; v < n; v++) { if (intersect(u, v) && dfs(v)) { return true; } } checkCircle(u); return false; } int main() { while (scanf("%d", &n) == 1) { flag = true; leftP = rightP = W; memset(vis, 0, sizeof(vis)); for (int i = 0; i < n; i++) { scanf("%lf%lf%lf", &x[i], &y[i], &r[i]); } for (int i = 0; i < n; i++) { if (y[i] + r[i] >= W && dfs(i)) { // 从上面开始dfs flag = false; break; } } if (flag) { printf("0.00 %.2lf %.2lf %.2lf ", leftP, W, rightP); } else { printf("IMPOSSIBLE "); } } return 0; }