#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int P=998244353;
const int N=1200103;
ll f[N],a[N],pw[N];

int qry(vector<int>&v) {
    cout<<"? "<<endl;
    for(auto x:v)cout<<x<<' ';
    cout<<endl;
    int x;cin>>x;
    return x;
}

void solve(){
    int n;
    cin>>n;
    random_device r;
    mt19937 gen(r());
    uniform_int_distribution <int>rng(0,1);
    for(int i=1;i<=n;i++)f[i]=rng(gen);
    int flg=1;
    int e=0;
    vector<int>v(n);
    iota(v.begin(),v.end(),1);
    e=qry(v);
    for(int i=1;i<30;i++) {
        vector<int>o[2];
        for(int j=1;j<=n;j++)o[f[i]].push_back(i);
        int g=e-qry(o[0])-qry(o[1]);
        flg&=g%2==0;
    }
    if(!flg) {
        cout<<"! NO"<<endl;
    }else {
        cout<<"! YES"<<endl;
    }
}

int main(){
    ios::sync_with_stdio(false);
    int T=1;
    //cin>>T;
    while(T--)solve();

}