結果

問題 No.376 立方体のN等分 (2)
ユーザー te-shte-sh
提出日時 2017-07-24 14:18:57
言語 D
(dmd 2.106.1)
結果
WA  
実行時間 -
コード長 1,453 bytes
コンパイル時間 886 ms
コンパイル使用メモリ 108,096 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 21:07:05
合計ジャッジ時間 6,184 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 1 ms
6,940 KB
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 AC 1 ms
6,940 KB
testcase_06 AC 1 ms
6,940 KB
testcase_07 AC 1 ms
6,944 KB
testcase_08 AC 9 ms
6,940 KB
testcase_09 AC 24 ms
6,944 KB
testcase_10 AC 32 ms
6,940 KB
testcase_11 AC 53 ms
6,940 KB
testcase_12 AC 65 ms
6,940 KB
testcase_13 AC 67 ms
6,940 KB
testcase_14 AC 75 ms
6,940 KB
testcase_15 AC 90 ms
6,944 KB
testcase_16 AC 92 ms
6,944 KB
testcase_17 AC 93 ms
6,944 KB
testcase_18 AC 102 ms
6,944 KB
testcase_19 AC 99 ms
6,944 KB
testcase_20 AC 101 ms
6,940 KB
testcase_21 AC 103 ms
6,944 KB
testcase_22 AC 108 ms
6,940 KB
testcase_23 AC 109 ms
6,944 KB
testcase_24 AC 111 ms
6,940 KB
testcase_25 AC 110 ms
6,944 KB
testcase_26 AC 112 ms
6,940 KB
testcase_27 AC 185 ms
6,944 KB
testcase_28 AC 217 ms
6,940 KB
testcase_29 AC 111 ms
6,940 KB
testcase_30 AC 111 ms
6,940 KB
testcase_31 AC 216 ms
6,944 KB
testcase_32 AC 185 ms
6,944 KB
testcase_33 AC 217 ms
6,944 KB
testcase_34 AC 217 ms
6,940 KB
testcase_35 AC 164 ms
6,940 KB
testcase_36 AC 218 ms
6,940 KB
testcase_37 AC 217 ms
6,940 KB
testcase_38 AC 138 ms
6,944 KB
testcase_39 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.conv, std.range, std.stdio, std.string;

void main()
{
  auto n = readln.chomp.to!long;

  auto calcP(long n)
  {
    auto fn = factors(n).assumeSorted, nc = n.ncbrt;
    auto p1 = fn.lowerBound(nc);
    auto p2 = fn.equalRange(nc);
    auto p3 = fn.upperBound(nc);

    long[] r;
    if (!p1.empty) r ~= p1.back;
    if (!p2.empty) r ~= p2.front;
    if (!p3.empty) r ~= p3.front;

    return r;
  }

  auto calcQ(long n)
  {
    auto fn = factors(n).assumeSorted, nc = n.nsqrt;
    auto q1 = fn.lowerBound(nc);
    auto q2 = fn.equalRange(nc);
    auto q3 = fn.upperBound(nc);

    long[] r;
    if (!q1.empty) r ~= q1.back;
    if (!q2.empty) r ~= q2.front;
    if (!q3.empty) r ~= q3.front;

    return r;
  }

  auto r = long.max;
  foreach (p; calcP(n))
    foreach (q; calcQ(n/p))
      r = min(r, p + q + n/p/q - 3);

  writeln(r, " ", n-1);
}

auto factors(long n)
{
  auto r = [1L, n];

  foreach (i; 2..n.nsqrt+1)
    if (n % i == 0) r ~= [i, n/i];

  r.sort();
  return r;
}

pure T nsqrt(T)(T n)
{
  import std.algorithm, std.conv, std.range, core.bitop;
  if (n <= 1) return n;
  T m = 1 << (n.bsr / 2 + 1);
  return iota(1, m).map!"a * a".assumeSorted!"a <= b".lowerBound(n).length.to!T;
}

pure T ncbrt(T)(T n)
{
  import std.algorithm, std.conv, std.range, core.bitop;
  if (n <= 1) return n;
  T m = 1 << (n.bsr / 3 + 1);
  return iota(1, m).map!"a * a * a".assumeSorted!"a <= b".lowerBound(n).length.to!T;
}
0