链接
一个模拟题,我们要从与上边界相切的圆开始考虑,然后考虑与这种圆相切的圆,如果这种圆和左边界有交点,就比较之后修改入点,和右边界有交点,就比较之后修改出点,因为他们上面都是出不去的
#include <bits/stdc++.h>
#define ll long long
#define pb push_back
#define inf 0x3f3f3f3f
#define rep(i,a,b) for(int i=a;i<b;i++)
#define rep1(i,a,b) for(int i=a;i>=b;i--)
#define rson rt<<1|1,m+1,r
#define lson rt<<1,l,m
using namespace std;
const int N=1000+100;
int vis[N],n;
double ans1,ans2;
struct node
{
int x,y,r;
}c[N];
bool can(node a,node b)//半段两个圆是否相交或相切
{
int dx = a.x - b.x;
int dy = a.y - b.y;
return dx*dx + dy*dy - (a.r + b.r)*(a.r + b.r) <= 0;
}
bool dfs(int u)
{
vis[u]=1;
if(c[u].y-c[u].r<=0) return 0;//与下边界有交点
if (c[u].x - c[u].r <= 0)ans1 = min(ans1, c[u].y - sqrt(c[u].r*c[u].r - c[u].x*c[u].x));//与左边界有交点,取纵坐标的较小者
if (1000 - c[u].x - c[u].r <= 0)ans2 = min(ans2, c[u].y - sqrt(c[u].r*c[u].r - (1000 - c[u].x)*(1000 - c[u].x)));//与右边界有交点,取纵坐标较小者
for (int i = 0; i < n; i++)
if (!vis[i]&&can(c[u], c[i])&&!dfs(i))//找到与它相切的圆继续判断
return false;
return true;
}
int main()
{
ios::sync_with_stdio(false);
while(cin>>n)
{
memset(vis,0,sizeof vis);
rep(i,0,n)
cin>>c[i].x>>c[i].y>>c[i].r;
ans1=ans2=1000;
int flag=1;
rep(i,0,n)
if(!vis[i]&&c[i].y+c[i].r>=1000&&!dfs(i))//从与上边界有交点的圆开始,如果没有则结果都是1000
{
cout<<"IMPOSSIBLE"<<endl;flag=0;
break;
}
if(flag) cout<<fixed<<setprecision(2)<<"0.00 "<<ans1<<" 1000.00 "<<ans2<<endl;
}
return 0;
}