結果

問題 No.109 N! mod M
ユーザー te-shte-sh
提出日時 2017-04-25 12:23:46
言語 D
(dmd 2.106.1)
結果
AC  
実行時間 1,321 ms / 5,000 ms
コード長 1,195 bytes
コンパイル時間 593 ms
コンパイル使用メモリ 105,088 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 18:52:07
合計ジャッジ時間 2,520 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 69 ms
6,944 KB
testcase_02 AC 42 ms
6,940 KB
testcase_03 AC 1 ms
6,944 KB
testcase_04 AC 17 ms
6,940 KB
testcase_05 AC 1,321 ms
6,940 KB
testcase_06 AC 11 ms
6,944 KB
testcase_07 AC 5 ms
6,944 KB
testcase_08 AC 1 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

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

void main()
{
  auto t = readln.chomp.to!size_t;

  foreach (_; 0..t) {
    auto rd = readln.split.to!(long[]), n = rd[0], m = rd[1];
    if (n >= m) {
      writeln(0);
    } else if (n == 0) {
      writeln(1 % m);
    } else if (n <= 10 ^^ 5) {
      auto r = 1L;
      foreach (i; 2..n+1)
        r = r * i % m;
      writeln(r);
    } else {
      auto minF = minFactor(m);
      if (minF == -1) {
        auto r = m - 1;
        foreach (i; n+1..m)
          r = (r * modInv(i, m) % m + m) % m;
        writeln(r);
      } else {
        if (n >= m / minF) {
          writeln(0);
        } else {
          auto r = 1L;
          foreach (i; 2..n+1)
            r = r * i % m;
          writeln(r);
        }
      }
    }
  }
}

pure T minFactor(T)(T x)
{
  auto m = x.to!real.sqrt.to!T;
  foreach (f; 2..m+1)
    if (x % f == 0) return f;
  return -1;
}

pure T modInv(T)(T x, T m)
{
  T a, b;
  exEuclid(x, m, a, b);
  return a;
}

pure T exEuclid(T)(T a, T b, ref T x, ref T y)
{
  auto g = a;
  x = 1;
  y = 0;
  if (b != 0) {
    g = exEuclid(b, a % b, y, x);
    y -= a / b * x;
  }
  return g;
}
0