結果

問題 No.93 ペガサス
コンテスト
ユーザー maine_honzuki
提出日時 2020-05-12 23:33:29
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 61 ms / 5,000 ms
コード長 1,487 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 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
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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;
}
0