結果

問題 No.117 組み合わせの数
ユーザー te-shte-sh
提出日時 2017-01-30 14:58:32
言語 D
(dmd 2.106.1)
結果
WA  
実行時間 -
コード長 1,617 bytes
コンパイル時間 7,969 ms
コンパイル使用メモリ 300,192 KB
実行使用メモリ 27,476 KB
最終ジャッジ日時 2024-06-12 06:42:26
合計ジャッジ時間 8,913 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.conv, std.range, std.stdio, std.string;
import std.regex;     // RegEx

const int p = 10 ^^ 9 + 7;

void main()
{
  auto bufFrac = new int[](2_000_000), sf = 1;
  bufFrac[0] = 1;

  auto bufInvFrac = new int[](2_000_000);
  bufInvFrac[] = -1;

  auto frac(int n) {
    if (sf <= n) {
      foreach (i; sf..n+1)
        bufFrac[i] = modMul(bufFrac[i - 1], i, p);
      sf = n + 1;
    }
    return bufFrac[n];
  }

  auto invFrac(int n) {
    if (bufInvFrac[n] == -1)
      bufInvFrac[n] = invMod(frac(n), p);
    return bufInvFrac[n];
  }

  writeln(invFrac(1));

  auto t = readln.chomp.to!size_t;
  foreach (_; t.iota) {
    auto rd = readln.chomp;
    auto m = rd.matchFirst(r"([CPH])\((\d+),(\d+)\)");
    auto n = m[2].to!int, k = m[3].to!int;

    auto r = 0;
    switch (m[1]) {
    case "C":
      if (n >= k) {
        r = frac(n);
        r = modMul(r, invFrac(k), p);
        r = modMul(r, invFrac(n - k), p);
      }
      break;
    case "P":
      if (n >= k) {
        r = frac(n);
        r = modMul(r, invFrac(n - k), p);
      }
      break;
    case "H":
      r = frac(n + k - 1);
      r = modMul(r, invFrac(k), p);
      r = modMul(r, invFrac(n - 1), p);
      break;
    default:
      assert(0);
    }

    writeln(r);
  }
}

int modMul(int a, int b, int mod)
{
  return ((a.to!long * b.to!long) % mod).to!int;
}

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;
}

T invMod(T)(T x, T m)
{
  T a, b;
  exEuclid(x, m, a, b);
  return ((a % m) + m) % m;
}
0