飛快秒掉了，沒報名痛失首殺，痛苦。
簡略題解：
考慮先隨機二元下標 \((x_0, y_0)\) 滿足刪去 \((x_0, y_0)\) 後查詢的中位數還是 \(\frac{n}{2}, \frac{n}{2} + 1\)，那麼這就說明 \(p_{x_0}, p_{y_0}\) 一定在中位數的兩邊。
那麼還剩下的 \(n - 2\) 個下標兩兩配對成 \(\frac{n - 2}{2}\) 組。
然後每一組與 \((x_0, y_0)\) 一起問，注意到只要二元組裡有中位數那麼問出來的就必定帶中位數（\(p_{x_0}, p_{y_0}\) 一定在中位數的兩邊，那麼加入中位數肯定被夾在中間，不管另一個在哪都肯定還在）。
於是規模就縮小到的 \(2\) 個二元組了，每個二元組各存在一箇中位數（可能只有 \(1\) 個二元組，但這個答案更簡單，就是這個二元組。）
那麼就只有 \(2\times 2 = 4\) 種情況，可以問出 \(3\) 種，都不是就是另一種（只是減掉了一次詢問，看起來不是很必要（？））。
於是後面部分的詢問次數就是 \(\frac{n - 3}{2} + 3 = \frac{n}{2} + 2\) 次，那麼就留下了 \(78 - \frac{n}{2}\) 次操作。
注意到隨機成功的機率就是 \(\frac{2(\frac{n}{2} - 1)^2}{n(n - 1)}\)，越小其實機率越低，但是越小的次數越多，算一下看起來就很能過.jpg。
跑了下 desmos，掛掉的機率 \(\le 8\times 10^{-9}\)，應該沒啥問題。
#include<bits/stdc++.h>
std::mt19937_64 eng(std::chrono::steady_clock::now().time_since_epoch().count());
constexpr int maxn = 1e2 + 10;
std::pair<int, int> query(std::vector<int> a) {
std::cout << "? " << a.size();
for (int x : a) std::cout << ' ' << x;
std::cout << std::endl;
int x, y;
std::cin >> x >> y;
return std::make_pair(x, y);
}
inline void solve() {
int n;
std::cin >> n;
int x0 = 0, y0 = 0;
do {
int x = eng() % n + 1, y;
do {
y = eng() % n + 1;
} while (y == x);
std::vector<int> vec;
for (int i = 1; i <= n; i++) {
if (i == x || i == y) continue;
vec.push_back(i);
}
if (query(vec) == std::make_pair(n / 2, n / 2 + 1)) {
x0 = x, y0 = y;
}
} while (! x0);
std::pair<int, int> c[2] = {{0, 0}, {0, 0}};
for (int i = 1, j = 0, k = 0; i <= n; i++) {
if (i == x0 || i == y0) continue;
if (j) {
auto [x, y] = query({i, j, x0, y0});
if (x == n / 2 || x == n / 2 + 1 || y == n / 2 || y == n / 2 + 1) {
c[k++] = {i, j};
}
j = 0;
} else {
j = i;
}
}
if (! c[1].first) {
std::cout << "! " << c[0].first << ' ' << c[0].second << std::endl;
} else if (query({c[0].first, c[1].first, x0, y0}) == std::make_pair(n / 2, n / 2 + 1)) {
std::cout << "! " << c[0].first << ' ' << c[1].first << std::endl;
} else if (query({c[0].first, c[1].second, x0, y0}) == std::make_pair(n / 2, n / 2 + 1)) {
std::cout << "! " << c[0].first << ' ' << c[1].second << std::endl;
} else if (query({c[0].second, c[1].first, x0, y0}) == std::make_pair(n / 2, n / 2 + 1)) {
std::cout << "! " << c[0].second << ' ' << c[1].first << std::endl;
} else {
std::cout << "! " << c[0].second << ' ' << c[1].second << std::endl;
}
}
int main() {
int T;
for (std::cin >> T; T--; ) {
solve();
}
return 0;
}