結果

問題 No.1611 Minimum Multiple with Double Divisors
ユーザー nebukuro09nebukuro09
提出日時 2021-07-21 22:23:30
言語 D
(dmd 2.105.2)
結果
WA  
実行時間 -
コード長 1,747 bytes
コンパイル時間 861 ms
コンパイル使用メモリ 102,284 KB
実行使用メモリ 6,472 KB
最終ジャッジ日時 2023-09-04 13:22:15
合計ジャッジ時間 10,369 ms
ジャッジサーバーID
(参考情報)
judge14 / judge13
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 859 ms
5,728 KB
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 AC 10 ms
4,380 KB
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
testcase_35 WA -
testcase_36 WA -
testcase_37 WA -
testcase_38 WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

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;

void main() {
    auto T = readln.chomp.to!int;

    long[] P;
    auto table = new bool[](10^^6);
    for (long i = 2; i * i <= 10L^^11; ++i) {
        if (table[i.to!int]) continue;
        P ~= i;
        for (long j = i * 2; j * j <= 10L ^^ 11; j += i) {
            table[j.to!int] = true;
        }
    }

    long[long] Q;

    while (T--) {
        auto X = readln.chomp.to!long;
        Q.clear();
        foreach (i; 0..P.length.to!int) {
            if (X % P[i] != 0) {
                long acm = 1;
                foreach (c; Q.values()) acm *= c + 1;
                long ans = min(X * P[i], X * dfs(1, X, acm, Q.keys().array, Q, P[i], acm));
                ans.writeln;
                break;
            } else {
                long count = 0;
                long Y = X;
                while (Y % P[i] == 0) {
                    Y /= P[i];
                    count += 1;
                }
                Q[P[i]] = count;
            }
        }
    }
}

long dfs(long cur, long X, long original, long[] P, long[long] count, long maxP, long acm) {
    if (cur >= maxP) {
        return 1L << 59;
    }
    if (original * 2 == acm) {
        return cur;
    }
    if (original * 2 < acm) {
        return 1L << 59;
    }
    long ans = 1L << 59;
    foreach (p; P) if (cur * p <= maxP) {
        acm /= count[p] + 1;
        count[p] += 1;
        acm *= count[p] + 1;
        ans = min(ans, dfs(cur * p, X, original, P, count, maxP, acm));
        acm /= count[p] + 1;
        count[p] -= 1;
        acm *= count[p] + 1;
    }
    return ans;
}
0