import os
import sys
from io import BytesIO, IOBase

BUFSIZE = 8192


class FastIO(IOBase):
    newlines = 0

    def __init__(self, file):
        self._fd = file.fileno()
        self.buffer = BytesIO()
        self.writable = "x" in file.mode or "r" not in file.mode
        self.write = self.buffer.write if self.writable else None

    def read(self):
        while True:
            b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE))
            if not b:
                break
            ptr = self.buffer.tell()
            self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr)
        self.newlines = 0
        return self.buffer.read()

    def readline(self):
        while self.newlines == 0:
            b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE))
            self.newlines = b.count(b"\n") + (not b)
            ptr = self.buffer.tell()
            self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr)
        self.newlines -= 1
        return self.buffer.readline()

    def flush(self):
        if self.writable:
            os.write(self._fd, self.buffer.getvalue())
            self.buffer.truncate(0), self.buffer.seek(0)


class IOWrapper(IOBase):
    def __init__(self, file):
        self.buffer = FastIO(file)
        self.flush = self.buffer.flush
        self.writable = self.buffer.writable
        self.write = lambda s: self.buffer.write(s.encode("ascii"))
        self.read = lambda: self.buffer.read().decode("ascii")
        self.readline = lambda: self.buffer.readline().decode("ascii")


sys.stdin, sys.stdout = IOWrapper(sys.stdin), IOWrapper(sys.stdout)
input = lambda: sys.stdin.readline().rstrip("\r\n")
##########################################################
from collections import Counter
# c=sorted((i,int(val))for i,val in enumerate(input().split()))
import heapq
import math
# c=sorted((i,int(val))for i,val in enumerate(input().split()))
#n=int(input())
#arr = list(map(int, input().split()))
#n,x= map(int, input().split())
def f(x):
    d=[]
    for i in range(1,int(n//2)+1):
        if x%i==0:
            d.append(i)
    #print(d)
    return d


for _ in range(int(input())):
    n=int(input())
    arr = list(map(int, input().split()))
    c=Counter(arr)

    ans=f(max(arr))
    ans.sort(reverse=True)
    p=[]
    var=max(arr)
    p.extend([var*c[var]])
    #print(p)
    c[var]=0
    #print(c)
    for i in ans:
        for j in c:
            if c[j]>0 and (i==j or j%i==0):
                p.extend(*[[j]*c[j]])
                c[j]=0
        #print("a",p)
    #rr=[]
    #for i in c:
        #if c[i]>0:
            #rr.extend([i*c[i]])
    #rr.sort(reverse=True)
    #p.extend(rr)
    print(*p)
    #var=p[0]
    #for i in  range(n):
        #print(math.gcd(var,p[i]))

            #print(i,end=" ")





















