結果
問題 | No.1704 Many Bus Stops (easy) |
ユーザー |
![]() |
提出日時 | 2021-10-10 12:56:23 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 829 ms / 2,000 ms |
コード長 | 1,041 bytes |
コンパイル時間 | 408 ms |
コンパイル使用メモリ | 82,296 KB |
実行使用メモリ | 77,704 KB |
最終ジャッジ日時 | 2024-09-14 07:21:11 |
合計ジャッジ時間 | 26,353 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 41 |
ソースコード
import sys#input = sys.stdin.readlineinput = sys.stdin.buffer.readlinedef mat_dot(A, B, MOD=pow(10,9)+7):#A:n*m行列, B:m*l行列n, m, l = len(A), len(A[0]), len(B[0])X = [[0]*l for i in range(n)]for i in range(n):for j in range(l):temp = 0for k in range(m):temp = (temp + A[i][k] * B[k][j]) % MODX[i][j] = tempreturn Xdef mat_pow(A, x, MOD):#A^xn = len(A)X = [[0] * n for i in range(n)]for i in range(n):X[i][i] = 1for i in range(x.bit_length()):if (x>>i)&1:X = mat_dot(X, A, MOD)A = mat_dot(A, A, MOD)return Xdef main():T = int(input()); MOD = pow(10,9) + 7inv3 = pow(3,MOD-2,MOD)X = [[inv3,0,0,1,0,0],[0,inv3,0,0,1,0],[0,0,inv3,0,0,1],[0,inv3,inv3,0,0,0],[inv3,0,inv3,0,0,0],[inv3,inv3,0,0,0,0]]for _ in range(T):N = int(input())v = [[1],[0],[0],[0],[0],[0]]XN = mat_pow(X,N,MOD)g = mat_dot(XN,v,MOD)#print(XN)print(g[0][0])if __name__ == '__main__':main()