結果

問題 No.1810 RGB Biscuits
ユーザー titia
提出日時 2022-01-14 21:47:40
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 122 ms / 2,000 ms
コード長 1,102 bytes
コンパイル時間 176 ms
コンパイル使用メモリ 82,236 KB
実行使用メモリ 77,032 KB
最終ジャッジ日時 2024-11-20 09:18:13
合計ジャッジ時間 2,689 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.readline

A,B=map(int,input().split())
mod=10**9+7


# 行列の計算（numpyを使えないとき,modを使用）
def prod(A,B,k,l,m):# A:k*l,B:l*m
    C=[[None for i in range(m)] for j in range(k)]

    for i in range(k):
        for j in range(m):
            ANS=0
            for pl in range(l):
                ANS=(ANS+A[i][pl]*B[pl][j])%mod

            C[i][j]=ANS

    return C

def plus(A,B,k,l):# a,B:k*l
    C=[[None for i in range(l)] for j in range(k)]

    for i in range(k):
        for j in range(l):
            C[i][j]=(A[i][j]+B[i][j])%mod

    return C

# 漸化式を行列累乗で求める（ダブリング）

POWA=[[[A,1],[B,0]]]

for i in range(65):
    POWA.append(prod(POWA[-1],POWA[-1],2,2,2)) # ベキを求めて


N=int(input())
for tests in range(N):
    t=int(input())

    n=t//2

    X=[[1,1]]
    while n:
        X=prod(X,POWA[n.bit_length()-1],1,2,2) # n乗の場合
        n-=1<<(n.bit_length()-1)

    if t%2==0:
        print(sum(X[0])%mod)
    else:
        X=X[0]
        print((X[0]+X[1]+X[0]*A+X[1]*B)%mod)
0