好不容易自己做出来个 (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;
}