結果
| 問題 |
No.140 みんなで旅行
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2021-03-18 16:44:04 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 2,488 ms / 5,000 ms |
| コード長 | 1,848 bytes |
| コンパイル時間 | 2,788 ms |
| コンパイル使用メモリ | 193,644 KB |
| 最終ジャッジ日時 | 2025-01-19 18:04:05 |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 19 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
const long long MOD = 1e9 + 7;
long long dp[2][560][560]; // dp[][n][k]: n個グループがあり、k個のグループは車を出せる
int main() {
ios_base::sync_with_stdio(0);
cin.tie(0);
dp[0][0][0] = 1;
int n;
cin >> n;
for (int i = 0; i < n; i++) {
int cur = i & 1;
int nxt = (i + 1) & 1;
for (int j = 0; j <= n; j++)
for (int k = 0; k <= n; k++)
dp[nxt][j][k] = 0;
for (int j = 0; j <= n; j++) {
for (int k = 0; k <= n; k++) {
// 女は既存グループに入る。男は女と同じグループに入る。
dp[nxt][j][k] += k * dp[cur][j][k];
dp[nxt][j][k] %= MOD;
dp[nxt][j][k+1] += (j-k) * dp[cur][j][k];
dp[nxt][j][k+1] %= MOD;
// 女は既存グループに入る。男は女と違う既存グループに入る。
dp[nxt][j][k] += j * (j - 1) * dp[cur][j][k];
dp[nxt][j][k] %= MOD;
// 女は既存のグループに入る。男は新しいグループを作る。
dp[nxt][j+1][k] += j * dp[cur][j][k];
dp[nxt][j+1][k] %= MOD;
// 女は新しいグループを作る。男は女と同じグループに入る。
dp[nxt][j+1][k+1] += dp[cur][j][k];
dp[nxt][j+1][k+1] %= MOD;
// 女は新しいグループを作る。男は既存のグループに入る。
dp[nxt][j+1][k] += j * dp[cur][j][k];
dp[nxt][j+1][k] %= MOD;
// 女は新しいグループを作る。男も新しいグループを作る。
dp[nxt][j+2][k] += dp[cur][j][k];
dp[nxt][j+2][k] %= MOD;
}
}
}
long long ret = 0;
for (int i = 1; i <= n; i++) {
for (int j = i; j <= n; j++) {
ret += dp[n&1][i][j];
ret %= MOD;
}
}
cout << ret << endl;
return 0;
}