/** * date : 2021-07-23 21:03:52 */ #define NDEBUG using namespace std; // intrinstic #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include // utility namespace Nyaan { using ll = long long; using i64 = long long; using u64 = unsigned long long; using i128 = __int128_t; using u128 = __uint128_t; template using V = vector; template using VV = vector>; using vi = vector; using vl = vector; using vd = V; using vs = V; using vvi = vector>; using vvl = vector>; template struct P : pair { template P(Args... args) : pair(args...) {} using pair::first; using pair::second; T &x() { return first; } const T &x() const { return first; } U &y() { return second; } const U &y() const { return second; } P &operator+=(const P &r) { first += r.first; second += r.second; return *this; } P &operator-=(const P &r) { first -= r.first; second -= r.second; return *this; } P &operator*=(const P &r) { first *= r.first; second *= r.second; return *this; } P operator+(const P &r) const { return P(*this) += r; } P operator-(const P &r) const { return P(*this) -= r; } P operator*(const P &r) const { return P(*this) *= r; } }; using pl = P; using pi = P; using vp = V; constexpr int inf = 1001001001; constexpr long long infLL = 4004004004004004004LL; template int sz(const T &t) { return t.size(); } template inline bool amin(T &x, U y) { return (y < x) ? (x = y, true) : false; } template inline bool amax(T &x, U y) { return (x < y) ? (x = y, true) : false; } template inline T Max(const vector &v) { return *max_element(begin(v), end(v)); } template inline T Min(const vector &v) { return *min_element(begin(v), end(v)); } template inline long long Sum(const vector &v) { return accumulate(begin(v), end(v), 0LL); } template int lb(const vector &v, const T &a) { return lower_bound(begin(v), end(v), a) - begin(v); } template int ub(const vector &v, const T &a) { return upper_bound(begin(v), end(v), a) - begin(v); } constexpr long long TEN(int n) { long long ret = 1, x = 10; for (; n; x *= x, n >>= 1) ret *= (n & 1 ? x : 1); return ret; } template pair mkp(const T &t, const U &u) { return make_pair(t, u); } template vector mkrui(const vector &v, bool rev = false) { vector ret(v.size() + 1); if (rev) { for (int i = int(v.size()) - 1; i >= 0; i--) ret[i] = v[i] + ret[i + 1]; } else { for (int i = 0; i < int(v.size()); i++) ret[i + 1] = ret[i] + v[i]; } return ret; }; template vector mkuni(const vector &v) { vector ret(v); sort(ret.begin(), ret.end()); ret.erase(unique(ret.begin(), ret.end()), ret.end()); return ret; } template vector mkord(int N, F f) { vector ord(N); iota(begin(ord), end(ord), 0); sort(begin(ord), end(ord), f); return ord; } template vector mkinv(vector &v) { int max_val = *max_element(begin(v), end(v)); vector inv(max_val + 1, -1); for (int i = 0; i < (int)v.size(); i++) inv[v[i]] = i; return inv; } } // namespace Nyaan // bit operation namespace Nyaan { __attribute__((target("popcnt"))) inline int popcnt(const u64 &a) { return _mm_popcnt_u64(a); } inline int lsb(const u64 &a) { return a ? __builtin_ctzll(a) : 64; } inline int ctz(const u64 &a) { return a ? __builtin_ctzll(a) : 64; } inline int msb(const u64 &a) { return a ? 63 - __builtin_clzll(a) : -1; } template inline int gbit(const T &a, int i) { return (a >> i) & 1; } template inline void sbit(T &a, int i, bool b) { if (gbit(a, i) != b) a ^= T(1) << i; } constexpr long long PW(int n) { return 1LL << n; } constexpr long long MSK(int n) { return (1LL << n) - 1; } } // namespace Nyaan // inout namespace Nyaan { template ostream &operator<<(ostream &os, const pair &p) { os << p.first << " " << p.second; return os; } template istream &operator>>(istream &is, pair &p) { is >> p.first >> p.second; return is; } template ostream &operator<<(ostream &os, const vector &v) { int s = (int)v.size(); for (int i = 0; i < s; i++) os << (i ? " " : "") << v[i]; return os; } template istream &operator>>(istream &is, vector &v) { for (auto &x : v) is >> x; return is; } void in() {} template void in(T &t, U &... u) { cin >> t; in(u...); } void out() { cout << "\n"; } template void out(const T &t, const U &... u) { cout << t; if (sizeof...(u)) cout << sep; out(u...); } void outr() {} template void outr(const T &t, const U &... u) { cout << t; outr(u...); } struct IoSetupNya { IoSetupNya() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(15); cerr << fixed << setprecision(7); } } iosetupnya; } // namespace Nyaan // debug namespace DebugImpl { template struct is_specialize : false_type {}; template struct is_specialize< U, typename conditional::type> : true_type {}; template struct is_specialize< U, typename conditional::type> : true_type {}; template struct is_specialize::value, void>> : true_type { }; void dump(const char& t) { cerr << t; } void dump(const string& t) { cerr << t; } void dump(const bool& t) { cerr << (t ? "true" : "false"); } template ::value, nullptr_t> = nullptr> void dump(const U& t) { cerr << t; } template void dump(const T& t, enable_if_t::value>* = nullptr) { string res; if (t == Nyaan::inf) res = "inf"; if constexpr (is_signed::value) { if (t == -Nyaan::inf) res = "-inf"; } if constexpr (sizeof(T) == 8) { if (t == Nyaan::infLL) res = "inf"; if constexpr (is_signed::value) { if (t == -Nyaan::infLL) res = "-inf"; } } if (res.empty()) res = to_string(t); cerr << res; } template void dump(const pair&); template void dump(const pair&); template void dump(const T& t, enable_if_t::value>* = nullptr) { cerr << "[ "; for (auto it = t.begin(); it != t.end();) { dump(*it); cerr << (++it == t.end() ? "" : ", "); } cerr << " ]"; } template void dump(const pair& t) { cerr << "( "; dump(t.first); cerr << ", "; dump(t.second); cerr << " )"; } template void dump(const pair& t) { cerr << "[ "; for (int i = 0; i < t.second; i++) { dump(t.first[i]); cerr << (i == t.second - 1 ? "" : ", "); } cerr << " ]"; } void trace() { cerr << endl; } template void trace(Head&& head, Tail&&... tail) { cerr << " "; dump(head); if (sizeof...(tail) != 0) cerr << ","; trace(forward(tail)...); } } // namespace DebugImpl #ifdef NyaanDebug #define trc(...) \ do { \ cerr << "## " << #__VA_ARGS__ << " = "; \ DebugImpl::trace(__VA_ARGS__); \ } while (0) #else #define trc(...) (void(0)) #endif // macro #define each(x, v) for (auto&& x : v) #define each2(x, y, v) for (auto&& [x, y] : v) #define all(v) (v).begin(), (v).end() #define rep(i, N) for (long long i = 0; i < (long long)(N); i++) #define repr(i, N) for (long long i = (long long)(N)-1; i >= 0; i--) #define rep1(i, N) for (long long i = 1; i <= (long long)(N); i++) #define repr1(i, N) for (long long i = (N); (long long)(i) > 0; i--) #define reg(i, a, b) for (long long i = (a); i < (b); i++) #define regr(i, a, b) for (long long i = (b)-1; i >= (a); i--) #define fi first #define se second #define ini(...) \ int __VA_ARGS__; \ in(__VA_ARGS__) #define inl(...) \ long long __VA_ARGS__; \ in(__VA_ARGS__) #define ins(...) \ string __VA_ARGS__; \ in(__VA_ARGS__) #define in2(s, t) \ for (int i = 0; i < (int)s.size(); i++) { \ in(s[i], t[i]); \ } #define in3(s, t, u) \ for (int i = 0; i < (int)s.size(); i++) { \ in(s[i], t[i], u[i]); \ } #define in4(s, t, u, v) \ for (int i = 0; i < (int)s.size(); i++) { \ in(s[i], t[i], u[i], v[i]); \ } #define die(...) \ do { \ Nyaan::out(__VA_ARGS__); \ return; \ } while (0) namespace Nyaan { void solve(); } int main() { Nyaan::solve(); } // struct SlopeTrick { private: using ll = long long; // x : 座標, c : 傾きの変化量 struct P { ll x, c; P(ll _x, ll _c) : x(_x), c(_c) {} friend bool operator<(const P& a, const P& b) { return a.x < b.x; } friend bool operator>(const P& a, const P& b) { return a.x > b.x; } }; ll addL, addR, min_y; priority_queue

L; priority_queue, greater<>> R; void pushL(ll x, ll c = 1) { L.emplace(x - addL, c); } void pushR(ll x, ll c = 1) { R.emplace(x - addR, c); } P getL() { return P{L.top().x + addL, L.top().c}; } P getR() { return P{R.top().x + addR, R.top().c}; } void popL() { L.pop(); } void popR() { R.pop(); } public: SlopeTrick() : addL(0), addR(0), min_y(0) {} void debug() { cerr << "addL, addR, min_y : "; cerr << addL << ", " << addR << ", " << min_y << endl; auto LL{L}; cerr << "L : "; while (!LL.empty()) { cerr << "( " << LL.top().x + addL << ", " << LL.top().c << " ), "; LL.pop(); } cerr << endl; auto RR{R}; cerr << "R : "; while (!RR.empty()) { cerr << "( " << RR.top().x + addR << ", " << RR.top().c << " ), "; RR.pop(); } cerr << endl; cerr << "min : "; cerr << "( " << get_min().first << ", " << get_min().second << " )" << endl; } // return {x, y} s.t. {argmin, min} pair get_min() { ll x = L.empty() ? R.empty() ? 0 : getR().x : getL().x; return {x, min_y}; } void shift_L(ll a) { addL += a; } void shift_R(ll a) { addR += a; } // f(x) <- f(x - a) void shift_x(ll a) { addL += a, addR += a; } // f(x) <- f(x) + a void shift_y(ll a) { min_y += a; } // add (x-a)_+ _____/ void add_xma(ll a, ll c = 1) { ll used = 0; while (used < c && !L.empty()) { auto [X, C] = getL(); if (X <= a) break; popL(); ll tmp = min(c - used, C); pushR(X, tmp); if (C != tmp) pushL(X, C - tmp); min_y += (X - a) * tmp; used += tmp; } if (used) pushL(a, used); if (c - used) pushR(a, c - used); } // add (a-x)_+ \_____ void add_amx(ll a, ll c = 1) { ll used = 0; while (used < c && !R.empty()) { auto [X, C] = getR(); if (X >= a) break; popR(); ll tmp = min(c - used, C); pushL(X, tmp); if (C != tmp) pushR(X, C - tmp); min_y += (a - X) * tmp; used += tmp; } if (used) pushR(a, used); if (c - used) pushL(a, c - used); } // add |x-a| \____/ void add_abs_xma(ll a, ll c = 1) { add_xma(a, c), add_amx(a, c); } // chmin right side \_/ -> \__ void chmin_right() { decltype(R){}.swap(R); } // chmin left side \_/ -> __/ void chmin_left() { decltype(L){}.swap(L); } // destructive merge void merge(SlopeTrick& r) { if (L.size() + R.size() < r.L.size() + r.R.size()) swap(*this, r); while (!r.L.empty()) { auto [x, c] = r.getL(); r.popL(); add_amx(x, c); } while (!r.R.empty()) { auto [x, c] = r.getR(); r.popR(); add_xma(x, c); } shift_y(r.min_y); } // eval f(x) O(|L| + |R|) ll eval(ll x) { ll y = min_y; auto LL{L}; while (!LL.empty()) { auto [X, C] = LL.top(); if ((X += addL) <= x) break; y += (X - x) * C; } auto RR{R}; while (!RR.empty()) { auto [X, C] = RR.top(); if ((X += addR) >= x) break; y += (x - X) * C; } return y; } }; using namespace Nyaan; // https://atcoder.jp/contests/abc127/tasks/abc127_f void ABC127F() { ini(Q); SlopeTrick st; rep(i, Q) { ini(c); if (c == 1) { inl(a, b); st.add_abs_xma(a); st.shift_y(b); } else { out(st.get_min()); } // st.debug(); } } // https://atcoder.jp/contests/dwango2016-prelims/tasks/dwango2016qual_e void dwacon2016qual_E() { ini(N, L); vvi ls(TEN(5) + 10); rep(i, N) { ini(t, l); ls[t].push_back(l); } SlopeTrick st; rep(i, sz(ls)) { if (ls.empty()) continue; st.chmin_right(); each(x, ls[i]) st.add_abs_xma(x); } out(st.get_min().second); } // https://atcoder.jp/contests/utpc2012/tasks/utpc2012_12 void utpc2012_L() { ini(N, C0); vl C(N); C[0] = C0; vvi g(N); rep1(i, N - 1) { ini(p, c); --p; g[p].push_back(i); C[i] = c; } auto dfs = [&](auto rc, int c) -> SlopeTrick { SlopeTrick st; each(d, g[c]) { SlopeTrick stc = rc(rc, d); stc.shift_x(1); st.merge(stc); } st.add_abs_xma(C[c]); st.chmin_right(); return move(st); }; out(dfs(dfs, 0).get_min().second); } void yuki1467() { ini(M, N); vl A(M), B(N); in(A, B); vl zip{A}; each(x, B) zip.push_back(x); zip = mkuni(zip); vl as(sz(zip)), bs(sz(zip)); each(x, A) as[lb(zip, x)]++; each(x, B) bs[lb(zip, x)]++; rep1(k, M) { SlopeTrick st; st.add_amx(0, TEN(10)); rep(i, sz(zip)) { st.shift_x(as[i] - bs[i] * k); st.chmin_right(); if (i != sz(zip) - 1) { st.add_abs_xma(0, zip[i + 1] - zip[i]); } else { st.add_xma(0, 1e10); } } out(st.get_min().second); } } void Nyaan::solve() { // ABC127F(); // dwacon2016qual_E(); // utpc2012_L(); yuki1467(); }