#include<bits/stdc++.h>
using namespace std;
int n,m,k;
int mat[2005][2005];
vector<int> vec[2005],edge[2005];
int deg[2005],vis[2005];
priority_queue<pair<int,int>,vector<pair<int,int> >,greater<pair<int,int> > >que;
stack<int> st;
int num[2005];
int col[2005];
void rebuild(int n){
	for(int i=1;i<=n;i++)deg[i]=edge[i].size(),que.push(make_pair(deg[i],i));
	while(!que.empty()){
		int x=que.top().second;que.pop();
		if(vis[x])continue;
		k=max(k,deg[x]+1);
		vis[x]=1;
		st.push(x);
		for(int y:edge[x])if(!vis[y]){
			deg[y]--;
			que.push(make_pair(deg[y],y));
		}
	}
	for(int i=1;i<=k;i++)vec[i].clear();
	while(!st.empty()){
		int x=st.top();st.pop();
		vis[x]=0;
		for(int i=1;i<=k;i++)num[i]=0;
		for(int y:edge[x])if(!vis[y])num[col[y]]++;
		for(int i=1;i<=k;i++)if(!num[i]){col[x]=i;vec[i].push_back(x);break;}
	}
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		int ins=0;
		for(int j=1;j<=k;j++){
			int flag=1,cnt=0;
			while(cnt<vec[j].size()){
				putchar('?');
				int tot=1;
				for(int x:vec[j])if(!mat[x][i])tot++;
				printf(" %d",tot);
				for(int x:vec[j])if(!mat[x][i])printf(" %d",x);
				printf(" %d\n",i);
				fflush(stdout);
				int x,y;scanf("%d%d",&x,&y);
				if(x<0&&y<0)break;
				else{
					flag=0;
					mat[x^y^i][i]=1;
					edge[i].push_back(x^y^i);
					edge[x^y^i].push_back(i);
					cnt++;m++;
				}
			}
			if(flag&&!ins)ins=1,vec[j].push_back(i);
		}
		rebuild(i);
	}
	printf("! %d\n",m);
	for(int i=1;i<=n;i++)
		for(int j=i+1;j<=n;j++)if(mat[i][j])printf("%d %d\n",i,j);
	fflush(stdout);
	return 0;
}