結果
問題 |
No.14 最小公倍数ソート
|
ユーザー |
|
提出日時 | 2024-08-04 16:48:49 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 78 ms / 5,000 ms |
コード長 | 1,143 bytes |
コンパイル時間 | 2,246 ms |
コンパイル使用メモリ | 207,860 KB |
最終ジャッジ日時 | 2025-02-23 21:01:03 |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
#include <bits/stdc++.h> using namespace std; int main() { ios_base::sync_with_stdio(false); cin.tie(nullptr); int N; cin >> N; vector<int> A(N); for(auto &a : A) cin >> a; vector<vector<int>> Ds(10001); for(int i=1; i<=10000; i++) for(int k=i; k<=10000; k+=i) Ds.at(k).push_back(i); vector<set<pair<int,int>>> S(10001); for(int i=1; i<N; i++) for(auto &d : Ds.at(A.at(i))) S.at(d).insert({A.at(i),i}); vector<int> rev(N),P(N); iota(rev.begin(),rev.end(),0); for(int i=1; i<N; i++){ int a = A.at(P.at(i-1)); int lcm = 1001001001,B = 1001001001,pos = -1; for(auto &d : Ds.at(a)){ if(S.at(d).size() == 0) continue; auto &[b,p] = (*S.at(d).begin()); int now = a*b/d; if(now < lcm) lcm = now,B = b,pos = p; else if(now == lcm && b < B) lcm = now,B = b,pos = p; } for(auto &d : Ds.at(A.at(pos))) S.at(d).erase({A.at(pos),pos}); int p1 = P.at(i); P.at(i) = pos; rev.at(p1) = rev.at(pos); } for(int i=0; i<N; i++) cout << A.at(P.at(i)) << (i==N-1?"\n":" "); }