/* O(X log log X + Σ_i σ(a[i])) where X = 10000 (a[i]の上限), σ(x) = xの約数の数 (なはず)! シンプルに O(N + X^{1+\epsilon}) ということもできるはず。 */ #define _CRT_SECURE_NO_WARNINGS #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define rep(i,n) for(int (i)=0;(i)<(int)(n);++(i)) #define rer(i,l,u) for(int (i)=(int)(l);(i)<=(int)(u);++(i)) #define reu(i,l,u) for(int (i)=(int)(l);(i)<(int)(u);++(i)) #if defined(_MSC_VER) || __cplusplus > 199711L #define aut(r,v) auto r = (v) #else #define aut(r,v) typeof(v) r = (v) #endif #define each(it,o) for(aut(it, (o).begin()); it != (o).end(); ++ it) #define all(o) (o).begin(), (o).end() #define pb(x) push_back(x) #define mp(x,y) make_pair((x),(y)) #define mset(m,v) memset(m,v,sizeof(m)) #define INF 0x3f3f3f3f #define INFL 0x3f3f3f3f3f3f3f3fLL using namespace std; typedef vector vi; typedef pair pii; typedef vector > vpii; typedef long long ll; typedef vector vl; typedef pair pll; typedef vector > vpll; typedef vector vs; typedef long double ld; template inline void amin(T &x, U y) { if(y < x) x = y; } template inline void amax(T &x, U y) { if(x < y) x = y; } vector smallestPrimeFactor; void getSmallestPrimeFactors(int n) { n = max(n, 1); if((int)smallestPrimeFactor.size() >= n+1) return; smallestPrimeFactor.assign(n+1, 0); smallestPrimeFactor[1] = 1; rer(i, 2, n) if(!smallestPrimeFactor[i]) { smallestPrimeFactor[i] = i; if((ll)i * i <= n) for(int j = i * i; j <= n; j += i) if(!smallestPrimeFactor[j]) smallestPrimeFactor[j] = i; } } typedef int FactorsInt; typedef vector > Factors; void primeFactors(FactorsInt x, Factors &out_v) { getSmallestPrimeFactors(x); out_v.clear(); while(x != 1) { int p = smallestPrimeFactor[x], k = 0; x /= p, k ++; while(x % p == 0) x /= p, k ++; out_v.push_back(make_pair(p, k)); } } void getDivisors(FactorsInt x, vector &out_v) { Factors fs; primeFactors(x, fs); out_v.assign(1, 1); rep(i, fs.size()) { for(int j = (int)out_v.size()-1; j >= 0; j --) { FactorsInt x = out_v[j]; rep(k, fs[i].second) { x *= fs[i].first; out_v.push_back(x); } } } // sort(all(out_v)); } struct Node { Node *next, *prev; Node(): next(0), prev(0) { } }; void insertToList(Node *&head, Node *t) { if(t->next = head) head->prev = t; head = t; } void removeFromList(Node *&head, Node *t) { if(t->prev) t->prev->next = t->next; if(t->next) t->next->prev = t->prev; if(head == t) head = t->next; t->prev = t->next = 0; } int main() { int N; scanf("%d", &N); vector a(N); rep(i, N) scanf("%d", &a[i]); const int X = 10000; getSmallestPrimeFactors(X); vector divisors(N); rep(i, N) getDivisors(a[i], divisors[i]); vector indices(X+1); rep(i, N) indices[a[i]].push_back(i); vector offsets(N); int numTotalDivisors = 0; rep(i, N) { offsets[i] = numTotalDivisors; numTotalDivisors += divisors[i].size(); } vector invoffsets(numTotalDivisors); rep(i, N) rep(j, divisors[i].size()) invoffsets[offsets[i] + j] = i; vector nodes(numTotalDivisors); vector heads(X+1, nullptr); for(int x = X; x >= 1; -- x) each(it, indices[x]) { int i = *it; rep(j, divisors[i].size()) { int d = divisors[i][j]; insertToList(heads[d], &nodes[offsets[i] + j]); } } int curi = 0; printf("%d", a[curi]); rep(k, N-1) { rep(j, divisors[curi].size()) { int d = divisors[curi][j]; removeFromList(heads[d], &nodes[offsets[curi] + j]); } pair p(mp(INF, INF), -1); rep(j, divisors[curi].size()) { int d = divisors[curi][j]; Node *t = heads[d]; if(t != nullptr) { int i = invoffsets[t - &nodes[0]]; amin(p, mp(mp(a[i] / d, a[i]), i)); } } curi = p.second; printf(" %d", a[curi]); } puts(""); return 0; }