結果

問題 No.1513 simple 門松列 problem
ユーザー ktr216
提出日時 2021-05-21 22:22:39
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,600 bytes
コンパイル時間 1,719 ms
コンパイル使用メモリ 174,696 KB
実行使用メモリ 137,892 KB
最終ジャッジ日時 2024-10-10 09:09:28
合計ジャッジ時間 6,451 ms
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2 TLE * 1
other -- * 18
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
using namespace std;

typedef long long ll;

int main() {
    int N, K;
    scanf("%d %d", &N, &K);
    // dp[i][j][k] i番目がj, i+1番目がkの門松列の数
    vector<vector<vector<ll>>> dp(N - 1, vector<vector<ll>>(K, vector<ll>(K)));
    vector<vector<vector<ll>>> dp2(N - 1, vector<vector<ll>>(K, vector<ll>(K)));
    rep(i, K) rep(j, K) {
        if (i != j) dp[0][i][j] = 1;
        if (i != j) dp2[0][i][j] = i + j;
    }
    ll MOD = 998244353;
    rep(i, N - 2) {
        rep(j, K) {
            rep(k, K) {
                if (j < k) {
                    rep(l, k) {
                        if (l == j) continue;
                        dp[i + 1][k][l] += dp[i][j][k];
                        dp[i + 1][k][l] %= MOD;
                        dp2[i + 1][k][l] += dp2[i][j][k] + l * dp[i][j][k];
                        dp2[i + 1][k][l] %= MOD;
                    }
                }
                if (j > k) {
                    for (int l = k + 1; l < K; l++) {
                        if (l == j) continue;
                        dp[i + 1][k][l] += dp[i][j][k];
                        dp[i + 1][k][l] %= MOD;
                        dp2[i + 1][k][l] += dp2[i][j][k] + l * dp[i][j][k];
                        dp2[i + 1][k][l] %= MOD;
                    }
                }
            }
        }
    }
    ll ans1 = 0, ans2 = 0;
    rep(j, K) rep(k, K) {
        ans1 = (ans1 + dp[N - 2][j][k]) % MOD;
        ans2 = (ans2 + dp2[N - 2][j][k]) % MOD;
    }
    printf("%lld %lld\n", ans1, ans2);
}
0