結果

問題 No.140 みんなで旅行
ユーザー nebukuro09
提出日時 2018-12-04 09:51:03
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 1,484 ms / 5,000 ms
コード長 1,240 bytes
コンパイル時間 1,040 ms
コンパイル使用メモリ 115,500 KB
実行使用メモリ 12,260 KB
最終ジャッジ日時 2024-06-13 02:00:46
合計ジャッジ時間 9,815 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

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, core.stdc.string;

immutable long MOD = 10^^9 + 7;
long[1150][560][2] dp;

void main() {
    auto N = readln.chomp.to!int;
    //auto dp = new long[][][](2, N+5, 2*N+5);
    memset(dp[0].ptr, 0, dp[0].sizeof);
    dp[0][0][0] = 1;

    //foreach (i; 0..N) {
    for (int i = 0; i < N; ++i) {
        int cur = i & 1;
        int tar = cur ^ 1;
        memset(dp[tar].ptr, 0, dp[tar].sizeof);
        //foreach (j; 0..N+1) if (j < i + 1) foreach (k; 0..2*N+1) {
        for (int j = 0; j < i + 1; ++j) for (int k = 0; j + k <= i * 2; ++k) {
            dp[tar][j][k+2] += dp[cur][j][k];
            dp[tar][j][k+1] += dp[cur][j][k] * 2 * (j+k);
            dp[tar][j][k]   += dp[cur][j][k] * ((j+k) * (j+k) - k);
            dp[tar][j+1][k] += dp[cur][j][k];
            if (k > 0) dp[tar][j+1][k-1] += dp[cur][j][k] * k;
        }
        for (int j = 0; j < i + 2; ++j) for (int k = 0; j + k <= (i + 1) * 2; ++k) {
            dp[tar][j][k] %= MOD;
        }
    }

    long ans = 0;
    foreach (i; 1..N+1) (ans += dp[N % 2][i][0]) %= MOD;
    ans.writeln;
}
0