#include <iostream>
#include <cstdio> 
#include <algorithm>
#include <vector>
using namespace std;
namespace Grader{
	int e[4005][4005];
	int _cnt;
	vector<int> ask(vector<int> a){
		int l=(int)a.size();
		_cnt+=l;
		vector<int> ans(l);
		for (int i=0;i<l;i++)
			for (int j=0;j<i;j++)
				ans[i]+=e[min(a[i],a[j])][max(a[i],a[j])];
		return ans;
	}
	int cnt[4005][4005];
	vector<pair<int,int>> qwq[10005];
	void work(int now,int l,int r){
		if (l==r){
			for (int i=0;i<(int)qwq[now].size();i++)
				cnt[l][qwq[now][i].first]=qwq[now][i].second;
			cnt[l][l]=ask({l,l})[1];
			return;
		}
		sort(qwq[now].begin(),qwq[now].end());
		int mid=(l+r)/2;
		vector<int> ovo;
		for (int i=l;i<=mid;i++)ovo.push_back(i);
		for (int i=0;i<(int)qwq[now].size();i++)ovo.push_back(qwq[now][i].first);
		vector<int> awa=ask(ovo);
		for (int i=0;i<(int)qwq[now].size();i++){
			int p=awa[mid-l+1+i];
			for (int j=0;j<i;j++)p-=cnt[qwq[now][j].first][qwq[now][i].first];
			if (p>0)qwq[now*2].push_back(make_pair(qwq[now][i].first,p));
			if (p<qwq[now][i].second)qwq[now*2+1].push_back(make_pair(qwq[now][i].first,qwq[now][i].second-p));
		}
		work(now*2+1,mid+1,r);
		ovo.clear(),awa.clear();
		for (int i=l;i<=r;i++)ovo.push_back(i);
		awa=ask(ovo);
		for (int i=0;i<r-mid;i++){
			int p=awa[mid-l+1+i];
			for (int j=0;j<i;j++)p-=cnt[mid+1+j][mid+1+i];
			if (p>0)qwq[now*2].push_back(make_pair(mid+1+i,p));
		}
		work(now*2,l,mid);
		return;
	} 
	void main(){
		freopen("7510.in","r",stdin);
		freopen("7510.out","w",stdout);
		int n,m;
		cin>>n>>m;
		for (int i=1;i<=m;i++){
			int a,b;
			cin>>a>>b;
			e[min(a,b)][max(a,b)]++;
		}
		work(1,1,n);
		cout<<_cnt<<endl;
		for (int i=1;i<=n;i++)
			for (int j=i;j<=n;j++)
				if (cnt[i][j]!=e[i][j])cout<<"QwQ"<<endl;
		return;
	}
}
vector<int> ask(vector<int> a){
	int l=(int)a.size();
	cout<<"? "<<l<<' ';
	for (int i=0;i<l;i++)cout<<a[i]<<' ';
	cout<<endl;
	vector<int> ans(l);
	for (int i=0;i<l;i++)cin>>ans[i];
	return ans;
}
int cnt[4005][4005];
vector<pair<int,int>> qwq[10005];
void work(int now,int l,int r){
	if (l==r){
		for (int i=0;i<(int)qwq[now].size();i++)
			cnt[l][qwq[now][i].first]=qwq[now][i].second;
		cnt[l][l]=ask({l,l})[1];
		return;
	}
	sort(qwq[now].begin(),qwq[now].end());
	int mid=(l+r)/2;
	vector<int> ovo;
	for (int i=l;i<=mid;i++)ovo.push_back(i);
	for (int i=0;i<(int)qwq[now].size();i++)ovo.push_back(qwq[now][i].first);
	vector<int> awa=ask(ovo);
	for (int i=0;i<(int)qwq[now].size();i++){
		int p=awa[mid-l+1+i];
		for (int j=0;j<i;j++)p-=cnt[qwq[now][j].first][qwq[now][i].first];
		if (p>0)qwq[now*2].push_back(make_pair(qwq[now][i].first,p));
		if (p<qwq[now][i].second)qwq[now*2+1].push_back(make_pair(qwq[now][i].first,qwq[now][i].second-p));
	}
	work(now*2+1,mid+1,r);
	ovo.clear(),awa.clear();
	for (int i=l;i<=r;i++)ovo.push_back(i);
	awa=ask(ovo);
	for (int i=0;i<r-mid;i++){
		int p=awa[mid-l+1+i];
		for (int j=0;j<i;j++)p-=cnt[mid+1+j][mid+1+i];
		if (p>0)qwq[now*2].push_back(make_pair(mid+1+i,p));
	}
	work(now*2,l,mid);
	return;
} 
int main(){
	int n;
	cin>>n;
	work(1,1,n);
	int m=0;
	for (int i=1;i<=n;i++)
		for (int j=i;j<=n;j++)
			for (int k=1;k<=cnt[i][j];k++)
				m++;
	cout<<"! "<<m<<' ';
	for (int i=1;i<=n;i++)
		for (int j=i;j<=n;j++)
			for (int k=1;k<=cnt[i][j];k++)
				cout<<i<<' '<<j<<' ';
	cout<<endl;
	return 0;
}