結果
問題 | No.336 門松列列 |
ユーザー |
![]() |
提出日時 | 2017-04-22 13:47:48 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.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 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 8 |
ソースコード
#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; }