結果

問題 No.1050 Zero (Maximum)
ユーザー ooaiu
提出日時 2024-12-16 13:38:55
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 13 ms / 2,000 ms
コード長 1,233 bytes
コンパイル時間 2,865 ms
コンパイル使用メモリ 255,016 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-16 13:38:59
合計ジャッジ時間 3,885 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff # 

#ifndef LOCAL
#include <bits/stdc++.h>

using namespace std;

#define debug(...) (void(0))
#else
#include "algo/debug.h"
#endif

#include <atcoder/modint>
using mint = atcoder::modint1000000007;

void solve() {
    int M, K;
    cin >> M >> K;
    vector dp(M, vector<mint>(M));
    for (int i = 0; i < M; i++) {
        for (int j = 0; j < M; j++) {
            dp[(i + j) % M][i] += 1;
            dp[(i * j) % M][i] += 1;
        }
    }
    vector ans(M, vector<mint>(M));
    for (int i = 0; i < M; i++) ans[i][i] = 1;
    using dat = vector<vector<mint>>;
    const auto op = [&](const dat& a, const dat& b) -> dat {
        dat res(M, vector<mint>(M));
        for (int i = 0; i < M; i++) {
            for (int k = 0; k < M; k++) {
                for (int j = 0; j < M; j++) {
                    res[i][j] += a[i][k] * b[k][j];
                }
            }
        }
        return res;
    };
    while (K) {
        if (K & 1) {
            ans = op(ans, dp);
        }
        dp = op(dp, dp);
        K >>= 1;
    }
    cout << ans[0][0].val() << endl;
}

int main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    int tt = 1;
    // std::cin >> tt;
    while (tt--) {
        solve();
    }
}
0