#include<bits/stdc++.h>
using namespace std;
#define For(i,l,r) for(int i=(int)(l);i<=(int)(r);i++)
#define pb push_back
const int N=4005;
int n,ans[N][N];
vector<int> ask(vector<int> v){
    cout<<"? "<<v.size();
    for(auto i:v)cout<<" "<<i;
    cout<<endl;
    vector<int> res=v;
    for(auto &i:res)cin>>i;
    return res;
}
vector<int> ask(vector<int> s1,vector<int> s2){

    vector<int> tmp=s1;
    for(auto i:s2)tmp.pb(i);
    auto res=ask(tmp);
    auto zzz=s2;
    For(i,0,s2.size()-1){
        int x=s2[i],tmp=res[i+s1.size()];
        For(j,0,i-1)if(res[s1.size()+j]==0)tmp-=ans[s2[j]][s2[i]];
        zzz[i]=tmp;
    }
    return zzz;
}
void solve(vector<int> s1,vector<int> s2,vector<int> s3){
    if(s1.size()==0||s2.size()==0)return;
    if(s1.size()==1){
        For(o,0,s2.size()-1)ans[s2[o]][s1[0]]=ans[s1[0]][s2[o]]=s3[o]; return;
    }
    vector<int> ls1,rs1,ls2,rs2,ls3,rs3;
    int mid=s1.size()/2;
    For(i,0,mid-1)ls1.pb(s1[i]); 
    For(i,mid,s1.size()-1)rs1.pb(s1[i]);
    auto t=ask(ls1,s2);
    For(i,0,s2.size()-1)if(t[i]==0){
        rs2.pb(s2[i]); rs3.pb(s3[i]);
    }else if(t[i]==s3[i]){
        ls2.pb(s2[i]); ls3.pb(s3[i]);
    }else{
        ls2.pb(s2[i]); ls3.pb(t[i]);
        rs2.pb(s2[i]); rs3.pb(s3[i]-t[i]);
    }
    solve(ls1,ls2,ls3);
    solve(rs1,rs2,rs3);
}
void solve(vector<int> v){
    if(v.size()<=1)return;
    auto res=ask(v);
    vector<int> s1,s2;
    For(i,0,v.size()-1){
        if(res[i]==0)s1.pb(v[i]); else s2.pb(v[i]);
    }
    solve(s2);
    solve(s1,s2,ask(s1,s2));
}
int main(){

    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    cin>>n;
    vector<int> v;
    For(i,1,n)v.pb(i);
    solve(v);
    For(i,1,n){
        ans[i][i]=ask({i,i})[1];
    }
    vector<pair<int,int>> sb;
    For(i,1,n)For(j,i,n)For(k,1,ans[i][j])sb.pb(make_pair(i,j));
    cout<<"! "<<sb.size();
    for(auto i:sb)cout<<" "<<i.first<<" "<<i.second;
    cout<<endl;
}