結果
| 問題 | No.140 みんなで旅行 | 
| コンテスト | |
| ユーザー |  kyuna | 
| 提出日時 | 2019-08-19 19:04:35 | 
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 2,581 ms / 5,000 ms | 
| コード長 | 811 bytes | 
| コンパイル時間 | 649 ms | 
| コンパイル使用メモリ | 76,528 KB | 
| 実行使用メモリ | 22,680 KB | 
| 最終ジャッジ日時 | 2024-10-04 20:22:26 | 
| 合計ジャッジ時間 | 14,053 ms | 
| ジャッジサーバーID (参考情報) | judge4 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 19 | 
ソースコード
#include <algorithm>
#include <iostream>
#include <vector>
using namespace std;
using ll = long long;
const int MOD = (int)1e9 + 7;
int main() {
    int n; cin >> n;
    vector<vector<ll>> dp(n * 2 + 1, vector<ll>(n * 2 + 1));
    dp[0][0] = 1;
	for (int t = 0; t < n; t++) {
        vector<vector<ll>> dp2(n * 2 + 1, vector<ll>(n * 2 + 1));
		for (int i = 0; i <= t * 2; i++) {
			for (int j = 0; i + j <= t * 2; j++) {
				(dp2[i][j + 2] += dp[i][j]) %= MOD;
				(dp2[i][j + 1] += dp[i][j] * (i + j) * 2) %= MOD;
				(dp2[i][j] += dp[i][j] * ((i + j) * (i + j) - j)) %= MOD;
				(dp2[i + 1][j] += dp[i][j]) %= MOD;
				if (j) (dp2[i + 1][j - 1] += dp[i][j] * j) %= MOD;
			}
		}
        swap(dp, dp2);
	}
	long long ans = 0;
	for (int i = 0; i <= n * 2; i++) (ans += dp[i][0]) %= MOD;
	cout << ans << endl;
}
            
            
            
        