結果

問題 No.14 最小公倍数ソート
ユーザー not_522not_522
提出日時 2015-08-20 13:49:25
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 97 ms / 5,000 ms
コード長 2,529 bytes
コンパイル時間 1,948 ms
コンパイル使用メモリ 161,012 KB
実行使用メモリ 8,572 KB
最終ジャッジ日時 2023-09-25 13:48:25
合計ジャッジ時間 4,203 ms
ジャッジサーバーID
(参考情報)
judge14 / judge12
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 10 ms
4,376 KB
testcase_04 AC 97 ms
8,564 KB
testcase_05 AC 48 ms
6,704 KB
testcase_06 AC 55 ms
6,712 KB
testcase_07 AC 63 ms
7,076 KB
testcase_08 AC 79 ms
7,420 KB
testcase_09 AC 87 ms
8,020 KB
testcase_10 AC 91 ms
8,168 KB
testcase_11 AC 91 ms
8,552 KB
testcase_12 AC 92 ms
8,276 KB
testcase_13 AC 94 ms
8,288 KB
testcase_14 AC 96 ms
8,272 KB
testcase_15 AC 97 ms
8,572 KB
testcase_16 AC 57 ms
6,360 KB
testcase_17 AC 48 ms
6,168 KB
testcase_18 AC 32 ms
5,116 KB
testcase_19 AC 74 ms
7,536 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#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