#include <bits/stdc++.h>

using i64 = long long;
using u64 = unsigned long long;
using u32 = unsigned;

void solve() {
    int n, m;
    std::cin >> n >> m;
    
    m = n - m;
    
    std::vector<int> a(n);
    for (int i = 0; i < n; i++) {
        std::cin >> a[i];
    }
    
    std::vector<int> ord(n);
    std::iota(ord.begin(), ord.end(), 0);
    std::sort(ord.begin(), ord.end(),
        [&](int i, int j) {
            if (a[i] != a[j]) {
                return a[i] < a[j];
            }
            return i < j;
        });
    
    if (m == 1) {
        std::cout << 1 << "\n";
        return;
    }
    
    i64 s1 = 0, s2 = 0;
    for (int i = n - m; i < n; i++) {
        s1 += a[ord[i]];
        s2 += a[ord[i]] * a[ord[i]];
    }
    i64 ans = s2 * m - s1 * s1;
     
    for (int i = 0; i < m; i++) {
        s1 += a[ord[i]];
        s2 += a[ord[i]] * a[ord[i]];
        s1 -= a[ord[i + n - m]];
        s2 -= a[ord[i + n - m]] * a[ord[i + n - m]];
        ans = std::max(ans, s2 * m - s1 * s1);
    }
    
    std::vector<int> p;
    
    s1 = 0, s2 = 0;
    for (int i = n - m; i < n; i++) {
        s1 += a[ord[i]];
        s2 += a[ord[i]] * a[ord[i]];
    }
    
    auto check = [&](int l) {
        std::vector<int> q;
        for (int i = 0; i < l; i++) {
            q.push_back(ord[i]);
        }
        int x = l + n - m;
        int y = x;
        while (y < n && a[ord[x]] == a[ord[y]]) {
            y++;
        }
        for (int i = y; i < n; i++) {
            q.push_back(ord[i]);
        }
        while (x > l && a[ord[x - 1]] == a[ord[y - 1]]) {
            x--;
            y--;
        }
        for (int i = x; i < y; i++) {
            q.push_back(ord[i]);
        }
        std::sort(q.begin(), q.end());
        if (p.empty() || p > q) {
            p = q;
        }
    };
    
    for (int i = 0; i <= m; i++) {
        if (ans == s2 * m - s1 * s1 && (i == 0 || a[ord[i - 1]] != a[ord[i - 1 + n - m]])) {
            check(i);
        }
        if (i < m) {
            s1 += a[ord[i]];
            s2 += a[ord[i]] * a[ord[i]];
            s1 -= a[ord[i + n - m]];
            s2 -= a[ord[i + n - m]] * a[ord[i + n - m]];
        }
    }
    for (int i = 0; i < m; i++) {
        std::cout << p[i] + 1 << " \n"[i == m - 1];
    }
}

int main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    
    int t;
    std::cin >> t;
    
    while (t--) {
        solve();
    }
    
    return 0;
}