結果

問題 No.93 ペガサス
コンテスト
ユーザー szkhts
提出日時 2026-03-13 11:40:11
言語 Python3
(3.14.3 + numpy 2.4.2 + scipy 1.17.0)
コンパイル:
python3 -mpy_compile _filename_
実行:
python3 _filename_
結果
TLE  
実行時間 -
コード長 1,906 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 460 ms
コンパイル使用メモリ 20,828 KB
実行使用メモリ 366,684 KB
最終ジャッジ日時 2026-03-13 11:40:44
合計ジャッジ時間 27,810 ms
ジャッジサーバーID
（参考情報） 		judge2_1 / judge3_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 15 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

N=int(input())
mod=1000000007
dp=[[[[0 for i in range(2)] for j in range(2)] for k in range(N+1)] for l in range(N+2)]
dp[1][0][0][0] = 1
dp[2][0][0][0] = 2

for i in range(2, N+1):
    for j in range(N+1):
        if dp[i][j][0][0] > 0:
            dp[i+1][j][0][0] += (i-1-j) * dp[i][j][0][0]
            dp[i+1][j][0][0] %= mod
            dp[i+1][j][0][1] += 2 * dp[i][j][0][0]
            dp[i+1][j][0][1] %= mod
            if j-1 >= 0:
                dp[i+1][j-1][0][0] += j * dp[i][j][0][0]
                dp[i+1][j-1][0][0] %= mod

        if dp[i][j][0][1] > 0:
            dp[i+1][j][1][0] += (i-2-j) * dp[i][j][0][1]
            dp[i+1][j][1][0] %= mod
            dp[i+1][j][1][1] += 2 * dp[i][j][0][1]
            dp[i+1][j][1][1] %= mod
            dp[i+1][j][0][0] += dp[i][j][0][1]
            dp[i+1][j][0][0] %= mod
            if j-1 >= 0:
                dp[i+1][j-1][1][0] += j * dp[i][j][0][1]
                dp[i+1][j-1][1][0] %= mod

        if dp[i][j][1][0]>0:
            dp[i+1][j+1][0][0] += (i-1-j)*dp[i][j][1][0]
            dp[i+1][j+1][0][0] %= mod
            dp[i+1][j+1][0][1] += dp[i][j][1][0]
            dp[i+1][j+1][0][1] %= mod
            dp[i+1][j][0][1] += dp[i][j][1][0]
            dp[i+1][j][0][1] %= mod
            if j-1 >= 0:
                dp[i+1][j][0][0] += j * dp[i][j][1][0]
                dp[i+1][j][0][0] %= mod

        if dp[i][j][1][1] > 0:
            dp[i+1][j+1][1][0] += (i-2-j) * dp[i][j][1][1]
            dp[i+1][j+1][1][0] %= mod
            dp[i+1][j+1][1][1] += dp[i][j][1][1]
            dp[i+1][j+1][1][1] %= mod
            dp[i+1][j][1][1] += dp[i][j][1][1]
            dp[i+1][j][1][1] %= mod
            dp[i+1][j+1][0][0] += dp[i][j][1][1]
            dp[i+1][j+1][0][0] %= mod
            if j-1 >= 0:
                dp[i+1][j][1][0] += j * dp[i][j][1][1]
                dp[i+1][j][1][0] %= mod

print(dp[N][0][0][0])
0