題目鏈接
The Olympic Games in Bercouver are in full swing now. Here everyone has their own objectives: sportsmen compete for medals, and sport commentators compete for more convenient positions to give a running commentary. Today the main sport events take place at three round stadiums, and the commentator’s objective is to choose the best point of observation, that is to say the point from where all the three stadiums can be observed. As all the sport competitions are of the same importance, the stadiums should be observed at the same angle. If the number of points meeting the conditions is more than one, the point with the maximum angle of observation is prefered.

Would you, please, help the famous Berland commentator G. Berniev to find the best point of observation. It should be noted, that the stadiums do not hide each other, the commentator can easily see one stadium through the other.

Input

The input data consists of three lines, each of them describes the position of one stadium. The lines have the format x,??y,??r, where (x,?y) are the coordinates of the stadium’s center (?-??103?≤?x,??y?≤?103), and r (1?≤?r??≤?103) is its radius. All the numbers in the input data are integer, stadiums do not have common points, and their centers are not on the same line.

Output

Print the coordinates of the required point with five digits after the decimal point. If there is no answer meeting the conditions, the program shouldn’t print anything. The output data should be left blank.
Examples

Input

0 0 10
60 0 10
30 30 10

Output

30.00000 0.00000

思路

從中心開始移動，判斷誤差，不斷逼近答案。

AC

#include<bits/stdc++.h> using namespace std; struct point{double x, y, r; }circle[4]; double eps = 1e-5; double dis[4]; double find_dis(point a, point b) {return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y)); } double solve(point x) {double error = 0;for (int i = 0; i < 3; i++) {dis[i] = find_dis(x, circle[i]) / circle[i].r;}// 如果是答案，三個dis值應該相同。誤差分析：差的平方，加上平方更加精確 for (int i = 0; i < 3; i++) {error += (dis[i] - dis[(i + 1) % 3]) * (dis[i] - dis[(i + 1) % 3]); }return error; } int dx[] = {0, 0, 1, -1}; int dy[] = {1, -1, 0, 0}; int main() {// 假設起始點是重心 point ans;ans.x = 0; ans.y = 0; ans.r = 0; // freopen("in.txt", "r", stdin);for (int i = 0; i < 3; i++) {double x, y, r;scanf("%lf%lf%lf", &x, &y, &r);circle[i].x = x;circle[i].y = y;circle[i].r = r;ans.x += x / 3;ans.y += y / 3;}// 每次移動的步數 double move = 1.0;double error = solve(ans);while (move >= eps ) { point new_point;int flag = 1;// 枚舉四個方向 for (int i = 0; i < 4; i++) {new_point.x = ans.x + dx[i] * move;new_point.y = ans.y + dy[i] * move;double temp = solve(new_point);if (temp < error) {flag = 0;error = temp;ans.x = new_point.x;ans.y = new_point.y; }}if (flag) move /= 2; }//存在答案輸出 if (error < eps)printf("%.5lf %.5lf\n", ans.x, ans.y);return 0; }

總結

以上是生活随笔為你收集整理的C. Commentator problem的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。