Paintball UVA - 11853（模拟，贪心）

链接
一个模拟题，我们要从与上边界相切的圆开始考虑，然后考虑与这种圆相切的圆，如果这种圆和左边界有交点，就比较之后修改入点，和右边界有交点，就比较之后修改出点，因为他们上面都是出不去的

#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;
}