def eratosthenes(n):
    res=[0 for i in range(n+1)]
    prime=set([])
    for i in range(2,n+1):
        if not res[i]:
            prime.add(i)
            for j in range(1,n//i+1):
                res[i*j]=1
    return prime

prime = eratosthenes(32000)

def factorization(n):
    res = []
    for p in prime:
        cnt = 0
        while n%p==0:
            cnt += 1
            n //= p
        if cnt:
            res.append((p,cnt))
    if n!=1:
        res.append((n,1))
    return res

def divisors_prime(prime):
    res = [1]
    for p,cnt in prime:
        newres = []
        for d in res:
            t = d
            for j in range(cnt+1):
                newres.append(t)
                t *= p
        res = newres
    return res

import sys,random,bisect
from collections import deque,defaultdict
from heapq import heapify,heappop,heappush
from itertools import permutations
from math import gcd,log

input = lambda :sys.stdin.readline().rstrip()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

M = 10**9

N = int(input())
A = li()

all_prime = {}
all_divisors = set()
for a in A:
    a_prime = factorization(a)
    for p,cnt in a_prime:
        if p not in all_prime:
            all_prime[p] = 0
        all_prime[p] = max(all_prime[p],pow(p,cnt))
    div = divisors_prime(a_prime)
    for d in div:
        all_divisors.add(d)

P = [all_prime[p] for p in all_prime]
P.sort()
n = len(P)

deq = deque([(1,-1)])
break_point = -1
while deq:
    v,k = deq.popleft()
    for i in range(k+1,n):
        tmp = v * P[i]
        if M < tmp or tmp not in all_divisors:
            break_point = tmp
            deq = []
            break
        deq.append((tmp,i))

if break_point==-1:
    exit(print(-1))

selected_P = []
for i in range(n):
    if break_point%P[i]==0:
        selected_P.append(P[i])

k = len(selected_P)
S = [[] for i in range(k)]
for i in range(N):
    for j in range(k):
        if A[i]%selected_P[j]:
            S[j].append(A[i])
            break

print(k)
for i in range(k):
    S[i] = [len(S[i])] + S[i]
    print(*S[i])