#ifdef __LOCAL #define _GLIBCXX_DEBUG #endif #include #include using namespace std; using namespace atcoder; typedef unsigned long long ull; typedef long long ll; const double PI = 3.14159265358979323846; #define rep(i, n) for (ll i = 0; i < (ll)(n); i++) // aよりもbが大きいならばaをbで更新する // (更新されたならばtrueを返す) template bool chmax(T& a, const T& b) { if (a < b) { a = b; // aをbで更新 return true; } return false; } // aよりもbが小さいならばaをbで更新する // (更新されたならばtrueを返す) template bool chmin(T& a, const T& b) { if (a > b) { a = b; // aをbで更新 return true; } return false; } // a^p // 2^3 = 2 * 2^2 // 2^2 = 2 * (2^1) // 2^1 = 2 ll powpow(ll a, ll p) { if (p == 0) return 1; if (p % 2 == 0) { ll half = powpow(a, p / 2); return half * half; } else { return a * powpow(a, p - 1); } } class combination { private: ll N; vector> C; public: explicit combination(ll _n) : N(_n) { C.resize(N + 1); rep(i, N + 1) C[i].resize(N + 1); C[0][0] = 1; rep(i, N) rep(j, i + 1) { C[i + 1][j + 1] += C[i][j]; C[i + 1][j] += C[i][j]; } } ll operator()(ll n, ll k) { if (n < k) return 0; if (n < 0 || k < 0) return 0; return C[n][k]; } }; // ref. https://drken1215.hatenablog.com/entry/2018/06/08/210000 class modcombination { private: ll mod; vector fac, finv, inv; public: modcombination(ll n, ll _mod) : mod(_mod) { fac.resize(n + 1); finv.resize(n + 1); inv.resize(n + 1); fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for (ll i = 2; i <= n; i++) { fac[i] = fac[i - 1] * i % mod; inv[i] = mod - inv[mod % i] * (mod / i) % mod; finv[i] = finv[i - 1] * inv[i] % mod; } } ll operator()(ll n, ll k) { if (n < k) return 0; if (n < 0 || k < 0) return 0; return fac[n] * (finv[k] * finv[n - k] % mod) % mod; } }; // cf. https://qiita.com/drken/items/a14e9af0ca2d857dad23 vector enum_divisors(ll n) { vector res; // sqrt(n)まで試し割り for (ll i = 1; i * i <= n; i++) { if (n % i == 0) { res.push_back(i); // 重複しないならばiの相方であるn/iも約数 // e.g. n=25のときのi=5は重複 if (n / i != i) res.push_back(n / i); } } sort(res.begin(), res.end()); return res; } // cf. https://qiita.com/drken/items/a14e9af0ca2d857dad23 map prime_factors(ll n) { map res; // sqrt(n)まで試し割り for (ll a = 2; a * a <= n; a++) { if (n % a != 0) continue; // nで割れる限り割る while (n % a == 0) { res[a]++; n /= a; } } if (n != 1) res[n]++; return res; } // p/q struct fraction { ll p, q; fraction(ll _p = 0, ll _q = 1) : p(_p), q(_q) { if (q == 0) { p = 1; return; } if (q < 0) { p = -p; q = -q; } ll g = gcd(p, q); p /= g; q /= g; } bool operator<(const fraction& other) const { return p * other.q < q * other.p; } bool operator<=(const fraction& other) const { return p * other.q <= q * other.p; } bool operator==(const fraction& other) const { return p == other.p && q == other.q; } }; // res[i][c] := i 文字目以降で最初に文字 c が登場する index (存在しないときは n) vector> calcNext(const string& S) { ll n = (ll)S.size(); vector> res(n + 1, vector(26, n)); for (ll i = n - 1; i >= 0; --i) { for (ll j = 0; j < 26; ++j) res[i][j] = res[i + 1][j]; res[i][S[i] - 'a'] = i; } return res; } // ref. https://algo-logic.info/bridge-lowlink/ struct LowLink { vector> G; vector ord, low; vector visited; vector> bridges; explicit LowLink(const vector>& _G) : G(_G) { visited.resize(G.size(), false); ord.resize(G.size(), 0); low.resize(G.size(), 0); ll k = 0; rep(i, (ll)G.size()) { if (visited[i]) continue; k = dfs(i, k); } } ll dfs(ll node, ll k, ll parent = -1) { visited[node] = true; ord[node] = k; low[node] = k; k++; for (auto g : G[node]) { if (!visited[g]) { k = dfs(g, k, node); low[node] = min(low[node], low[g]); if (ord[node] < low[g]) { bridges.emplace_back(node, g); } } else if (g != parent) { low[node] = min(low[node], ord[g]); } } return k; } }; // 3x3(アフィン変換) struct matrix { vector> a; matrix(const vector>& _a = {{1, 0, 0}, {0, 1, 0}, {0, 0, 1}}) : a(_a) {} matrix operator*(const matrix& x) { matrix res({{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}); rep(i, 3) rep(j, 3) rep(k, 3) { // cout << i << "," << j << " <- " << a[i][k] << "*" << x.a[k][j] << endl; res.a[i][j] += a[i][k] * x.a[k][j]; } return res; } vector operator*(const vector& x) { matrix other({{x[0], 0, 0}, {x[1], 0, 0}, {1, 0, 0}}); auto res = *this * other; return vector({res.a[0][0], res.a[1][0], res.a[2][0]}); } }; ll inverse_number(const vector& a) { ll n = a.size(); ll ret = 0; fenwick_tree t(n); rep(i, n) { ret += t.sum(a[i] + 1, n); t.add(a[i], 1); } return ret; } ll swap_distance(const vector& a, const vector& b) { if (a.size() != b.size()) { return -1; } ll n = a.size(); using P = pair; vector

_a(n), _b(n); rep(i, n) { _a[i] = {a[i], i}; _b[i] = {b[i], i}; } sort(_a.begin(), _a.end()); sort(_b.begin(), _b.end()); vector p(n); rep(i, n) { auto [va, ai] = _a[i]; auto [vb, bi] = _b[i]; if (va != vb) { return -1; } p[ai] = bi; } return inverse_number(p); } // {g,x,y}: ax+by = g tuple extgcd(ll a, ll b) { if (b == 0) { return {a, 1, 0}; } auto [g, x, y] = extgcd(b, a % b); return {g, y, x - (a / b) * y}; } int main() { cin.tie(nullptr); ios_base::sync_with_stdio(false); cout << fixed << setprecision(15); ll n; cin >> n; vector a(n), b(n); rep(i, n) cin >> a[i]; rep(i, n) cin >> b[i]; vector sa(n); rep(i, n) { sa[i] = abs(a[i]); } sort(sa.begin(), sa.end()); auto f = [&](ll x) -> ll { ll ret = 0; rep(i, n) { ret += b[i] * abs(x - a[i]); } return ret; }; ll left = -1e9; ll right = 1e9; while (abs(left - right) > 2) { ll m1 = (left * 2 + right) / 3; ll m2 = (left + right * 2) / 3; if (f(m1) > f(m2)) { left = m1; } else { right = m2; } } ll ans = f(left); ll j = left; for (ll i = left + 1; i <= right; i++) { if (ans > f(i)) { ans = f(i); j = i; } } cout << j << " " << ans << endl; return 0; }