結果

問題 No.14 最小公倍数ソート
ユーザー not_522not_522
提出日時 2015-08-20 13:49:25
言語 C++11
(gcc 11.4.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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 1 ms
5,376 KB
testcase_03 AC 9 ms
5,376 KB
testcase_04 AC 80 ms
8,704 KB
testcase_05 AC 39 ms
6,528 KB
testcase_06 AC 44 ms
6,784 KB
testcase_07 AC 53 ms
7,296 KB
testcase_08 AC 62 ms
7,808 KB
testcase_09 AC 73 ms
8,320 KB
testcase_10 AC 73 ms
8,320 KB
testcase_11 AC 75 ms
8,448 KB
testcase_12 AC 78 ms
8,320 KB
testcase_13 AC 76 ms
8,448 KB
testcase_14 AC 75 ms
8,448 KB
testcase_15 AC 77 ms
8,448 KB
testcase_16 AC 46 ms
6,784 KB
testcase_17 AC 39 ms
6,272 KB
testcase_18 AC 26 ms
5,376 KB
testcase_19 AC 61 ms
7,680 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