#include <bits/stdc++.h>
using namespace std;

struct line {
	double a, b, c;
	int id;

	double y(double x) {
		return (c-a*x) / b;
	}

	double x(double y) {
		return (c-b*y) / a;
	}
};

int n;
long long h;
line a[100005];

template<class T, int size>
struct fenwick {
	T a[size];

	/* precondition: pos > 0 */
	void add(int pos, const T& val) {
		while (pos < size) {
			a[pos] += val;
			pos += pos & -pos;
		}
	}

	T sum(int pos) {
		T ret = T();
		while (pos > 0) {
			ret += a[pos];
			pos -= pos & -pos;
		}
		return ret;
	}

	void clear() {
		for (int i=0; i<size; i++) {
			a[i] = T();
		}
	}
};

fenwick<int, 100005> tree;

int main() {
	ios_base::sync_with_stdio(false);
	cin.tie(nullptr);
	cout.tie(nullptr);
	cerr.tie(nullptr);

	cin >> n;
	for (int i=1; i<=n; i++) {
		cin >> a[i].a >> a[i].b >> a[i].c;
	}

	h = (n*(n-1ll) / 2 + 1) / 2;

	cerr << "h: " << h << '\n';

	const double neginf = -1e30;

	double solx, soly;

	{
		sort(a+1, a+n+1, [&](line u, line v) { return u.y(neginf) < v.y(neginf); });
		for (int i=1; i<=n; i++) {
			a[i].id = i;
		}

		double l = -1e9, r = 1e9;

		for (int i=0; i<72; i++) {
			double mid = (l+r) / 2;

			sort(a+1, a+n+1, [&](line u, line v) { return u.y(mid) < v.y(mid); });

			long long invs = 0;

			tree.clear();

			for (int i=1; i<=n; i++) {
				tree.add(a[i].id, 1);
				invs += i - tree.sum(a[i].id);
			}

			if (invs < h) {
				l = mid;
			} else {
				r = mid;
			}
		}

		solx = (l+r) / 2;
	}

	{
		sort(a+1, a+n+1, [&](line u, line v) { return u.x(neginf) < v.x(neginf); });
		for (int i=1; i<=n; i++) {
			a[i].id = i;
		}

		double l = -1e9, r = 1e9;

		for (int i=0; i<72; i++) {
			double mid = (l+r) / 2;

			sort(a+1, a+n+1, [&](line u, line v) { return u.x(mid) < v.x(mid); });

			long long invs = 0;

			tree.clear();

			for (int i=1; i<=n; i++) {
				tree.add(a[i].id, 1);
				invs += i - tree.sum(a[i].id);
			}

			if (invs < h) {
				l = mid;
			} else {
				r = mid;
			}
		}

		soly = (l+r) / 2;;
	}

	cout.precision(20);
	cout << fixed << solx << ' ' << fixed << soly << '\n';		
	
}
### Prompt

Your task is to create a Cpp solution to the following problem:

There are N lines in the xy-plane. The i-th line is represented by A_ix+B_iy=C_i. Any two lines among the N+2 lines, the above N lines plus the x-axis and y-axis, cross each other at exactly one point.

For each pair 1 \leq i < j \leq N, there is a car at the cross point of the i-th and j-th lines. Even where three or more lines intersect at a point, a car is individually placed for each pair of lines. That is, there will be k(k-1)/2 cars placed at the intersection of k lines.

Those cars are already very old, and can only be moved parallel to the x-axis or y-axis.

Takahashi will hold an exhibition of antique cars at a place on the xy-plane. In order to avoid damaging the half-broken cars too much, he will select the place of the exhibition so that the total distance covered will be minimized when all the cars are moved to the place. If such a place is not uniquely determined, among the places that satisfy the condition above, the place with the minimum x-coordinate will be selected. If the place is still not uniquely determined, among the places that satisfy the two conditions above, the place with the minimum y-coordinate will be selected.

Find the place of the exhibition that will be selected.

Constraints

* 2 \leq N \leq 4 × 10^4
* 1 \leq |A_i|,|B_i| \leq 10^4(1 \leq i \leq N)
* 0 \leq |C_i| \leq 10^4(1 \leq i \leq N)
* No two given lines are parallel.
* All input values are integers.

Inputs

Input is given from Standard Input in the following format:


N
A_1 B_1 C_1
:
A_N B_N C_N


Outputs

Print the x-coordinate and y-coordinate of the place of the exhibition that will be selected, in this order, with a space in between. The output will be judged as correct when the absolute or relative error is at most 10^{-9}.

Examples

Input

3
1 1 1
2 -1 2
-1 2 2


Output

1.000000000000000 1.000000000000000


Input

4
1 1 2
1 -1 0
3 -1 -2
1 -3 4


Output

-1.000000000000000 -1.000000000000000


Input

7
1 7 8
-2 4 9
3 -8 -5
9 2 -14
6 7 5
-8 -9 3
3 8 10


Output

-1.722222222222222 1.325000000000000

### Response

```cpp
#include <bits/stdc++.h>
using namespace std;

struct line {
	double a, b, c;
	int id;

	double y(double x) {
		return (c-a*x) / b;
	}

	double x(double y) {
		return (c-b*y) / a;
	}
};

int n;
long long h;
line a[100005];

template<class T, int size>
struct fenwick {
	T a[size];

	/* precondition: pos > 0 */
	void add(int pos, const T& val) {
		while (pos < size) {
			a[pos] += val;
			pos += pos & -pos;
		}
	}

	T sum(int pos) {
		T ret = T();
		while (pos > 0) {
			ret += a[pos];
			pos -= pos & -pos;
		}
		return ret;
	}

	void clear() {
		for (int i=0; i<size; i++) {
			a[i] = T();
		}
	}
};

fenwick<int, 100005> tree;

int main() {
	ios_base::sync_with_stdio(false);
	cin.tie(nullptr);
	cout.tie(nullptr);
	cerr.tie(nullptr);

	cin >> n;
	for (int i=1; i<=n; i++) {
		cin >> a[i].a >> a[i].b >> a[i].c;
	}

	h = (n*(n-1ll) / 2 + 1) / 2;

	cerr << "h: " << h << '\n';

	const double neginf = -1e30;

	double solx, soly;

	{
		sort(a+1, a+n+1, [&](line u, line v) { return u.y(neginf) < v.y(neginf); });
		for (int i=1; i<=n; i++) {
			a[i].id = i;
		}

		double l = -1e9, r = 1e9;

		for (int i=0; i<72; i++) {
			double mid = (l+r) / 2;

			sort(a+1, a+n+1, [&](line u, line v) { return u.y(mid) < v.y(mid); });

			long long invs = 0;

			tree.clear();

			for (int i=1; i<=n; i++) {
				tree.add(a[i].id, 1);
				invs += i - tree.sum(a[i].id);
			}

			if (invs < h) {
				l = mid;
			} else {
				r = mid;
			}
		}

		solx = (l+r) / 2;
	}

	{
		sort(a+1, a+n+1, [&](line u, line v) { return u.x(neginf) < v.x(neginf); });
		for (int i=1; i<=n; i++) {
			a[i].id = i;
		}

		double l = -1e9, r = 1e9;

		for (int i=0; i<72; i++) {
			double mid = (l+r) / 2;

			sort(a+1, a+n+1, [&](line u, line v) { return u.x(mid) < v.x(mid); });

			long long invs = 0;

			tree.clear();

			for (int i=1; i<=n; i++) {
				tree.add(a[i].id, 1);
				invs += i - tree.sum(a[i].id);
			}

			if (invs < h) {
				l = mid;
			} else {
				r = mid;
			}
		}

		soly = (l+r) / 2;;
	}

	cout.precision(20);
	cout << fixed << solx << ' ' << fixed << soly << '\n';		
	
}
```