#include using namespace std; using V = pair; V operator+(const V &l, const V &r) { return {l.first + r.first, l.second + r.second}; } V operator-(const V &l, const V &r) { return {l.first - r.first, l.second - r.second}; } template struct BinaryIndexedTree0 { const int N; vector val; BinaryIndexedTree0(int N) : N(N), val(N) {} T sum(int i) const { T s = {0, 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 と、その座圧版を全部求めておく int num_medians = n - k + 1; vector medians(num_medians); { 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 の i 番目には、今注目している長さ k の区間で、座圧後に i になる値の、個数と合計を入れている。 // その range-sum ができるようになっている BinaryIndexedTree0 ft(m); int64_t ans = 1e18; for (int i = 0; i < n; i++) { // 値を入れる O(log N) int a = A[i]; ft.add(compress(a), {1, a}); // 範囲外に出たものを取り除く O(log N) if (i >= k) { int c = A[i-k]; ft.add(compress(c), {-1, -c}); } if (i >= k - 1) { int64_t median = medians[i-k+1]; int median_c = compress(median); // 区間内で median 未満の個数と総和を取得 V lower = ft.sum(median_c); // 区間内で median 以上の個数と総和を取得 V higher = ft.sum(median_c, m); int64_t inc_cost = median * lower.first - lower.second; int64_t dec_cost = higher.second - median * higher.first; 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; }