結果
問題 | No.336 門松列列 |
ユーザー | 👑 はまやんはまやん |
提出日時 | 2017-04-22 13:47:48 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 164 ms / 2,000 ms |
コード長 | 1,431 bytes |
コンパイル時間 | 1,815 ms |
コンパイル使用メモリ | 165,060 KB |
実行使用メモリ | 5,152 KB |
最終ジャッジ日時 | 2023-09-02 18:45:27 |
合計ジャッジ時間 | 4,360 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge11 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 156 ms
4,904 KB |
testcase_01 | AC | 121 ms
5,104 KB |
testcase_02 | AC | 120 ms
4,872 KB |
testcase_03 | AC | 121 ms
4,876 KB |
testcase_04 | AC | 164 ms
4,896 KB |
testcase_05 | AC | 120 ms
4,892 KB |
testcase_06 | AC | 125 ms
4,948 KB |
testcase_07 | AC | 136 ms
4,908 KB |
testcase_08 | AC | 148 ms
4,880 KB |
testcase_09 | AC | 164 ms
4,984 KB |
testcase_10 | AC | 164 ms
5,152 KB |
testcase_11 | AC | 164 ms
4,952 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; #define rep(i,a,b) for(int i=a;i<b;i++) int mod = 1000000007; int add(int x, int y) { return (x += y) >= mod ? x - mod : x; } template<class... T> int add(int x, T... y) { return add(x, add(y...)); } int mul(int x, int y) { return 1LL * x * y % mod; } template<class... T> int mul(int x, T... y) { return mul(x, mul(y...)); } int sub(int x, int y) { return add(x, mod - y); } int modpow(int a, long long b) { int ret = 1; while (b > 0) { if (b & 1) ret = 1LL * ret * a % mod; a = 1LL * a * a % mod; b >>= 1; } return ret; } int modinv(int a) { return modpow(a, mod - 2); } int fac[201010], ifac[201010]; void initfac() { fac[0] = ifac[0] = 1; rep(i, 1, 201010) fac[i] = 1LL * i * fac[i - 1] % mod; rep(i, 1, 201010) ifac[i] = modinv(fac[i]); } int aCb(int a, int b) { if (b < 0 || a < b) return 0; return (1LL * fac[a] * ifac[a - b] % mod) * ifac[b] % mod; } //----------------------------------------------------------------------------------- int N; int dp[2020]; //----------------------------------------------------------------------------------- int main() { initfac(); cin >> N; if (N <= 2) { printf("0\n"); return 0; } dp[0] = dp[1] = dp[2] = 1; dp[3] = 2; rep(i, 3, N) for (int j = 0; j <= i;j+=2) { dp[i + 1] = add(dp[i + 1], mul(aCb(i, j), dp[j], dp[i - j])); } cout << mul(dp[N],2) << endl; }