#include <bits/stdc++.h>

typedef long long ll;
typedef unsigned long long ull;

#define FOR(i,a,b) for(int (i)=(a);i<(b);i++)
#define REP(i,n) FOR(i,0,n)
#define RANGE(vec) (vec).begin(),(vec).end()

using namespace std;

template<typename T>
T gcd(T a, T b) {
    if ( a < b )
        std::swap(a,b);
    if ( b == 0 )
        return a;
    return gcd(b, a%b);
}
template<typename T>
T lcm(T a, T b) {
        return a*b/gcd(a,b);
}

class LCMsort {
public:
    void solve(void) {
            int N;
            cin>>N;
            vector<int> a(N,0);
            REP(i,N)
                cin>>a[i];
            vector<bool> used(N,false);

            // sort を進めるにつれソート対象列が短くなるので、実際にソートするのでなく、lcm が最も小さい
            // ものを取り出していくだけでよい。
            // O(N^2*A)
            REP(pivot,N)
            {
                cout<<a[pivot]<<" ";
                int mx = (1<<30);
                int mi = -1;
                // O(N*A) pivot が更新されるにつれこのループ速度は高速化されていくはず。
                FOR(i, pivot+1, N)
                {
                    int k = lcm(a[pivot], a[i]);
                    if (k < mx || (k == mx && a[i] < a[mi]))
                    {
                        mx = k;
                        mi = i;
                    }
                }
                if (mi < 0)
                    break;
                // pivot を入れ替えることで FOR(i, pivot+1, N) のループ回数を減らせる
                swap(a[pivot+1], a[mi]);
            }
            cout<<endl;
    }
};

#if 1
int main(int argc, char *argv[])
{
        ios::sync_with_stdio(false);
        auto obj = new LCMsort();
        obj->solve();
        delete obj;
        return 0;
}
#endif