#ifndef LOCAL
#define FAST_IO
#endif

// ============
#include <algorithm>
#include <array>
#include <bitset>
#include <cassert>
#include <cmath>
#include <iomanip>
#include <iostream>
#include <list>
#include <map>
#include <numeric>
#include <queue>
#include <random>
#include <set>
#include <stack>
#include <string>
#include <tuple>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>

#define OVERRIDE(a, b, c, d, ...) d
#define REP2(i, n) for (i32 i = 0; i < (i32)(n); ++i)
#define REP3(i, m, n) for (i32 i = (i32)(m); i < (i32)(n); ++i)
#define REP(...) OVERRIDE(__VA_ARGS__, REP3, REP2)(__VA_ARGS__)
#define PER(i, n) for (i32 i = (i32)(n) - 1; i >= 0; --i)
#define ALL(x) begin(x), end(x)

using namespace std;

using u32 = unsigned int;
using u64 = unsigned long long;
using i32 = signed int;
using i64 = signed long long;
using f64 = double;
using f80 = long double;

template <typename T>
using Vec = vector<T>;

template <typename T>
bool chmin(T &x, const T &y) {
    if (x > y) {
        x = y;
        return true;
    }
    return false;
}
template <typename T>
bool chmax(T &x, const T &y) {
    if (x < y) {
        x = y;
        return true;
    }
    return false;
}

template <typename T>
Vec<tuple<i32, i32, T>> runlength(const Vec<T> &a) {
    if (a.empty()) {
        return Vec<tuple<i32, i32, T>>();
    }
    Vec<tuple<i32, i32, T>> ret;
    i32 prv = 0;
    REP(i, 1, a.size()) {
        if (a[i - 1] != a[i]) {
            ret.emplace_back(prv, i, a[i - 1]);
            prv = i;
        }
    }
    ret.emplace_back(prv, (i32)a.size(), a.back());
    return ret;
}

#ifdef INT128

using u128 = __uint128_t;
using i128 = __int128_t;

istream &operator>>(istream &is, i128 &x) {
    i64 v;
    is >> v;
    x = v;
    return is;
}
ostream &operator<<(ostream &os, i128 x) {
    os << (i64)x;
    return os;
}
istream &operator>>(istream &is, u128 &x) {
    u64 v;
    is >> v;
    x = v;
    return is;
}
ostream &operator<<(ostream &os, u128 x) {
    os << (u64)x;
    return os;
}

#endif

[[maybe_unused]] constexpr i32 INF = 1000000100;
[[maybe_unused]] constexpr i64 INF64 = 3000000000000000100;
struct SetUpIO {
    SetUpIO() {
#ifdef FAST_IO
        ios::sync_with_stdio(false);
        cin.tie(nullptr);
#endif
        cout << fixed << setprecision(15);
    }
} set_up_io;
// ============

#ifdef DEBUGF
#else
#define DBG(x) (void)0
#endif

// ============

#include <cassert>
#include <utility>
#include <functional>

template <typename F, typename I, typename Comp = std::less<>>
auto golden_search(const F &f, I mn, I mx, Comp comp = Comp()) -> std::pair<I, decltype(f(mn))> {
    using V = decltype(f(mn));
    assert(mn <= mx);
    I x0, x1, x2;
    {
        x0 = mn - 1;
        I s = 1, t = 1;
        while (x0 + s + t <= mx) {
            I u = s + t;
            s = t;
            t = u;
        }
        x1 = x0 + s + t;
        x2 = x0 + s;
    }
    V fx2 = f(x2);
    V fx3 = 0;
    while (2 * x2 != x0 + x1) {
        I x3 = x0 + x1 - x2;
        if (x3 > mx) {
            x1 = x0;
            x0 = x3;
        } else {
            fx3 = f(x3);
            if (comp(fx2, fx3)) {
                x1 = x0;
                x0 = x3;
            } else {
                x0 = x2;
                x2 = x3;
                fx2 = fx3;
            }
        }
    }
    return std::pair<I, V>(x2, fx2);
}
// ============
// ============

#include <algorithm>
#include <vector>

template <typename T>
class CoordinateCompression {
    std::vector<T> data;

    int size_sum() {
        return 0;
    }

    template <typename... Tail>
    int size_sum(const std::vector<T> &head, const Tail &...tail) {
        return (int) head.size() + size_sum(tail...);
    }

    void push() {}

    template <typename... Tail>
    void push(const std::vector<T> &head, const Tail &...tail) {
        for (const T &ele : head) {
            data.emplace_back(ele);
        }
        push(tail...);
    }

    void compress() {}

    template <typename... Tail>
    void compress(std::vector<T> &head, Tail &...tail) {
        for (T &ele : head) {
            ele = (T) (std::lower_bound(data.begin(), data.end(), ele) - data.begin());
        }
        compress(tail...);
    }

public:
    template <typename... V>
    CoordinateCompression(V &...v) {
        data.reserve(size_sum(v...));
        push(v...);
        std::sort(data.begin(), data.end());
        data.erase(std::unique(data.begin(), data.end()), data.end());
        compress(v...);
    }

    const T &operator[](const T &ele) const {
        return data[ele];
    }

    int size() const {
        return data.size();
    }
    
    bool contains(const T &ele) const {
        auto it = std::lower_bound(data.begin(), data.end(), ele);
        return it != data.end() && *it == ele;
    }
    
    T cc(const T &ele) const {
        return (T) (std::lower_bound(data.begin(), data.end(), ele) - data.begin());
    }
};
// ============

using usize = i32;

std::vector<std::vector<usize>> calc_inv(std::vector<usize> p) {
    usize n = p.size();
    
    std::vector<std::vector<usize>> sum(n + 1, std::vector<usize>(n + 1, 0));
    for (usize i = 0; i < n; ++i) {
        sum[i + 1] = sum[i];
        for (usize j = p[i]; j < n; ++j) {
            ++sum[i + 1][j + 1];
        }
    }
    
    std::vector<std::vector<usize>> inv(n, std::vector<usize>(n + 1, 0));
    for (usize i = 0; i < n; ++i) {
        usize s = 0;
        for (usize j = i; j < n; ++j) {
            s += sum[j][n] - sum[i][n] - sum[j][p[j] + 1] + sum[i][p[j] + 1];
            inv[i][j + 1] = s;
        }
    }
    return inv;
}

void solve() {
    i32 n, k;
    cin >> n >> k;
    Vec<i32> a(n);
    REP(i, n) {
        cin >> a[i];
    }
    
    if (k == n) {
        cout << 0 << '\n';
        return;
    }
    // divide into k segments
    ++k;
    
    CoordinateCompression<i32> cc(a);
    DBG(a);
    
    Vec<Vec<usize>> invs = calc_inv(a);
    DBG(invs);
    
    const auto relaxed = [&](i64 lambda) -> i64 {
        Vec<i64> dp(n + 1, INF64);
        dp[0] = 0;
        REP(i, 1, n + 1) {
            REP(j, i) {
                chmin(dp[i], dp[j] + invs[j][i] + lambda);
            }
        }
        return dp[n] - lambda * k;
    };
    i64 ans = golden_search(relaxed, (i64)(-n * n), (i64) (n * n), greater{}).second;
    cout << ans << '\n';
}

int main() {
    int t;
    std::cin >> t;
    for (int i = 0; i < t; ++i) {
        solve();
    }
}