結果

問題 No.14 最小公倍数ソート
コンテスト
ユーザー nebukuro09
提出日時 2018-09-29 02:03:57
言語 D
(dmd 2.109.1)
結果
WA  
実行時間 -
コード長 2,098 bytes
コンパイル時間 1,949 ms
コンパイル使用メモリ 162,976 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-06-13 01:42:00
合計ジャッジ時間 2,977 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3 WA * 17
権限があれば一括ダウンロードができます
コンパイルメッセージ
/home/linuxbrew/.linuxbrew/opt/dmd/include/dlang/dmd/std/numeric.d(2999): Warning: cannot inline function `std.numeric.gcdImpl!uint.gcdImpl`

ソースコード

diff #
raw source code 

import std.stdio, std.array, std.string, std.conv, std.algorithm;
import std.typecons, std.range, std.random, std.math, std.container;
import std.numeric, std.bigint, core.bitop, std.bitmanip;

immutable int MAX = 100001;

void main() {
    auto N = readln.chomp.to!int;
    auto A = readln.split.map!(to!int).array;
    int[][int] P;
    RedBlackTree!(int, "a < b", true)[int] B;

    auto rbt = new RedBlackTree!(int, "a < b", true)();
    foreach (i; 1..N) rbt.insert(A[i]);

    foreach (i; 0..N) {
        if (A[i] in P) continue;
        int x = A[i];
        for (int j = 2; j * j <= A[i]; ++j) {
            if (x % j == 0) P[A[i]] ~= j;
            while (x % j == 0) x /= j;
        }
        if (x > 1) P[A[i]] ~= x;
    }

    foreach (i; 1..N) {
        if (A[i] == 1) continue;
        foreach (p; P[A[i]]) {
            if (p !in B)
                B[p] = new RedBlackTree!(int, "a < b", true)();
            B[p].insert(A[i]);
        }
    }

    int[] ans = [A[0]];
    int cur = A[0];

    foreach (i; 0..N-1) {
        if ((1 in rbt) || cur == 1) {
            int x = rbt.front;
            rbt.removeKey(x);
            if (x in P) foreach (p; P[x]) B[p].removeKey(x);
            ans ~= x;
            cur = x;
            continue;
        }
        int min_lcm = 1 << 29;
        int min_num = -1;
        foreach (p; P[cur]) {
            if (B[p].empty) continue;
            int x = B[p].front;
            int lcm = x * cur / gcd(x, cur);
            if (lcm < min_lcm) {
                min_lcm = lcm;
                min_num = x;
            } else if (lcm == min_lcm && x < min_num) {
                min_num = x;
            }
        }
        int x = rbt.front;
        int lcm = x * cur / gcd(x, cur);
        if (lcm < min_lcm) {
            min_lcm = lcm;
            min_num = x;
        } else if (lcm == min_lcm && x < min_num) {
            min_num = x;
        }

        x = min_num;
        ans ~= x;
        if (x in P) foreach (p; P[x]) B[p].removeKey(x);
        rbt.removeKey(x);
        cur = x;
    }

    ans.map!(to!string).join(" ").writeln;
}
0