#define PROBLEM "https://yukicoder.me/problems/no/705" #include template T monge_shortest_path(int n, Func cost_function) { assert(n > 0); std::vector dist(n); std::vector amin(n, 0); auto update = [&](int i, int k) { if (i <= k) return; T c = dist[k] + cost_function(k, i); if (c < dist[i]) { dist[i] = c; amin[i] = k; } }; auto rec_solve = [&](auto &&self, int l, int r) -> void { if (r - l == 1) return; int m = (l + r) / 2; for (int k = amin[l]; k <= amin[r]; k++) { update(m, k); } self(self, l, m); for (int k = l + 1; k <= m; k++) { update(r, k); } self(self, m, r); }; dist[0] = T(); for (int i = 1; i < n; i++) { dist[i] = cost_function(0, i); } rec_solve(rec_solve, 0, n - 1); return dist.back(); } void solve() { int n; std::cin >> n; std::vector A(n); std::vector X(n); std::vector Y(n); for (auto &a : A) std::cin >> a; for (auto &x : X) std::cin >> x; for (auto &y : Y) std::cin >> y; auto f = [&](int l, int r) { assert(l < r); long long dx = abs(A[r - 1] - X[l]); long long dy = Y[l]; return dx * dx * dx + dy * dy * dy; }; long long ans = monge_shortest_path(n + 1, f); std::cout << ans << std::endl; } int main() { std::cin.tie(0)->sync_with_stdio(0); int t = 1; // cin >> t; while (t--) solve(); } // #define PROBLEM "https://yukicoder.me/problems/no/705" // // #include // // #include "graph/monge_shortest_path.hpp" // // void solve() { // int n; // std::cin >> n; // std::vector A(n); // std::vector X(n); // std::vector Y(n); // for (auto &a : A) std::cin >> a; // for (auto &x : X) std::cin >> x; // for (auto &y : Y) std::cin >> y; // // auto f = [&](int l, int r) { // assert(l < r); // long long dx = abs(A[r - 1] - X[l]); // long long dy = Y[l]; // return dx * dx * dx + dy * dy * dy; // }; // // long long ans = monge_shortest_path(n + 1, f); // std::cout << ans << std::endl; // } // // int main() { // std::cin.tie(0)->sync_with_stdio(0); // // int t = 1; // // cin >> t; // while (t--) solve(); // }