好不容易自己做出来个 (2700) ...
首先我们对题目中的 (n_1n_2 + s_1s_2 - n_1s_2 - n_2s_1) 因式分解，得到 (F = (n_1 - s_1)(n_2 - s_2))
也就是说，(F eq 0) 当且仅当 (n_1 eq s_1) , (n_2 eq s_2) 。
然后，题目中说“至少有 (2) 个 (N / S) 的磁铁”，先考虑只有 (2) 个 (N / S) 的情况（一个子问题），我们必须得在 (n) 次操作内找出它们。
我们根据 (F eq 0) 的条件思考，要构造出“一个磁铁在 (L) 集合中，另一个磁铁在 (R) 集合中”的情况。在排除各种千奇百怪的方法后，主要有以下 3 种：
1.对于每个 (i) ，询问：(L[1,2,...,i-1], R[i+1,i+2,...,n])
2.对于每个 (i) ，询问：(L[i], R[1,2,...,i-1,i+1,...,n])
3.对于每个 (i) ，询问：(L[1,2,...,i-1], R[i])
然后推广到一般情况，1 可以扔掉，主要剩下后两种。
其实我一开始想的是 2，想了很长时间。中间还挺顺的，看似可以轻松出解，结果 (|num_S - num_N| = 1) 的情况实在做不了... （读者可以自己推一下试试）
然后灵光一现，想到了 3。
我们从左到右询问，第一次发现 (F eq 0) ，说明当前位置肯定是第二个 (N / S) 的位置。然后继续向后问，可以轻松地把后面的问出来。
咦，怎么才 (n - 1) 次操作，题目中的限制不是 (n + lfloor log n floor) 吗？是不是我错了？但我感觉很对啊！（当时想到这里，我还真的看了一下评测记录，看看大家是不是都带 (log)）
结果写程序的时候才发现，第一个 (N / S) 的位置还没求呢！当然，解决这个问题也很简单，二分一下就好了。这样询问次数就对上了。
这样就 A 了。
code(有点丑):
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define mit map<int,int>::iterator
#define sit set<int>::iterator
#define itrm(g,x) for(mit g=x.begin();g!=x.end();g++)
#define itrs(g,x) for(sit g=x.begin();g!=x.end();g++)
#define ltype int
#define rep(i,j,k) for(ltype(i)=(j);(i)<=(k);(i)++)
#define rap(i,j,k) for(ltype(i)=(j);(i)<(k);(i)++)
#define per(i,j,k) for(ltype(i)=(j);(i)>=(k);(i)--)
#define pii pair<int,int>
#define fi first
#define se second
#define mpr make_pair
#define pb push_back
#define fastio ios::sync_with_stdio(false)
const int inf=0x3f3f3f3f,mod=1000000007;
const double pi=3.1415926535897932,eps=1e-6;
int T,n;
vector<int> ok;
bool bad[2005];
int ask(int x,int l,int r){
printf("? %d %d
",1,r-l+1);
printf("%d
",x);
rep(i,l,r) printf("%d ",i);
puts("");
fflush(stdout);
int a;scanf("%d",&a);return a;
}
void solve()
{
rep(i,1,n) bad[i] = 0;
ok.clear();
scanf("%d",&n);
int lim = 0;
rep(i,2,n) {
lim = i;
int x = ask(i,1,i-1);
if(x) break;
}
ok.pb(lim);
int l = 1, r = lim - 1;
while(l < r) {
int mid = (l + r) >> 1;
if(ask(lim, 1, mid)) r = mid;
else l = mid + 1;
}
ok.pb(l);
rep(i,lim+1,n) if(ask(lim,i,i)) ok.pb(i);
rap(i,0,ok.size()) bad[ok[i]] = 1;
printf("! %d ",n - ok.size());
rep(i,1,n) if(!bad[i]) printf("%d ",i);
puts("");fflush(stdout);
}
int main()
{
scanf("%d",&T);
while(T--) solve();
return 0;
}