結果

問題 No.1810 RGB Biscuits
ユーザー trineutrontrineutron
提出日時 2022-01-14 21:55:27
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 32 ms / 2,000 ms
コード長 1,574 bytes
コンパイル時間 3,257 ms
コンパイル使用メモリ 256,716 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-04-30 14:27:04
合計ジャッジ時間 4,154 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 2 ms
6,944 KB
testcase_05 AC 30 ms
6,940 KB
testcase_06 AC 32 ms
6,944 KB
testcase_07 AC 32 ms
6,940 KB
testcase_08 AC 32 ms
6,940 KB
testcase_09 AC 32 ms
6,940 KB
testcase_10 AC 21 ms
6,944 KB
testcase_11 AC 4 ms
6,940 KB
testcase_12 AC 2 ms
6,944 KB
testcase_13 AC 2 ms
6,940 KB
testcase_14 AC 7 ms
6,940 KB
testcase_15 AC 2 ms
6,944 KB
testcase_16 AC 3 ms
6,940 KB
testcase_17 AC 17 ms
6,940 KB
testcase_18 AC 17 ms
6,940 KB
testcase_19 AC 2 ms
6,940 KB
testcase_20 AC 2 ms
6,944 KB
testcase_21 AC 6 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/modint>

using namespace std;
using namespace atcoder;

using mint = modint1000000007;
using vec = vector<mint>;
using matrix = vector<vec>;

matrix operator*(const matrix &a, const matrix &b)
{
    const int n = a.size();
    matrix c(n, vec(n));
    for (int i = 0; i < n; i++)
    {
        for (int k = 0; k < n; k++)
        {
            for (int j = 0; j < n; j++)
            {
                c.at(i).at(j) += a.at(i).at(k) * b.at(k).at(j);
            }
        }
    }
    return c;
}

matrix matpow(matrix a, int64_t k)
{
    const int n = a.size();
    matrix res(n, vec(n));
    for (int i = 0; i < n; i++)
    {
        res.at(i).at(i) = 1;
    }
    while (k)
    {
        if (k & 1)
        {
            res = res * a;
        }
        a = a * a;
        k >>= 1;
    }
    return res;
}

int main()
{
    int64_t a, b, n;
    cin >> a >> b >> n;
    for (int i = 0; i < n; i++)
    {
        int64_t t;
        cin >> t;
        matrix mat(3, vec(3)), fst(3, vec(3));
        mat[0][0] = a;
        mat[0][1] = 1;
        mat[1][0] = b;
        mat[2][0] = 1;
        fst[0][0] = 1;
        fst[1][1] = 1;
        fst[0][2] = a;
        fst[1][2] = b;
        matrix newmat = matpow(mat, t / 2);
        if (t % 2)
        {
            newmat = newmat * fst;
        }
        mint ans = 0;
        for (int j = 0; j < 2; j++)
        {
            for (int k = 0; k < 3; k++)
            {
                ans += newmat[j][k];
            }
        }
        cout << ans.val() << endl;
    }
    return 0;
}
0