結果

問題 No.181 A↑↑N mod M
ユーザー te-shte-sh
提出日時 2017-05-09 12:06:18
言語 D
(dmd 2.106.1)
結果
AC  
実行時間 2 ms / 5,000 ms
コード長 1,326 bytes
コンパイル時間 961 ms
コンパイル使用メモリ 103,416 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 19:05:04
合計ジャッジ時間 2,093 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

diff # 

import std.algorithm, std.conv, std.range, std.stdio, std.string;
import std.typecons;  // Tuple, Nullable, BigFlags
import std.math;      // math functions

void main()
{
  auto rd = readln.split.to!(long[]);
  writeln(calc(rd[0], rd[1], rd[2]));
}

auto calc(long a, long n, long m)
{
  if (m == 1) return 0L;
  if (n == 0) return 1L;

  auto c = calcCyc(a, m);
  auto t = modTet(a, n - 1, m);

  if (t < 0)
    return modPow(a, m + calc(a, n - 1, c) - m % c, m);
  else
    return modPow(a, t, m);
}

auto calcCyc(long a, long m)
{
  auto bi = new size_t[](m);
  bi[] = size_t.max;
  auto ci = new long[](m*2);
  ci[0] = 1;
  foreach (i; 1..m*2) {
    ci[i] = ci[i-1] * a % m;
    if (bi[ci[i]] != size_t.max)
      return i - bi[ci[i]];
    else
      bi[ci[i]] = i;
  }
  assert(0);
}

auto modPow(long a, long n, long m)
{
  if (m == 1) return 0L;
  if (n == 0) return 1L;

  auto r = 1L, s = a;
  while (n > 0) {
    if (n & 1)
      r = (r * s) % m;
    s = (s * s) % m;
    n >>= 1;
  }

  return r;
}

auto modTet(long a, long n, long m)
{
  long r;
  if (n == 0) {
    return 1L;
  } else if (n == 1) {
    r = a;
  } else if (n == 2) {
    if (a >= 5) return -1L;
    r = a ^^ a;
  } else if (n == 3) {
    if (a >= 3) return -1L;
    r = a ^^ (a ^^ a);
  } else {
    return -1L;
  }
  return r >= m ? -1L : r;
}
0