#include #define FOR(v, a, b) for(int v = (a); v < (b); ++v) #define FORE(v, a, b) for(int v = (a); v <= (b); ++v) #define REP(v, n) FOR(v, 0, n) #define REPE(v, n) FORE(v, 0, n) #define REV(v, a, b) for(int v = (a); v >= (b); --v) #define ALL(x) (x).begin(), (x).end() #define ITR(it, c) for(auto it = (c).begin(); it != (c).end(); ++it) #define RITR(it, c) for(auto it = (c).rbegin(); it != (c).rend(); ++it) #define EXIST(c,x) ((c).find(x) != (c).end()) #define LLI long long int #define fst first #define snd second #ifndef M_PI #define M_PI 3.14159265358979323846 #endif #ifndef M_E #define M_E 2.71828182845904523536 #endif #ifdef DEBUG #include #else #define dump(x) #endif #define gcd __gcd using namespace std; template constexpr T lcm(T m, T n){return m/gcd(m,n)*n;} template void join(ostream &ost, I s, I t, string d=" "){for(auto i=s; i!=t; ++i){if(i!=s)ost< istream& operator>>(istream &is, vector &v){for(auto &a : v) is >> a; return is;} template istream& operator>>(istream &is, pair &p){is >> p.first >> p.second; return is;} template T& chmin(T &a, const U &b){return a = (a<=b?a:b);} template T& chmax(T &a, const U &b){return a = (a>=b?a:b);} template void fill_array(T (&a)[N], const U &v){fill((U*)a, (U*)(a+N), v);} template class LiChaoSegmentTree{ using line = pair; struct node{ line l; node *left, *right; node(const line &l): l(l), left(nullptr), right(nullptr){} }; vector xs; int n; node *root; T e; function cmp; function chm; public: LiChaoSegmentTree(const vector &xs, const bool min_mode, const T &e): xs(xs), n(xs.size()), root(nullptr), e(e){ if(min_mode){ cmp = [](const T &a, const T &b){return a < b;}; chm = [](const T &a, const T &b){return min(a,b);}; }else{ cmp = [](const T &a, const T &b){return a > b;}; chm = [](const T &a, const T &b){return max(a,b);}; } } T apply(const line &l, const T &x) const{ return l.fst*x + l.snd; } protected: node* add(node *t, const int left, const int right, line l){ // [l,r] if(!t) return new node(l); if(cmp(apply(t->l, xs[left]), apply(l, xs[left])) and cmp(apply(t->l, xs[right]), apply(l, xs[right]))){ return t; } if(!cmp(apply(t->l, xs[left]), apply(l, xs[left])) and !cmp(apply(t->l, xs[right]), apply(l, xs[right]))){ t->l = l; return t; } int mid = (left + right) / 2; if(!cmp(apply(t->l, xs[mid]), apply(l, xs[mid]))) swap(t->l, l); if(!cmp(apply(t->l, xs[left]), apply(l, xs[left]))){ t->left = add(t->left, left, mid, l); }else{ t->right = add(t->right, mid+1, right, l); } return t; } public: void add_line(T a, T b){ root = add(root, 0, n-1, make_pair(a,b)); } protected: T query(node *t, const int left, const int right, const int i) const{ if(!t) return e; if(right == left) return apply(t->l, xs[i]); int mid = (left + right) / 2; if(i <= mid) return chm(apply(t->l, xs[i]), query(t->left, left, mid, i)); else return chm(apply(t->l, xs[i]), query(t->right, mid+1, right, i)); } public: T query(const T &x) const{ return query(root, 0, n-1, lower_bound(ALL(xs), x) - xs.begin()); } }; int main(){ cin.tie(0); ios::sync_with_stdio(false); int n; while(cin >> n){ vector a(n), x(n), y(n); cin >> a >> x >> y; LiChaoSegmentTree cht(a,true,LLONG_MAX); vector dp(n); dp[0] = y[0]*y[0] + (x[0]-a[0])*(x[0]-a[0]); cht.add_line(-2*x[1], dp[0]+y[1]*y[1]+x[1]*x[1]); FOR(i,1,n){ dp[i] = y[0]*y[0] + (a[i]-x[0])*(a[i]-x[0]); chmin(dp[i], cht.query(a[i]) + a[i]*a[i]); if(i