結果

問題 No.767 配られたジャパリまん
ユーザー ikdikd
提出日時 2018-12-25 16:57:11
言語 D
(dmd 2.106.1)
結果
AC  
実行時間 575 ms / 2,000 ms
コード長 1,801 bytes
コンパイル時間 678 ms
コンパイル使用メモリ 105,584 KB
実行使用メモリ 13,568 KB
最終ジャッジ日時 2024-06-13 02:25:33
合計ジャッジ時間 2,832 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,940 KB
testcase_03 AC 1 ms
6,944 KB
testcase_04 AC 3 ms
6,940 KB
testcase_05 AC 70 ms
6,944 KB
testcase_06 AC 1 ms
5,376 KB
testcase_07 AC 1 ms
5,376 KB
testcase_08 AC 1 ms
5,376 KB
testcase_09 AC 13 ms
5,376 KB
testcase_10 AC 39 ms
5,376 KB
testcase_11 AC 27 ms
5,376 KB
testcase_12 AC 18 ms
5,376 KB
testcase_13 AC 35 ms
5,376 KB
testcase_14 AC 28 ms
5,376 KB
testcase_15 AC 33 ms
5,376 KB
testcase_16 AC 24 ms
5,376 KB
testcase_17 AC 11 ms
5,376 KB
testcase_18 AC 13 ms
5,376 KB
testcase_19 AC 32 ms
5,376 KB
testcase_20 AC 575 ms
13,568 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

void main() {
  import std.stdio, std.string, std.conv, std.algorithm;
  import core.bitop : popcnt;

  const mod = 10 ^^ 8 + 7;
  long powmod(long b, long e) {
    if (e == 0)
      return 1;
    else if (e == 1)
      return b;
    else if (e & 1)
      return b * powmod(b, e - 1) % mod;
    else
      return powmod(b * b % mod, e / 2);
  }

  int h, w, k;
  rd(h, w, k);
  struct P {
    int y, x, idx;
  }

  auto fac = new long[](h + w + 1), inv = new long[](h + w + 1);
  fac[0] = inv[0] = 1;
  foreach (i; 1 .. (h + w + 1)) {
    fac[i] = fac[i - 1] * i % mod;
    inv[i] = powmod(fac[i], mod - 2);
  }
  long cmb(int n, int r) {
    if (r < 0 || n < r)
      return 0;
    else
      return fac[n] * inv[r] % mod * inv[n - r] % mod;
  }

  auto abi = new P[](k);
  foreach (i; 0 .. k) {
    int y, x;
    rd(y, x);
    abi[i] = P(y, x, i);
  }
  abi.sort!((p, q) => (p.y == q.y ? p.x < q.x : p.y < q.y));
  auto dp = new long[](1 << k);
  fill(dp, 1);
  foreach (bits; 0 .. (1 << k)) {
    int cur_x = 0, cur_y = 0;
    foreach (p; abi) {
      if (bits & (1 << p.idx)) {
        int dy = p.y - cur_y, dx = p.x - cur_x;
        dp[bits] = dp[bits] * cmb(dy + dx, dy) % mod;
        cur_y = p.y;
        cur_x = p.x;
      }
    }
    dp[bits] = dp[bits] * cmb((h - cur_y) + (w - cur_x), h - cur_y) % mod;
  }
  foreach (i; 0 .. k) {
    foreach (bits; 0 .. (1 << k)) {
      if ((bits & (1 << i)) == 0) {
        dp[bits ^ (1 << i)] = (dp[bits ^ (1 << i)] + (mod - dp[bits])) % mod;
      }
    }
  }
  foreach (bits, x; dp) {
    if (popcnt(bits) & 1)
      writeln(mod - x);
    else
      writeln(x);
  }
}

void rd(T...)(ref T x) {
  import std.stdio, std.string, std.conv;

  auto l = readln.split;
  assert(l.length == x.length);
  foreach (i, ref e; x)
    e = l[i].to!(typeof(e));
}
0