結果

問題 No.93 ペガサス
コンテスト
ユーザー vjudge1
提出日時 2025-09-06 18:32:07
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 1,253 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,702 ms
コンパイル使用メモリ 161,560 KB
実行使用メモリ 22,784 KB
最終ジャッジ日時 2025-09-06 18:32:10
合計ジャッジ時間 2,690 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 WA * 2
other AC * 2 WA * 14
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>

using namespace std;

const int MAXN = 1e3 + 10, P = 1e9 + 7;

int n;
long long dp[MAXN][MAXN][2][2];

int main(){
  cin >> n;
  dp[1][0][0][0] = 1, dp[2][0][0][0] = 2;
  for(int i = 2; i <= n; i++){
    for(int j = 0; j <= i; j++){
      (dp[i + 1][j + 1][0][1] += dp[i][j][0][0] * 2) %= P;
      (dp[i + 1][j - 1][0][0] += j * dp[i][j][0][0]) %= P;
      (dp[i + 1][j][0][0] += (i + 1 - j - 2) * dp[i][j][0][0]) %= P;

      (dp[i + 1][j + 1][1][1] += dp[i][j][0][1] * 2) %= P;
      (dp[i + 1][j - 1][1][0] += (j - 1) * dp[i][j][0][1]) %= P;
      (dp[i + 1][j - 1][0][0] += dp[i][j][0][1]) %= P;
      (dp[i + 1][j][1][0] += (i + 1 - (j - 1) - 1 - 2) * dp[i][j][0][1]) %= P;

      (dp[i + 1][j + 1][1][1] += dp[i][j][1][1]) %= P;
      (dp[i + 1][j + 1][1][1] += dp[i][j][1][1]) %= P;
      (dp[i + 1][j - 1][1][1] += (j - 1) * dp[i][j][1][1]) %= P;
      (dp[i + 1][j - 1][0][1] += dp[i][j][1][1]) %= P;
      (dp[i + 1][j][1][0] += (i + 1 - (j - 1) - 1 - 2) * dp[i][j][1][1]) %= P;

      (dp[i + 1][j + 1][0][1] += dp[i][j][1][0] * 2) %= P;
      (dp[i + 1][j - 1][0][0] += j * dp[i][j][1][0]) %= P;
      (dp[i + 1][j][1][0] += (i + 1 - j - 2) * dp[i][j][1][0]) %= P;
    }
  }
  cout << dp[n][0][0][0];
  return 0;
}
0