#include #define rep(i, n) for (int i = 0; i < int(n); i++) #define rrep(i, n) for (int i = int(n) - 1; i >= 0; i--) #define reps(i, n) for (int i = 1; i <= int(n); i++) #define rreps(i, n) for (int i = int(n); i >= 1; i--) #define repc(i, n) for (int i = 0; i <= int(n); i++) #define rrepc(i, n) for (int i = int(n); i >= 0; i--) #define repi(i, a, b) for (int i = int(a); i < int(b); i++) #define repic(i, a, b) for (int i = int(a); i <= int(b); i++) #define each(x, y) for (auto &x : y) #define all(a) (a).begin(), (a).end() #define bit(b) (1ll << (b)) using namespace std; using i32 = int; using i64 = long long; using u64 = unsigned long long; using f80 = long double; using vi32 = vector; using vi64 = vector; using vu64 = vector; using vf80 = vector; using vstr = vector; inline void yes() { cout << "Yes" << '\n'; exit(0); } inline void no() { cout << "No" << '\n'; exit(0); } inline i64 gcd(i64 a, i64 b) { if (min(a, b) == 0) return max(a, b); if (a % b == 0) return b; return gcd(b, a % b); } inline i64 lcm(i64 a, i64 b) { return a / gcd(a, b) * b; } inline u64 xorshift() { static u64 x = 88172645463325252ull; x = x ^ (x << 7); return x = x ^ (x >> 9); } template class pqasc : public priority_queue, greater> {}; template class pqdesc : public priority_queue, less> {}; template inline void amax(T &x, T y) { if (x < y) x = y; } template inline void amin(T &x, T y) { if (x > y) x = y; } template inline T exp(T x, i64 n, T e = 1) { T r = e; while (n > 0) { if (n & 1) r *= x; x *= x; n >>= 1; } return r; } template inline T bis(T ok, T ng, function f, T eps = 1) { while (abs(ok - ng) > eps) { T mi = (ok + ng) / 2; (f(mi) ? ok : ng) = mi; } return ok; } template istream& operator>>(istream &is, vector &v) { each(x, v) is >> x; return is; } template ostream& operator<<(ostream &os, vector &v) { rep(i, v.size()) { if (i) os << ' '; os << v[i]; } return os; } template istream& operator>>(istream &is, pair &p) { is >> p.first >> p.second; return is; } template ostream& operator<<(ostream &os, pair &p) { os << p.first << ' ' << p.second; return os; } void solve(); int main() { ios::sync_with_stdio(0); cin.tie(0); cout << fixed << setprecision(16); solve(); return 0; } struct ConvexHullTrick { using F = pair; deque deq; #define a first #define b second #define lines deq.size() bool invalid(F f0, F f1, F f2) { return (__int128_t) (f1.a - f0.a) * (f2.b - f1.b) >= (__int128_t) (f1.b - f0.b) * (f2.a - f1.a); } i64 f(F f, i64 x) { return f.a * x + f.b; } void add(F f) { while (lines >= 2 && invalid(deq[lines - 2], deq[lines - 1], f)) { deq.pop_back(); } deq.push_back(f); } i64 min(i64 x) { while (lines >= 2 && f(deq[0], x) >= f(deq[1], x)) { deq.pop_front(); } return f(deq[0], x); } #undef a #undef b #undef lines }; i64 dp[300001]; void solve() { int n; cin >> n; vi64 A(n), x(n), y(n); cin >> A >> x >> y; rep(i, 300001) dp[i] = 1e18; dp[0] = 0; ConvexHullTrick cht; rep(i, n) { i64 a = -2 * x[i]; i64 b = x[i] * x[i] + y[i] * y[i] + dp[i]; cht.add(make_pair(a, b)); dp[i + 1] = cht.min(A[i]) + A[i] * A[i]; } cout << dp[n] << endl; }