#include<bits/stdc++.h>
using namespace std;
int n,m,ct[5010][5010];
#define vi vector<int>
#define pb push_back
int All;
vi Ask(vi a){
    printf("? %d ",a.size());
    for(int x:a)printf("%d ",x);puts("");fflush(stdout);
    int sz=a.size();
    vi b;
    for(int i=0;i<sz;i++){int e;scanf("%d",&e);b.pb(e);}
    All+=sz;
    if(All>176000)assert(0);
    return b;
}
#define pi pair<int,int>
#define mp make_pair
#define fi first
#define se second
void cdq(vi a,vector<pi>b){
    if(b.empty())return;
    int sa=a.size();
    if(!sa)return;
    if(sa==1){
        for(pi o:b)ct[a[0]][o.fi]=ct[o.fi][a[0]]=o.se;
        return;
    }
    vi cl,cr;int mi=(sa-1)/2;
    for(int i=0;i<=mi;i++)cl.pb(a[i]);
    for(int i=mi+1;i<sa;i++)cr.pb(a[i]);
    vi a_=cl;int oi=cl.size(),bs=b.size();
    for(pi x:b)a_.pb(x.fi);
    vi E=Ask(a_);
    vector<pi>tl,tr;
    for(int i=0;i<bs;i++){
        int z=E[oi+i],x=b[i].fi;
        for(int j=0;j<i;j++)if(!E[oi+j])z-=ct[x][b[j].fi];
        if(z)tl.pb(mp(x,z));
        if(z!=b[i].se)tr.pb(mp(x,b[i].se-z));
    }
    cdq(cl,tl),cdq(cr,tr);
}
void sol(vi V){
    if(V.size()<=1)return;
    vi A,B,E=Ask(V);
    int sz=V.size();
    for(int i=0;i<sz;i++){
        if(E[i])B.pb(V[i]);
        else A.pb(V[i]);
    }
    sol(B);V.clear();
    for(int x:A)V.pb(x);for(int x:B)V.pb(x);
    int sa=A.size(),sb=B.size();
    E=Ask(V);
    vector<pi>C;
    for(int x:B)C.pb(mp(x,E[sa++]));
    cdq(A,C);
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        vi zh=Ask({i,i});
        ct[i][i]=zh[1];
    }
    vector<int>V;for(int i=1;i<=n;i++)V.pb(i);
    sol(V);
    int ps=0;
    for(int i=1;i<=n;i++)for(int j=i;j<=n;j++)ps+=ct[i][j];
    printf("! %d ",ps);
    for(int i=1;i<=n;i++)for(int j=i;j<=n;j++)while(ct[i][j]--)printf("%d %d ",i,j);
    puts("");fflush(stdout);
    return 0;
}