結果

問題 No.660 家を通り過ぎないランダムウォーク問題
ユーザー 0x19f0x19f
提出日時 2018-03-03 12:57:22
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 89 ms / 2,000 ms
コード長 1,205 bytes
コンパイル時間 1,544 ms
コンパイル使用メモリ 148,188 KB
実行使用メモリ 12,392 KB
最終ジャッジ日時 2023-09-06 05:04:58
合計ジャッジ時間 3,892 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,388 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 2 ms
4,380 KB
testcase_07 AC 2 ms
4,384 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 2 ms
4,380 KB
testcase_10 AC 1 ms
4,388 KB
testcase_11 AC 2 ms
4,384 KB
testcase_12 AC 2 ms
4,380 KB
testcase_13 AC 2 ms
4,380 KB
testcase_14 AC 2 ms
4,380 KB
testcase_15 AC 2 ms
4,384 KB
testcase_16 AC 2 ms
4,380 KB
testcase_17 AC 1 ms
4,380 KB
testcase_18 AC 1 ms
4,384 KB
testcase_19 AC 2 ms
4,380 KB
testcase_20 AC 1 ms
4,384 KB
testcase_21 AC 2 ms
4,380 KB
testcase_22 AC 1 ms
4,380 KB
testcase_23 AC 2 ms
4,380 KB
testcase_24 AC 1 ms
4,380 KB
testcase_25 AC 2 ms
4,380 KB
testcase_26 AC 2 ms
4,380 KB
testcase_27 AC 2 ms
4,380 KB
testcase_28 AC 2 ms
4,384 KB
testcase_29 AC 4 ms
4,384 KB
testcase_30 AC 4 ms
4,384 KB
testcase_31 AC 3 ms
4,384 KB
testcase_32 AC 5 ms
4,380 KB
testcase_33 AC 5 ms
4,380 KB
testcase_34 AC 6 ms
4,380 KB
testcase_35 AC 16 ms
4,384 KB
testcase_36 AC 16 ms
4,596 KB
testcase_37 AC 16 ms
4,604 KB
testcase_38 AC 39 ms
6,964 KB
testcase_39 AC 41 ms
7,184 KB
testcase_40 AC 42 ms
7,240 KB
testcase_41 AC 69 ms
10,032 KB
testcase_42 AC 69 ms
10,192 KB
testcase_43 AC 80 ms
11,400 KB
testcase_44 AC 89 ms
12,392 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define REP(i, a, n) for(ll i = ((ll) a); i < ((ll) n); i++)
#define MOD 1000000007LL
using namespace std;
typedef long long ll;

class ModuloCombination {
private:
  vector<ll> fact, inv_fact;

  static ll modulo_power(ll a, ll n) {
    if(n == 0) return 1;
    if(n % 2 == 0) return modulo_power((a * a) % MOD, n / 2);
    return (a * modulo_power(a, n - 1)) % MOD;
  }

public:
  ModuloCombination(ll N): fact(N + 1), inv_fact(N + 1) {
    fact[0] = 1;
    REP(i, 1, N + 1) fact[i] = (i * fact[i - 1]) % MOD;
    REP(i, 0, N + 1) inv_fact[i] = modulo_power(fact[i], MOD - 2);
  }

  ll operator()(ll n, ll r) {
    return fact[n] * inv_fact[r] % MOD * inv_fact[n - r] % MOD;
  }
};

int main(void) {
  ll N;
  cin >> N;

  ModuloCombination comb(N * 2);

  vector<ll> catalan(N + 1);
  catalan[0] = 1;
  REP(i, 1, N + 1) catalan[i] = (comb(2 * i, i) - comb(2 * i, i - 1) + MOD) % MOD;

  vector<ll> sum(N + 2, 0);
  REP(i, 0, N + 1) sum[i + 1] = (sum[i] + catalan[i]) % MOD;

  ll ans = 0;
  REP(i, 0, N / 2 + 1) {
    ll path = comb(N + i * 2, N + i);
    ans = (ans + path) % MOD;
    ans = (ans - (path * sum[N / 2 - i] * 2) % MOD + MOD) % MOD;
  }
  cout << ans << endl;
}
0