//#define _GLIBCXX_DEBUG #include using namespace std; #define endl '\n' #define lfs cout<= (ll)(n); i--) using ll = long long; using ld = long double; const ll MOD1 = 1e9+7; const ll MOD9 = 998244353; const ll INF = 1e18; using P = pair; template using PQ = priority_queue; template using QP = priority_queue,greater>; templatebool chmin(T1 &a,T2 b){if(a>b){a=b;return true;}else return false;} templatebool chmax(T1 &a,T2 b){if(avoid ans(bool x,T1 y,T2 z){if(x)cout<void anss(T1 x,T2 y,T3 z){ans(x!=y,x,z);}; templatevoid debug(const T &v,ll h,ll w,string sv=" "){for(ll i=0;ivoid debug(const T &v,ll n,string sv=" "){if(n!=0)cout<void debug(const vector&v){debug(v,v.size());} templatevoid debug(const vector>&v){for(auto &vv:v)debug(vv,vv.size());} templatevoid debug(stack st){while(!st.empty()){cout<void debug(queue st){while(!st.empty()){cout<void debug(deque st){while(!st.empty()){cout<void debug(PQ st){while(!st.empty()){cout<void debug(QP st){while(!st.empty()){cout<void debug(const set&v){for(auto z:v)cout<void debug(const multiset&v){for(auto z:v)cout<void debug(const array &a){for(auto z:a)cout<void debug(const map&v){for(auto z:v)cout<<"["<vector>vec(ll x, ll y, T w){vector>v(x,vector(y,w));return v;} ll gcd(ll x,ll y){ll r;while(y!=0&&(r=x%y)!=0){x=y;y=r;}return y==0?x:y;} vectordx={1,-1,0,0,1,1,-1,-1};vectordy={0,0,1,-1,1,-1,1,-1}; templatevector make_v(size_t a,T b){return vector(a,b);} templateauto make_v(size_t a,Ts... ts){return vector(a,make_v(ts...));} templateostream &operator<<(ostream &os, const pair&p){return os << p.first << " " << p.second;} templateostream &operator<<(ostream &os, const vector &v){for(auto &z:v)os << z << " ";cout<<"|"; return os;} templatevoid rearrange(vector&ord, vector&v){ auto tmp = v; for(int i=0;ivoid rearrange(vector&ord,Head&& head, Tail&&... tail){ rearrange(ord, head); rearrange(ord, tail...); } template vector ascend(const vector&v){ vectorord(v.size());iota(ord.begin(),ord.end(),0); sort(ord.begin(),ord.end(),[&](int i,int j){return v[i] vector descend(const vector&v){ vectorord(v.size());iota(ord.begin(),ord.end(),0); sort(ord.begin(),ord.end(),[&](int i,int j){return v[i]>v[j];}); return ord; } ll FLOOR(ll n,ll div){assert(div>0);return n>=0?n/div:(n-div+1)/div;} ll CEIL(ll n,ll div){assert(div>0);return n>=0?(n+div-1)/div:n/div;} ll digitsum(ll n){ll ret=0;while(n){ret+=n%10;n/=10;}return ret;} ll modulo(ll n,ll d){return (n%d+d)%d;}; templateT min(const vector&v){return *min_element(v.begin(),v.end());} templateT max(const vector&v){return *max_element(v.begin(),v.end());} templateT acc(const vector&v){return accumulate(v.begin(),v.end(),T(0));}; templateT reverse(const T &v){return T(v.rbegin(),v.rend());}; //mt19937 mt(chrono::steady_clock::now().time_since_epoch().count()); int popcount(ll x){return __builtin_popcountll(x);}; int poplow(ll x){return __builtin_ctzll(x);}; int pophigh(ll x){return 63 - __builtin_clzll(x);}; templateT poll(queue &q){auto ret=q.front();q.pop();return ret;}; templateT poll(priority_queue &q){auto ret=q.top();q.pop();return ret;}; templateT poll(QP &q){auto ret=q.top();q.pop();return ret;}; templateT poll(stack &s){auto ret=s.top();s.pop();return ret;}; ll MULT(ll x,ll y){if(LLONG_MAX/x<=y)return LLONG_MAX;return x*y;} ll POW2(ll x, ll k){ll ret=1,mul=x;while(k){if(mul==LLONG_MAX)return LLONG_MAX;if(k&1)ret=MULT(ret,mul);mul=MULT(mul,mul);k>>=1;}return ret;} ll POW(ll x, ll k){ll ret=1;for(int i=0;i struct edge { int to; T cost; int id; edge():id(-1){}; edge(int to, T cost = 1, int id = -1):to(to), cost(cost), id(id){} operator int() const { return to; } }; template using Graph = vector>>; template Graphrevgraph(const Graph &g){ Graphret(g.size()); for(int i=0;i Graph readGraph(int n,int m,int indexed=1,bool directed=false,bool weighted=false){ Graph ret(n); for(int es = 0; es < m; es++){ int u,v; T w=1; cin>>u>>v;u-=indexed,v-=indexed; if(weighted)cin>>w; ret[u].emplace_back(v,w,es); if(!directed)ret[v].emplace_back(u,w,es); } return ret; } template Graph readParent(int n,int indexed=1,bool directed=true){ Graphret(n); for(int i=1;i>p; p-=indexed; ret[p].emplace_back(i); if(!directed)ret[i].emplace_back(p); } return ret; } template vector>monotone_minima(int h, int w,const F &f){ vector>ret(h); auto dc=[&](auto &&self,int xl,int xr,int yl,int yr)->void{ if(xl >= xr)return; int mid = (xl + xr) / 2; ret[mid].first = -1; for(int y = yl; y < yr; y++){ T val = f(mid,y); if(ret[mid].first == -1 || ret[mid].second > val){ ret[mid] = make_pair(y,val); } } self(self,xl,mid,yl,ret[mid].first+1); self(self,mid+1,xr,ret[mid].first,yr); }; dc(dc,0,h,0,w); return ret; } template vector>smawk(int h, int w,const F &f){ vector>ret(h); auto rec=[&](auto &&self,const vector&row,const vector&col)->void{ if(row.empty())return; vectornr,nc; for(int i = 1; i < row.size(); i += 2){ nr.push_back(row[i]); } for(auto c:col){ while(!nc.empty() && f(row[nc.size()-1],nc.back()) > f(row[nc.size()-1],c)){ nc.pop_back(); } if(nc.size() < row.size()){ nc.push_back(c); } } self(self,nr,nc); int idx = 0; for(int i = 0; i < row.size(); i += 2){ int lim = nc.back(); if(i + 1 < row.size())lim = ret[row[i + 1]].first; ret[row[i]].first = -1; while(1){ T val = f(row[i], nc[idx]); if(ret[row[i]].first == -1 || ret[row[i]].second > val){ ret[row[i]] = make_pair(nc[idx], val); } if(nc[idx]==lim)break; idx++; } } }; vectorrow(h),col(w); iota(row.begin(),row.end(),0); iota(col.begin(),col.end(),0); rec(rec,row,col); return ret; } int main(){ cin.tie(nullptr); ios_base::sync_with_stdio(false); ll res=0,buf=0; bool judge = true; ll n;cin>>n; vectora(n),x(n),y(n); rep(i,0,n)cin>>a[i]; rep(i,0,n)cin>>x[i]; rep(i,0,n)cin>>y[i]; vectordp(n+1,INF); dp[0]=0; auto th=[&](ll k){ return k*k*k; }; auto rec=[&](auto &&f,ll l,ll r)->void{ if(r-l<=1)return; ll m=(l+r)/2; f(f,l,m); auto cost=[&](ll to,ll from){ from+=l,to+=m; return dp[from]+th(abs(a[to-1]-x[from]))+th(y[from]); }; auto minima=smawk(r-m,m-l,cost); //auto minima=monotone_minima(r-m,m-l,cost); rep(i,m,r)chmin(dp[i],minima[i-m].se); f(f,m,r); }; rec(rec,0,n+1); cout<