結果

問題 No.14 最小公倍数ソート
コンテスト
ユーザー not_522
提出日時 2015-08-20 13:49:25
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 80 ms / 5,000 ms
コード長 2,529 bytes
コンパイル時間 1,480 ms
コンパイル使用メモリ 173,920 KB
実行使用メモリ 8,704 KB
最終ジャッジ日時 2024-07-18 10:51:13
合計ジャッジ時間 3,350 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>

using namespace std;

class Prime {
private:
  vector<long long> div;
  
public:
  Prime(long long n = 0) {
    for (long long i = 0; i <= n; ++i) div.emplace_back(i);
    for (long long i = 2; i <= n / i; ++i) if (div[i] == i) {
      for (long long j = i * i; j <= n; j += i) div[j] = i;
    }
  }
  
  bool isPrime(long long n) const {
    if (n <= 1) return false;
    if (n < (long long)div.size()) return div[n] == n;
    for (long long i = 2; i <= n / i; ++i) if (n % i == 0) return false;
    return true;
  }
  
  vector<long long> factor(long long n) const {
    vector<long long> res;
    for (long long i = 2; i <= n / i && n >= (long long)div.size(); ++i) {
      while (n % i == 0) {
        res.emplace_back(i);
        n /= i;
      }
    }
    if (n >= max((long long)div.size(), 2ll)) {
      res.emplace_back(n);
    } else {
      while (n > 1) {
        res.emplace_back(div[n]);
        n /= div[n];
      }
    }
    sort(res.begin(), res.end());
    return res;
  }
  
  vector<long long> primeFactor(long long n) const {
    vector<long long> res;
    for (long long i = 2; i <= n / i && n >= (long long)div.size(); ++i) {
      if (n % i) continue;
      res.emplace_back(i);
      while (n % i == 0) n /= i;
    }
    if (n >= max((long long)div.size(), 2ll)) {
      res.emplace_back(n);
    } else {
      while (n > 1) {
        long long d = div[n];
        res.emplace_back(d);
        while (n % d == 0) n /= d;
      }
    }
    sort(res.begin(), res.end());
    return res;
  }

  vector<long long> divisor(long long n) const {
    vector<long long> res;
    for (long long i = 1; i <= n / i; ++i) {
      if (n % i == 0) {
        res.emplace_back(i);
        if (i != n / i) res.emplace_back(n / i);
      }
    }
    sort(res.begin(), res.end());
    return res;
  }
};

int main() {
  Prime prime(10000);
  int n;
  cin >> n;
  vector<int> a(n);
  for (int& i : a) cin >> i;
  multiset<int> s(a.begin() + 1, a.end());
  vector<multiset<int>> v(10001);
  for (int i : s) {
    for (auto d : prime.divisor(i)) v[d].insert(i);
  }
  int p = a[0];
  while (!s.empty()) {
    int next = -1, lcm = numeric_limits<int>::max();
    for (auto d : prime.divisor(p)) {
      if (v[d].empty()) continue;
      int l = p * *(v[d].begin()) / d;
      if (lcm > l) {
        next = *(v[d].begin());
        lcm = l;
      }
    }
    cout << p << " ";
    p = next;
    s.erase(s.find(p));
    for (auto d : prime.divisor(p)) v[d].erase(v[d].find(p));
  }
  cout << p << endl;
}
0