結果

問題 No.840 ほむほむほむら
コンテスト
ユーザー t33f
提出日時 2019-06-26 23:34:35
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 150 ms / 4,000 ms
コード長 1,619 bytes
コンパイル時間 640 ms
コンパイル使用メモリ 66,036 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-26 06:49:42
合計ジャッジ時間 2,517 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstring>
#include <iostream>
using namespace std;
int n, k, s;
constexpr int M = 998244353;

inline int index(int x, int y, int z) { return x + k * (y + k * z); }

void mmul(long long *X, const long long *A, const long long *B) {
    static long long *tmp = nullptr;
    if (tmp == nullptr) tmp = new long long[s*s];
    for (int i = 0; i < s; i++)
        for (int j = 0; j < s; j++) {
            long long sum = 0;
            for (int k = 0; k < s; k++) sum += A[i*s+k] * B[k*s+j] % M;
            tmp[i*s+j] = sum % M;
        }
    memcpy(X, tmp, sizeof(*X)*s*s);
}

void mpow(long long *A, int k) {
    static long long *tmp = nullptr, *P = nullptr;
    if (tmp == nullptr) tmp = new long long[s*s];
    if (P == nullptr) P = new long long[s*s];
    for (int i = 0; i < s; i++)
        for (int j = 0; j < s; j++)
            tmp[s*i + j] = (i==j), P[s*i+j] = A[s*i+j];
    while (k > 0) {
        if (k % 2) mmul(tmp, tmp, P);
        mmul(P, P, P);
        k /= 2;
    }
    memcpy(A, tmp, sizeof(*A)*s*s);
}

int main() {
    cin >> n >> k;
    s = k*k*k;
    long long A[s][s];
    memset(A, 0, sizeof A);
    for (int x = 0; x < k; x++)
        for (int y = 0; y < k; y++)
            for (int z = 0; z < k; z++) {
                int i = index(x, y, z);
                A[index((x+1)%k, y, z)][i]++;
                A[index(x, (y+x)%k, z)][i]++;
                A[index(x, y, (z+y)%k)][i]++;
            }

    mpow(&A[0][0], n);
    long long ans = 0;
    for (int x = 0; x < k; x++)
        for (int y = 0; y < k; y++)
            ans += A[index(x, y, 0)][0] % M;

    cout << ans%M << endl;
}
0