#define _USE_MATH_DEFINES #include "bits/stdc++.h" using namespace std; #define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i)) #define rep(i,j) FOR(i,0,j) #define each(x,y) for(auto &(x):(y)) #define mp make_pair #define MT make_tuple #define all(x) (x).begin(),(x).end() #define debug(x) cout<<#x<<": "<<(x)<; using vi = vector; using vll = vector; template struct ConvexHullTrick { deque > lines; Compare comp; ConvexHullTrick() :comp(Compare()) { } void addLine(Val slope, Val constant) { pair line(slope, constant); int n; while ((n = (int)lines.size()) >= 2) { Val a1, a2, a3, b1, b2, b3; tie(a1, b1) = lines[n - 2]; tie(a2, b2) = lines[n - 1]; tie(a3, b3) = line; if ((a1 - a2)*(b3 - b1) >(a1 - a3)*(b2 - b1))break; else lines.pop_back(); } lines.push_back(line); } Val query(Val x) { while ((int)lines.size() >= 2 && comp(f(x, lines[1]), f(x, lines[0]))) { lines.pop_front(); } return f(x, lines[0]); } Val f(Val x, pair line) { return x * line.first + line.second; } }; const int MAX_N = 300005; int N, A[MAX_N], X[MAX_N], Y[MAX_N]; ll f[MAX_N]; void solve() { cin >> N; for (int *x : { A,X,Y }) { rep(i, N)cin >> x[i]; } ConvexHullTrick> cht; rep(i, N) { ll slope = -2ll * X[i]; ll constant = f[i] + (ll)X[i] * X[i] + (ll)Y[i] * Y[i]; cht.addLine(slope, constant); f[i + 1] = (ll)A[i] * A[i] + cht.query(A[i]); } cout << f[N] << endl; } int main() { ios::sync_with_stdio(false); cin.tie(0); cout << fixed << setprecision(15); solve(); return 0; }