#include using namespace std; template struct BinaryIndexedTree0 { const int N; vector val; BinaryIndexedTree0(int N) : N(N), val(N) {} T sum(int i) const { T s = 0; --i; if (i >= N) i = N - 1; for (; i >= 0; i = (i & (i + 1)) - 1) s = s + val[i]; return s; } T sum(int l, int r) const { return sum(r) - sum(l); } void add(int i, T x) { for (; i < N; i |= i + 1) val[i] = val[i] + x; } void set(int i, T x) { T c = sum(i, i + 1); add(i, x - c); } }; int64_t solve(const int n, const int k, const vector &A) { if (k == 1) return 0; // 座圧の前準備 vector compressor(A); sort(compressor.begin(), compressor.end()); compressor.erase(unique(compressor.begin(), compressor.end()), compressor.end()); auto compress = [&compressor](int i) { return lower_bound(compressor.begin(), compressor.end(), i) - compressor.begin(); }; // 座圧した結果出てくる数字の種類数 int m = compressor.size(); // median と、その座圧版を全部求めておく vector medians(n - k + 1); { multiset S(&A[0], &A[k]); auto it = next(S.begin(), k / 2); for (int i = k; i <= n; i++) { medians[i - k] = *it; S.insert(A[i]); if (A[i] < *it) it--; if (A[i-k] <= *it) it++; S.erase(S.lower_bound(A[i-k])); } } // bit で、今注目している長さ k の区間で、座圧後に i になる値の、個数と合計を管理する BinaryIndexedTree0 ft_cnt(m); BinaryIndexedTree0 ft_sum(m); int64_t ans = 1e18; for (int i = 0; i < n; i++) { // O(n) 回ループ // 値を入れる int a = A[i]; int c = compress(a); // O(log N) ft_cnt.add(c, 1); // O(log N) ft_sum.add(c, a); // O(log N) // 範囲外に出たものを取り除く O(log N) if (i >= k) { int r = A[i-k]; int c = compress(r); // O(log N) ft_cnt.add(c, -1); // O(log N) ft_sum.add(c, -r); // O(log N) } if (i >= k - 1) { int64_t median = medians[i-k+1]; int median_c = compress(median); // O(log N) // 区間内で median 未満の個数と総和を取得 int lower_cnt = ft_cnt.sum(median_c); // O(log N) int64_t lower_sum = ft_sum.sum(median_c); // O(log N) // 区間内で median 以上の個数と総和を取得 int higher_cnt = ft_cnt.sum(median_c, m); // O(log N) int64_t higher_sum = ft_sum.sum(median_c, m); // O(log N) int64_t inc_cost = median * lower_cnt - lower_sum; int64_t dec_cost = higher_sum - median * higher_cnt; int64_t cand = inc_cost + dec_cost; ans = min(ans, cand); } } return ans; } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); int n, k; cin >> n >> k; vector A(n); for (int i = 0; i < n; i++) cin >> A[i]; cout << solve(n, k, A) << endl; return 0; }