結果

問題 No.1105 Many Triplets
コンテスト
ユーザー lam6er
提出日時 2025-03-20 21:06:20
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 43 ms / 2,000 ms
コード長 1,264 bytes
コンパイル時間 164 ms
コンパイル使用メモリ 81,788 KB
実行使用メモリ 62,120 KB
最終ジャッジ日時 2025-03-20 21:06:27
合計ジャッジ時間 2,161 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

MOD = 10**9 + 7

def multiply_matrix(A, B):
    res = [[0] * 3 for _ in range(3)]
    for i in range(3):
        for j in range(3):
            for k in range(3):
                res[i][j] += A[i][k] * B[k][j]
                res[i][j] %= MOD
    return res

def matrix_power(mat, power):
    result = [[1 if i == j else 0 for j in range(3)] for i in range(3)]
    while power > 0:
        if power % 2 == 1:
            result = multiply_matrix(result, mat)
        mat = multiply_matrix(mat, mat)
        power = power // 2
    return result

def multiply_vector(mat, vec):
    res = [0, 0, 0]
    for i in range(3):
        res[i] = (mat[i][0] * vec[0] + mat[i][1] * vec[1] + mat[i][2] * vec[2]) % MOD
    return res

def main():
    import sys
    input = sys.stdin.read().split()
    N = int(input[0])
    a = int(input[1]) % MOD
    b = int(input[2]) % MOD
    c = int(input[3]) % MOD
    vec = [a, b, c]
    
    if N == 1:
        print(a, b, c)
        return
    
    M = [
        [1, (MOD - 1) % MOD, 0],
        [0, 1, (MOD - 1) % MOD],
        [(MOD - 1) % MOD, 0, 1]
    ]
    
    power = N - 1
    M_exp = matrix_power(M, power)
    final = multiply_vector(M_exp, vec)
    print(final[0], final[1], final[2])

if __name__ == '__main__':
    main()
0