#include<bits/stdc++.h>
#define pb push_back
#define P make_pair
#define fi first
#define se second
#define bit(s,x) (((s)>>(x))&1)
#define pnp(s) __builtin_popcountll(s)
#define rep(i,a,b) for(int i=(a);i<=(b);++i)
#define per(i,a,b) for(int i=(a);i>=(b);--i)
using namespace std;
typedef long long ll; 
typedef unsigned long long ull;
typedef pair<int,int> pii;
inline ll read(){
	ll x=0,f=1,c=getchar();
	while(c<'0'||c>'9')f=(c=='-'?-1:1),c=getchar();
	while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
const int N=2e3+3;
int n,col[N];
inline int ask(vector<int>vec,int i){
	vec.pb(i);
	printf("? %d ",vec.size());
	for(int x:vec)printf("%d ",x);
	printf("\n");fflush(stdout);
	int x=read(),y=read();
	if(x==i)return y;
	return x;
}
vector<int>G[N],vec[N];
mt19937 gen(114514);
int p[N];
int main() {
    //freopen("data.in", "r", stdin);
    //freopen("myans.out","w",stdout);
	n=read();
	rep(i,1,n)p[i]=i;
	shuffle(p+1,p+n+1,gen);
	int x=p[1];
	vec[col[x]=1].pb(x);
	rep(I,2,n){
		int i=p[I];
		rep(j,1,n){
			if(!vec[j].size())break;
			vector<int>tmp=vec[j];
			while(tmp.size()){
				int p=ask(tmp,i);
				if(p==-1)break;
				G[i].pb(p);
				rep(k,0,tmp.size()-1)if(tmp[k]==p){
					swap(tmp[k],tmp.back());tmp.pop_back();
					break;
				}
			}
		}
		static bool a[N];
		rep(j,1,n)a[j]=0;
		for(int p:G[i])a[col[p]]=1;
		for(col[i]=1;a[col[i]];++col[i]);
		vec[col[i]].pb(i);
	}
	int m=0;rep(i,1,n)m+=G[i].size();
	printf("! %d\n",m);
	rep(i,1,n)for(int j:G[i])printf("%d %d\n",i,j);
	fflush(stdout);
    return 0;
}