#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
int gcd(int a, int b){
	if(b==0) return a;
	return gcd(b, a%b);
}
int main()
{
	int n;
	cin>>n;
	int a[10000];
	bool ok[10001]={};
	vector<int> d[10001];
	multiset<int> st[10001];
	for(int i=0; i<n; i++){
		cin>>a[i];
		for(int j=1; j*j<=a[i]; j++){
			if(a[i]%j==0){
				if(!ok[a[i]]) d[a[i]].push_back(j);
				st[j].insert(a[i]);
				if(j*j<a[i]){
					if(!ok[a[i]]) d[a[i]].push_back(a[i]/j);
					st[a[i]/j].insert(a[i]);
				}
			}
		}
		ok[a[i]]=1;
	}
	int ans[10000];
	ans[0]=a[0];
	for(auto x:d[a[0]]){
		auto itr=st[x].lower_bound(a[0]);
		st[x].erase(itr);
	}
	for(int i=1; i<n; i++){
		int mnl=10001, mna;
		for(auto x:d[ans[i-1]]){
		    if(st[x].empty()) continue;
			int a1=*(st[x].begin());
			if(mnl>a1/x){
				mnl=a1/x;
				mna=a1;
			}
		}
		ans[i]=mna;
		for(auto x:d[mna]){
			auto itr=st[x].lower_bound(mna);
			st[x].erase(itr);
		}
	}
	for(int i=0; i<n; i++){
		cout<<ans[i];
		if(i<n-1) cout<<" ";
	}
	cout<<endl;
	return 0;
}