結果
| 問題 | 
                            No.93 ペガサス
                             | 
                    
| コンテスト | |
| ユーザー | 
                             maine_honzuki
                         | 
                    
| 提出日時 | 2020-05-12 23:33:29 | 
| 言語 | C++14  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 61 ms / 5,000 ms | 
| コード長 | 1,487 bytes | 
| コンパイル時間 | 1,859 ms | 
| コンパイル使用メモリ | 176,200 KB | 
| 実行使用メモリ | 5,376 KB | 
| 最終ジャッジ日時 | 2024-09-14 02:58:19 | 
| 合計ジャッジ時間 | 2,958 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge2 / judge3 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 4 | 
| other | AC * 16 | 
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const ll mod = 1e9 + 7;
int main() {
    int N;
    cin >> N;
    vector<vector<ll>> dp(1, vector<ll>{1, 0});
    for (ll i = 1; i < N; i++) {
        ll n = dp.size();
        vector<vector<ll>> nxt(n + 1, vector<ll>(2, 0));
        if (i != N / 2) {
            for (ll a = 0; a <= n; a++) {
                if (a < n)
                    nxt[a][1] += dp[a][1];
                if (a > 0) {
                    nxt[a][1] += 2 * dp[a - 1][0];
                    nxt[a][1] += dp[a - 1][1];
                }
                nxt[a][1] %= mod;
                if (a + 1 < n) {
                    nxt[a][0] += (a + 1) * dp[a + 1][0];
                    nxt[a][0] += a * dp[a + 1][1];
                }
                if (a < n) {
                    nxt[a][0] += max(i - a - 1, 0ll) * dp[a][0];
                    nxt[a][0] += (i - a) * dp[a][1];
                }
                nxt[a][0] %= mod;
            }
        } else {
            for (ll a = 0; a <= n; a++) {
                if (a + 1 < n) {
                    nxt[a][0] += (a + 1) * dp[a + 1][0];
                    nxt[a][0] += (a + 1) * dp[a + 1][1];
                }
                if (a < n) {
                    nxt[a][0] += (i - a + 1) * dp[a][0];
                    nxt[a][0] += (i - a + 1) * dp[a][1];
                }
                nxt[a][0] %= mod;
            }
        }
        dp = nxt;
    }
    cout << dp[0][0] << endl;
}
            
            
            
        
            
maine_honzuki