#pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include constexpr long long INF = 1LL << 60; double PI = acos(-1.0); #define rep(i, n) for (ll i = 0; i < (n); ++i) #define rep1(i, n) for (ll i = 1; i <= (n); ++i) #define rrep(i, n) for (ll i = (n - 1); i >= 0; --i) #define perm(c) sort(ALL(c));for(bool c##p=1;c##p;c##p=next_permutation(ALL(c))) #define ALL(obj) (obj).begin(), (obj).end() #define RALL(obj) (obj).rbegin(), (obj).rend() #define pb push_back #define to_s to_string #define len(v) (ll)v.size() #define UNIQUE(v) v.erase(unique(v.begin(), v.end()), v.end()) #define print(x) cout << (x) << '\n' #define drop(x) cout << (x) << '\n', exit(0) #define debug(x) cout << #x << ": " << (x) << '\n' using namespace std; typedef long long ll; typedef unsigned long long ull; typedef long double ld; typedef pair P; typedef tuple tpl; typedef vector vec; typedef vector> vec2; typedef vector>> vec3; template inline bool chmax(S &a, const T &b) { if (a inline bool chmin(S &a, const T &b) { if (b ostream &operator << (ostream &os, const pair< T1, T2 > &p) { os << p.first << " " << p.second; return os; } template< typename T1, typename T2 > istream &operator >> (istream &is, pair< T1, T2 > &p) { is >> p.first >> p.second; return is; } template< typename T1, typename T2, typename T3 > ostream &operator << (ostream &os, const tuple< T1, T2, T3 > &t) { os << get<0>(t) << " " << get<1>(t) << " " << get<2>(t); return os; } template< typename T1, typename T2, typename T3 > istream &operator >> (istream &is, tuple< T1, T2, T3 > &t) { is >> get<0>(t) >> get<1>(t) >> get<2>(t); return is; } template< typename T > ostream &operator << (ostream &os, const vector< T > &v){ for (int i = 0; i < (int)v.size(); ++i) { os << v[i] << (i + 1 != v.size() ? " " : ""); } return os; } template< typename T > istream &operator >> (istream &is, vector< T > &v){ for(T &in : v) is >> in; return is; } template< typename T > ostream &operator << (ostream &os, const set< T > &st){ int ct = 0; for(auto& s : st) cout << s << (++ct != st.size() ? " " : ""); return os; } template constexpr set &operator|= (set &st1, const set &st2) { for(auto& s : st2) st1.insert(s); return st1; } template constexpr set &operator-= (set &st1, const set &st2) { for(auto& s : st2) if(st1.count(s)) st1.erase(s); return st1; } template constexpr set &operator&= (set &st1, const set &st2) { auto itr = st1.begin(); while(itr != st1.end()){ if(!st2.count(*itr)) itr = st1.erase(itr); else ++itr; } return st1; } template constexpr set operator| (const set &st1, const set &st2) { set res = st1; res |= st2; return res; } template constexpr set operator- (const set &st1, const set &st2) { set res = st1; res -= st2; return res; } template constexpr set operator& (const set &st1, const set &st2) { set res = st1; res &= st2; return res; } /*--------------------------------- Tools ------------------------------------------*/ template< typename T > vector cumsum(const vector &X){ vector res(X.size() + 1, 0); for(int i = 0; i < X.size(); ++i) res[i + 1] += res[i] + X[i]; return res; } template< typename S, typename T, typename F> pair bisearch(S left, T right, F f) { while(abs(right - left) > 1){ T mid = (right + left) / 2; if(f(mid)) right = mid; else left = mid; } return {left, right}; } template< typename S, typename T, typename F> double trisearch(S left, T right, F f, int maxLoop = 90){ double low = left, high = right; while(maxLoop--){ double mid_left = high / 3 + low * 2 / 3; double mid_right = high * 2 / 3 + low / 3; if(f(mid_left) >= f(mid_right)) low = mid_left; else high = mid_right; } return (low + high) * 0.5; } template< typename F > ll findMinimum(ll L, ll R, F f) { //[L, R) ll lo = L - 1, hi = R - 1; while (lo + 1 != hi) { ll mi = (lo + hi) / 2; if (f(mi) <= f(mi + 1)) hi = mi; else lo = mi; } return hi; } // #include // using namespace atcoder; constexpr long long MOD = 1000000007; /*------------------------------- Main Code Here -----------------------------------------*/ int main() { ll N, K; cin >> N >> K; vec A(N); cin >> A; auto f = [&](ll i) { vec ret; rep(j, N) ret.pb(abs(A[i] - A[j])); sort(ALL(ret)); //下からK個はpositive ll ans = 0; rep(j, N){ if(j < K) ans += ret[j]; else ans -= ret[j]; } return ans; }; ll bin = 30; ll step = devc(N, bin); ll ans = INF; rep(i, step){ ll ret = findMinimum(i * bin, min((i + 1) * bin, N), f); chmin(ans, f(ret)); } print(ans); return 0; }