#include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define int long long int MOD = 1000000007; template struct segment_tree { using M = monoid; using T = typename M::value_type; std::size_t sz; std::vector x; segment_tree(std::size_t n = 0) { sz = 1; while (sz < n) sz *= 2; x.assign(sz * 2, M::id()); initialize(); } template segment_tree(iterator first, iterator last) { sz = 1; std::size_t n = std::distance(first, last); while (sz < n) sz *= 2; x.assign(sz * 2, M::id()); std::copy(first, last, x.begin() + sz); initialize(); } void fill(const T &val) { std::fill(x.begin() + sz, x.end(), val); initialize(); } void initialize() { for (int i = (int)sz - 1; i >= 1; --i) { x[i] = M::op(x[i * 2 + 0], x[i * 2 + 1]); } } T accumulate(std::size_t l, std::size_t r) const { T al = M::id(), ar = M::id(); for (l += sz, r += sz; l < r; l /= 2, r /= 2) { if (l & 1) al = M::op(al, x[l++]); if (r & 1) ar = M::op(x[--r], ar); } return M::op(al, ar); } void update(std::size_t i, const T &val) { x[i += sz] = val; while (i > 1) { x[i / 2] = M::op(x[i], x[i ^ 1]); i /= 2; } } T operator[](std::size_t i) const { return x[sz + i]; } }; template struct min_monoid { using value_type = T; static constexpr T id() { return std::numeric_limits::max(); } static T op(const T &a, const T &b) { return std::min(a, b); } }; template struct max_monoid { using value_type = T; static constexpr value_type id() { return std::numeric_limits::min(); } static value_type op(const value_type &a, const value_type &b) { return std::max(a, b); } }; template struct sum_monoid { using value_type = T; static constexpr value_type id() { return 0; } static value_type op(const value_type &a, const value_type &b) { return a + b; } }; template using rmq = segment_tree>; template using rsq = segment_tree>; signed main() { cin.tie(0); ios::sync_with_stdio(false); int N; cin >> N; vector A(N); vector X(N); vector Y(N); vector dp(N + 1, 0); int res = 0; for (int i = 0; i < N; i++) { cin >> A[i]; A[i]++; } for (int i = 0; i < N; i++) { cin >> X[i]; X[i]++; } for (int i = 0; i < N; i++) { cin >> Y[i]; //Y[i]++; } rmq sg1(N); rmq sg2(N); dp[0] = 0; sg1.update(0, dp[0] - X[0] + Y[0]); sg2.update(0, dp[0] + X[0] + Y[0]); X.push_back(1 << 30); //sg1.update() for (int i = 1; i <= N; i++) { int t = lower_bound(X.begin(), X.end(), A[i - 1]) - X.begin(); int a; if (i <= t) { a = sg1.accumulate(0, i) + A[i - 1]; } else { if (t != 0) { a = sg1.accumulate(0, t) + A[i - 1]; } else { a = (int)1 << 60; } a = min(a, sg2.accumulate(t, i) - A[i - 1]); } dp[i] = a; if (i < N) { sg1.update(i, dp[i] - X[i] + Y[i]); sg2.update(i, dp[i] + X[i] + Y[i]); } } cout << dp[N] << endl; }