結果
問題 | No.336 門松列列 |
ユーザー |
|
提出日時 | 2017-11-16 10:11:49 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 1,164 ms / 2,000 ms |
コード長 | 1,284 bytes |
コンパイル時間 | 2,828 ms |
コンパイル使用メモリ | 115,212 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-12 01:13:39 |
合計ジャッジ時間 | 10,983 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 8 |
ソースコード
import std.stdio, std.array, std.string, std.conv, std.algorithm;import std.typecons, std.range, std.random, std.math, std.container;import std.numeric, std.bigint, core.bitop;immutable long MOD = 10^^9+7;void main() {auto N = readln.chomp.to!int;if (N <= 2) {writeln(0);return;}auto dp = new long[][][](N+1, 2, 2); // i, 左端, 右端dp[0][0][1] = 1;dp[0][1][0] = 1;dp[1][0][0] = 1;dp[2][1][0] = 1;dp[2][0][1] = 1;dp[3][0][0] = 2;dp[3][1][1] = 2;foreach (i; 4..N+1) {long comb = 1;foreach (j; 0..i) {int k = i - j - 1;(dp[i][0][0] += dp[j][0][0] * dp[k][0][0] % MOD * comb % MOD) %= MOD;(dp[i][1][0] += dp[j][1][0] * dp[k][0][0] % MOD * comb % MOD) %= MOD;(dp[i][0][1] += dp[j][0][0] * dp[k][0][1] % MOD * comb % MOD) %= MOD;(dp[i][1][1] += dp[j][1][0] * dp[k][0][1] % MOD * comb % MOD) %= MOD;comb = comb * k % MOD * powmod(j + 1, MOD-2, MOD) % MOD;}}writeln((dp[N][0].sum + dp[N][1].sum) % MOD);}long powmod(long a, long x, long m) {long ret = 1;while (x) {if (x % 2) ret = ret * a % m;a = a * a % m;x /= 2;}return ret;}