結果

問題 No.376 立方体のN等分 (2)
コンテスト
ユーザー te-sh
提出日時 2017-07-24 14:38:34
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 489 ms / 5,000 ms
コード長 977 bytes
コンパイル時間 894 ms
コンパイル使用メモリ 105,988 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 21:07:30
合計ジャッジ時間 6,940 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 38
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

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

  auto r = long.max, fn = factors(n);
  foreach (p; fn) {
    if (n / p < p * p) break;
    auto m = n / p;
    foreach (qb; fn) {
      if (qb % p != 0) continue;
      auto q = qb / p;
      if (m / q < q) break;
      r = min(r, p + q + m/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