#include<bits/stdc++.h>
using namespace std;
//#define int long long
#define double long double
#define mp(x,y) make_pair(x,y)
#define lb(x) (x&(-x))
#define fi first
#define se second
#define pii pair<int,int>
#define piii pair<int,pair<int,int> >
int read(){
	int x=0,f=1;
    char ch=getchar();
	while(ch<'0'||ch>'9'){
    	if(ch=='-')f=-1;
		ch=getchar();
    }
	while(ch>='0'&&ch<='9'){
    	x=x*10+ch-48;
		ch=getchar();
    }
	return x*f;
}
void writ(int x){
    if(x<0){
		putchar('-'),x=-x;
	}
    if(x>9){
		writ(x/10);
	}
    putchar(x%10 +'0');
}
void write(int x,char p=' '){
    writ(x);
    putchar(p);
}
string sread(){
    string t="";
    char c=getchar();
    while(c==' '||c=='\t'||c=='\r'||c=='\n'||c==0||c==EOF){
        c=getchar();
    }
    while(!(c==' '||c=='\t'||c=='\r'||c=='\n'||c==0||c==EOF)){
        t+=c,c=getchar();
    }
    return t;
}
const int QMR=1000000007,LSJ=998244353,inf=2e9+31,N=4005;
//int cal(int a,int b=QMR-2){
//	int base=a,ret=1;
//	while(b){
//		if(b&1){
//			ret=ret*base%QMR;
//		}
//		base=base*base%QMR;
//		b>>=1;
//	}
//	return ret;
//}
mt19937 rnd(time(0)^clock());
int n,a[N];
//map<basic_string<short>,basic_string<short> >h;
//void query(basic_string<short>t){
////	if(h.find(t)!=h.end())return h[t];
//	cout<<"? ";
//	write(t.size());
//	for(int i:t)write(i);
//	cout<<endl;
//}
int cnt=0;
basic_string<short>ans[N],b[N],qmr;
void work(int k,int l,int r,basic_string<short>s,basic_string<short>c){
	if(s.size()==0)return;;
	if(l==r){
		s=b[k][l]+s;
		for(int i=1;i<s.size();i++){
			for(int j=1;j<=c[i-1];j++){
				ans[s[i]]+=b[k][l];
			}
		}
		return;
	}
	int mid=(l+r)>>1;
	basic_string<short>L,R,ll,rr,Cl,Cr/*,CrCl3*/;
	for(int i=l;i<=mid;i++){
		L+=b[k][i];
	}
	for(int i=mid+1;i<=r;i++){
		R+=b[k][i];
	}
	cout<<"? "<<L.size()+s.size()<<" ";
	for(int i:L)cout<<i<<" ";
	for(int i:s)cout<<i<<" ";
	cout<<endl;
	basic_string<short>lsj;
	for(int i=0;i<L.size()+s.size();i++){
		int t=read();lsj+=t;
	}
	for(int i=L.size();i<lsj.size();i++){
		if(lsj[i]-ans[s[i-L.size()]].size()>0){
			ll+=s[i-L.size()];
			Cl+=lsj[i]-ans[s[i-L.size()]].size();
		}
		if(c[i-L.size()]-(lsj[i]-ans[s[i-L.size()]].size())>0){
			rr+=s[i-L.size()];
			Cr+=c[i-L.size()]-(lsj[i]-ans[s[i-L.size()]].size());
		}
	}
	work(k,l,mid,ll,Cl);
	ll.clear();Cl.clear();
	work(k,mid+1,r,rr,Cr);rr.clear();Cr.clear();
}
signed main(){
    //freopen("qmr.in","r",stdin);
    //freopen("qmr.out","w",stdout);
    ios::sync_with_stdio(0);cin.tie(0);
    n=read();
    basic_string<short>g;
    for(int i=1;i<=n;i++)g+=i;
    while(g.size()){
    	cnt++;
		cout<<"? "<<g.size()<<" ";
		for(int i:g)cout<<i<<" ";
		cout<<endl;
    	basic_string<short>tt,c,d;
		for(int i=0;i<g.size();i++){
			int t=read();tt+=t;
		}
    	for(int i=0;i<tt.size();i++){
    		if(tt[i]!=0)c+=g[i];
    		else b[cnt]+=g[i];
		}
		g=c;
	}
	basic_string<short>tmp=b[cnt];
	for(int i=cnt-1;i>=1;i--){
//		query(b[i]+b[i]+tmp);
		cout<<"? "<<b[i].size()+tmp.size()<<" ";
		for(int j:b[i])cout<<j<<" ";
		for(int j:tmp)cout<<j<<" ";
		cout<<endl;
		qmr.clear();
		for(int j=0;j<b[i].size()+tmp.size();j++){
			int t=read();qmr+=t;
		}
		for(int j=0;j<b[i].size();j++){
//			cout<<j<<" "<<qmr[j]<<endl;
			a[b[i][j]]+=qmr[j];
		}
		basic_string<short>pmt,cnt;
		for(int j=b[i].size();j<qmr.size();j++){
			if(qmr[j]-ans[j].size()>0){
				pmt+=tmp[j-b[i].size()];
				cnt+=qmr[j]-ans[j].size();
			}
		}
//		for(int j:pmt)cout<<j<<" ";
//		cout<<endl;
		work(i,0,b[i].size()-1,pmt,cnt);
		tmp+=b[i];b[i].clear();
	}
//	if(n==4000)return 0;
	int m=0;
	for(int i=1;i<=n;i++){
		cout<<"? 2 "<<i<<" "<<i<<endl;
		int t=read();a[i]=read();
		for(int j=1;j<=a[i];j++)ans[i]+=i;
		m+=ans[i].size();
	}
	cout<<"! "<<m<<" ";
	for(int i=1;i<=n;i++){
		for(int j:ans[i])cout<<i<<" "<<j<<" ";
	}
	cout<<endl;
	return 0;
}