数式で表すと、適当な平行移動と反転をして
　min. Σ| (x-x_i)^2 + (y-y_i)^2 - r^2 |
　s.t. x ≧ 0, y ≧ 0, r ≧ 0
を解く、という問題に落ちるんだけど、この目的関数は
凸でないし、滑らかでもないので非常に面倒な問題になるはず。

ただ、関数が綺麗だから逐次二次緩和くらいで十分解けそう。
とても実装する気にはならないけど。