結果

問題 No.109 N! mod M
ユーザー te-sh
提出日時 2017-04-25 12:21:22
言語 D
(dmd 2.109.1)
結果
WA  
実行時間 -
コード長 1,148 bytes
コンパイル時間 879 ms
コンパイル使用メモリ 105,220 KB
実行使用メモリ 6,940 KB
最終ジャッジ日時 2024-06-12 18:51:59
合計ジャッジ時間 2,808 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 8 WA * 1
権限があれば一括ダウンロードができます

ソースコード

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 <= 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