from fractions import gcd
from datetime import date, timedelta
from heapq import*
import math
from collections import defaultdict, Counter, deque
import sys
from bisect import *
import itertools
import copy
sys.setrecursionlimit(10 ** 7)
MOD = 10 ** 9 + 7


def main():
    n = int(input())
    ansl = []
    for i in range(1,n+1):
        q = []
        for j in range(1,n+1):
            if i != j:
                q.append(j)
        print("?" + " " + str(len(q)))
        print(*q)
        v = int(input())
        if v == 0:
            ansl.append(i)
    
    print("!" + " " + str(len(ansl)))
    print(*ansl)




    


if __name__ == '__main__':
    main()