結果
問題 | No.336 門松列列 |
ユーザー | はまやんはまやん |
提出日時 | 2017-04-22 13:47:48 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 157 ms / 2,000 ms |
コード長 | 1,431 bytes |
コンパイル時間 | 1,756 ms |
コンパイル使用メモリ | 167,976 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-12 01:13:18 |
合計ジャッジ時間 | 4,226 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 147 ms
6,812 KB |
testcase_01 | AC | 113 ms
6,816 KB |
testcase_02 | AC | 108 ms
6,940 KB |
testcase_03 | AC | 109 ms
6,940 KB |
testcase_04 | AC | 152 ms
6,940 KB |
testcase_05 | AC | 110 ms
6,944 KB |
testcase_06 | AC | 113 ms
6,944 KB |
testcase_07 | AC | 128 ms
6,944 KB |
testcase_08 | AC | 139 ms
6,940 KB |
testcase_09 | AC | 157 ms
6,940 KB |
testcase_10 | AC | 151 ms
6,944 KB |
testcase_11 | AC | 152 ms
6,944 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; }