/*
cat <<EOF >mistaken-paste
*/

#pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
#define _USE_MATH_DEFINES

#include <stdio.h>
#include <stdint.h>
#include <stdlib.h>
#include <string.h>
#include <stdbool.h>
#include <math.h>
#include <time.h>

#define BIG 2000000007
#define VERYBIG 2000000000000007LL

#define MOD 1000000007
#define FOD  998244353
typedef uint64_t ull;
typedef  int64_t sll;

#define N_MAX 1000000
#define M_MAX 200000

#ifdef __cplusplus
#include <queue>
#include <stack>
#include <tuple>
#include <set>
#include <map>
#include <string>
#include <algorithm>
#include <functional>
#include <array>

using std::queue;
using std::priority_queue;
using std::stack;
using std::tuple;
using std::set;
using std::map;
using std::vector;
using std::greater;
using std::pair;
using std::string;

#endif

typedef struct {
	int32_t a;
	int32_t b;
} hw;

typedef struct {
	sll a;
	sll b;
} hwll;

typedef struct {
	sll a;
	sll b;
	sll c;
} hwllc;

typedef struct {
	hwll a;
	hwll b;
} linell;

ull n, m;
ull h, w;
ull k;
ull q;
sll va, vb, vc, vd, ve, vf;
ull ua, ub, uc, ud, ue, uf;
long double vra, vrb, vrc;
double vda, vdb, vdc;
char ch, dh;

ull umin (ull x, ull y) {
	return (x < y) ? x : y;
}

ull umax (ull x, ull y) {
	return (x > y) ? x : y;
}

sll smin (sll x, sll y) {
	return (x < y) ? x : y;
}

sll smax (sll x, sll y) {
	return (x > y) ? x : y;
}

ull gcd (ull x, ull y) {
	if (x < y) {
		return gcd(y, x);
	} else if (y == 0) {
		return x;
	} else {
		return gcd(y, x % y);
	}
}

ull bitpow (ull a, ull x, ull modulo) {
	ull result = 1;
	while (x) {
		if (x & 1) {
			result *= a;
			result %= modulo;
		}
		x /= 2;
		a = (a * a) % modulo;
	}
	return result;
}

ull divide (ull a, ull b, ull modulo) {
	return (a * bitpow(b, modulo - 2, modulo)) % modulo;
}

ull udiff (ull a, ull b) {
	if (a >= b) {
		return a - b;
	} else {
		return b - a;
	}
}

sll sdiff (sll a, sll b) {
	if (a >= b) {
		return a - b;
	} else {
		return b - a;
	}
}

int bitcount (ull n) {
	int result = 0;
	while (n) {
		if (n & 1) result++;
		n /= 2;
	}
	return result;
}

// double distance (sll x1, sll y1, sll x2, sll y2) {
// 	double xdist2, ydist2, origindist, dist;

// 	xdist2 = (x1 - x2) * (x1 - x2);
// 	ydist2 = (y1 - y2) * (y1 - y2);
// 	return sqrt(xdist2 + ydist2);
// }

int32_t pullcomp (const void *left, const void *right) {
	ull l = *(ull*)left;
	ull r = *(ull*)right;
	if (l < r) {
		return -1;
	}
	if (l > r) {
		return +1;
	}
	return 0;
}

int32_t psllcomp (const void *left, const void *right) {
	sll l = *(sll*)left;
	sll r = *(sll*)right;
	if (l < r) {
		return -1;
	}
	if (l > r) {
		return +1;
	}
	return 0;
}

int32_t pcharcomp (const void *left, const void *right) {
	char l = *(char*)left;
	char r = *(char*)right;
	if (l < r) {
		return -1;
	}
	if (l > r) {
		return +1;
	}
	return 0;
}

int32_t pdoublecomp (const void *left, const void *right) {
	double l = *(double*)left;
	double r = *(double*)right;
	if (l < r) {
		return -1;
	}
	if (l > r) {
		return +1;
	}
	return 0;
}

int32_t pstrcomp (const void *left, const void *right) {
	char* l = *(char**)left;
	char* r = *(char**)right;

	return strcmp(l, r);
}

int32_t phwllABcomp (const void *left, const void *right) {
	hwll l = *(hwll*)left;
	hwll r = *(hwll*)right;
	if (l.a < r.a) {
		return -1;
	}
	if (l.a > r.a) {
		return +1;
	}
	if (l.b < r.b) {
		return -1;
	}
	if (l.b > r.b) {
		return +1;
	}
	return 0;
}

int32_t phwllREVcomp (const void *left, const void *right) {
	hwll l = *(hwll*)left;
	hwll r = *(hwll*)right;
	if (l.b < r.b) {
		return -1;
	}
	if (l.b > r.b) {
		return +1;
	}
	if (l.a < r.a) {
		return -1;
	}
	if (l.a > r.a) {
		return +1;
	}
	return 0;
}

int32_t ptriplecomp (const void *left, const void *right) {
	hwllc l = *(hwllc*)left;
	hwllc r = *(hwllc*)right;

	if (l.a < r.a) {
		return -1;
	}
	if (l.a > r.a) {
		return +1;
	}
	if (l.b < r.b) {
		return -1;
	}
	if (l.b > r.b) {
		return +1;
	}
	if (l.c < r.c) {
		return -1;
	}
	if (l.c > r.c) {
		return +1;
	}
	return 0;
}

bool isinrange (sll left, sll x, sll right) {
	return (left <= x && x <= right);
}

bool isinrange_soft (sll left, sll x, sll right) {
	return (left <= x && x <= right) || (left >= x && x >= right);
}

sll a[N_MAX];
// sll a[3001][3001];
sll b[N_MAX];
// sll b[3001][3001];
sll c[N_MAX];
// sll d[N_MAX];
// sll e[N_MAX];
char s[N_MAX + 1];
// char s[3010][3010];
char t[N_MAX + 1];
// char t[3010][3010];
// hwll xy[N_MAX];
// sll table[1000][1000];

set<ull> g[2002];
bool iscoprime[2002][2002];
clock_t tot;

ull iswatched[2002];
void toggle (ull v, ull x) {
	iswatched[v] = x;

	for (ull u : g[v]) {
		if (iswatched[u] ^ x) toggle(u, x);
	}
}
clock_t togrec (clock_t t) {
	static clock_t total = 0;
	total += t;
	return total;
}

typedef pair<ull, ull> P;
void merge(vector<vector<ull>> &x, vector<ull> &v) {
	ull curr[x.size()];
	priority_queue<P, vector<P> > que;
	for (ull i = 0; i < x.size(); i++) {
		curr[i] = 0;
		que.push(P(x[i][0], i));
	}
	while (!que.empty()) {
		P p = que.top();
		que.pop();

		v.push_back(p.first);
		ull i = p.second;
		if (++curr[i] < x[i].size()) {
			que.push(P(x[i][curr[i]], i));
		}
	}

	// return v;
}

clock_t ft = 0;
// bool iscome[2002];
vector<ull> func (ull v) { // This can be O(NlogN)
	vector<ull> r;
	
	// clock_t times = 0;
	// times -= clock();

	// fprintf(stderr, "%lld.. (%d)\n", v, clock() - tot);
	// fflush(stderr);
	// if (iscome[v]) {
	// 	fprintf(stderr, "came %lld!\n", v);
	// 	abort();
	// }
	// iscome[v] = true;

	priority_queue<ull, vector<ull>, greater<ull> > que;
	// clock_t start = clock();
	// toggle(v, 0);
	// clock_t end = clock();
	// fprintf(stderr, "toggle(%lld, 0): %d\n", v, end - start);
	// togrec(end - start);
	iswatched[v] = 1;

	for (ull u : g[v]) {
		g[u].erase(v);
		que.push(u);
	}

	vector<vector<ull>> vs;
	while (!que.empty()) {
		ull u = que.top();
		que.pop();
		if (iswatched[u]) continue;

		// clock_t start = clock();
		// toggle(u, 1);
		// clock_t end = clock();
		// fprintf(stderr, "toggle(%lld, 1): %d\n", u, end - start);
		// fflush(stderr);
		// togrec(end - start);

		// times += clock();
		vs.push_back(func(u));
		// times -= clock();
	}

	r.push_back(v);
	merge(vs, r);

	// times += clock();
	// fprintf(stderr, "%lld: %d\n", v, times);
	// fflush(stderr);
	// ft += times;

	return r;
}

ull solve () {
	sll i, j, ki, li;
	ull result = 0;
	// sll result = 0;
	double dresult = 0;
	// ull maybe = 0;
	sll maybe = 0;
	// ull sum = 0;
	sll sum = 0;
	sll item;
	ull *dpcell;

	qsort(a, n, sizeof(ull), pullcomp);

	for (i = 0; i < n; i++) {
		for (j = i + 1; j < n; j++) {
			if (gcd(a[i], a[j]) > 1) {
				g[i].insert(j);
				g[j].insert(i);
			}
		}
	}

	vector<vector<ull>> finals;

	// tot = clock();

	for (i = 0; i < n; i++) {
		if (iswatched[i]) continue;

		// clock_t start = clock();
		// toggle(i, 1);
		// clock_t end = clock();
		// fprintf(stderr, "toggle(%lld, 1): %d\n", i, end - start);
		// togrec(end - start);

		finals.push_back(func(i));
		// for (j = 0; j < finals[finals.size() - 1].size(); j++) printf("%llu ", finals[finals.size() - 1][j]);
		// 	puts("");
		// fflush(stdout);
	}
	vector<ull> r;
	merge(finals, r);
	for (i = 0; i < n; i++) {
		printf("%llu%c", a[r[i]], (i == n - 1 ? '\n' : ' '));
	}
	// fprintf(stderr, "togrec: %d\n", togrec(0));
	// fprintf(stderr, "ft: %d\n", ft);

	// printf("%lld\n", result);
	// printf("%.15lf\n", dresult);
	// puts(s);

	return 0;

	success:
	// puts("YES");
	puts("Yes");
	// printf("%llu\n", result);
	// puts("0");
	// puts("Yay!");
	return 0;

	fail:
	// puts("NO");
	puts("No");
	// puts("0");
	// puts("-1");
	// puts("-1 -1 -1");
	// puts(":(");
	return 1;
}

int32_t main (void) {
	int32_t i, j;
	int32_t x, y;

	scanf("%llu", &n, &m);
	// scanf("%llu", &k, &n, &m);
	// scanf("%llu%llu", &h, &w);
	// scanf("%lld", &va, &vb, &vc, &vd);
	// scanf("%llu%llu", &ua, &ub, &uc, &ud);
	// scanf("%llu", &q);
	// scanf("%s", s);
	// scanf("%s", t);
	for (i = 0; i < n; i++) {
		// scanf("%lld%lld", &xy[i].a, &xy[i].b);
		scanf("%lld", &a[i]);
		// scanf("%lld", &b[i]);
		// scanf("%lld", &c[i]);
		// scanf("%lld", &d[i]);
		// a[i]--;
		// b[i]--;
		// c[i]--;
	}
	// scanf("%llu", &m, &k);
	// scanf("%llu", &q);
	// for (i = 0; i < n; i++) {
	// 	scanf("%lld", &b[i]);
	// }

	// for (i = 0; i < h; i++) {
	// 	for (j = 0; j < w; j++) {
	// 		scanf("%llu", &table[i][j]);
	// 	}
	// }
	// for (i = 0; i < h; i++) {
	// 	scanf("%s", &s[i]);
	// }

	solve();

	return 0;
}


Measured memory_bytes: 137084