結果
問題 | No.1073 無限すごろく |
ユーザー | roaris |
提出日時 | 2020-06-05 23:04:07 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 103 ms / 2,000 ms |
コード長 | 1,201 bytes |
コンパイル時間 | 292 ms |
コンパイル使用メモリ | 86,860 KB |
実行使用メモリ | 77,216 KB |
最終ジャッジ日時 | 2023-08-22 17:09:28 |
合計ジャッジ時間 | 4,938 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 92 ms
71,396 KB |
testcase_01 | AC | 91 ms
71,252 KB |
testcase_02 | AC | 103 ms
76,988 KB |
testcase_03 | AC | 91 ms
71,396 KB |
testcase_04 | AC | 91 ms
71,180 KB |
testcase_05 | AC | 90 ms
71,476 KB |
testcase_06 | AC | 89 ms
71,032 KB |
testcase_07 | AC | 90 ms
71,252 KB |
testcase_08 | AC | 91 ms
71,488 KB |
testcase_09 | AC | 91 ms
71,312 KB |
testcase_10 | AC | 90 ms
71,332 KB |
testcase_11 | AC | 89 ms
71,392 KB |
testcase_12 | AC | 91 ms
71,488 KB |
testcase_13 | AC | 99 ms
77,208 KB |
testcase_14 | AC | 97 ms
77,132 KB |
testcase_15 | AC | 99 ms
77,184 KB |
testcase_16 | AC | 101 ms
77,104 KB |
testcase_17 | AC | 100 ms
77,148 KB |
testcase_18 | AC | 98 ms
77,160 KB |
testcase_19 | AC | 100 ms
76,996 KB |
testcase_20 | AC | 99 ms
76,976 KB |
testcase_21 | AC | 100 ms
77,000 KB |
testcase_22 | AC | 98 ms
77,136 KB |
testcase_23 | AC | 99 ms
77,132 KB |
testcase_24 | AC | 100 ms
77,156 KB |
testcase_25 | AC | 99 ms
77,124 KB |
testcase_26 | AC | 101 ms
77,012 KB |
testcase_27 | AC | 99 ms
77,176 KB |
testcase_28 | AC | 100 ms
77,120 KB |
testcase_29 | AC | 100 ms
77,136 KB |
testcase_30 | AC | 99 ms
77,016 KB |
testcase_31 | AC | 99 ms
77,216 KB |
testcase_32 | AC | 99 ms
77,212 KB |
ソースコード
import sys input = sys.stdin.readline from collections import * def f(x, y): return x*pow(y, MOD-2, MOD)%MOD def mat_mul(A, B): C = [[0]*len(B[0]) for _ in range(len(A))] for i in range(len(A)): for j in range(len(B[0])): for k in range(len(A[0])): C[i][j] += A[i][k]*B[k][j] C[i][j] %= MOD return C def mat_pow(A, n): res = [[0]*len(A) for _ in range(len(A))] for i in range(len(A)): res[i][i] = 1 while n>0: if n&1: res = mat_mul(res, A) n >>= 1 A = mat_mul(A, A) return res N = int(input()) MOD = 10**9+7 x = [[-1] for _ in range(6)] x[0][0] = f(1, 6) x[1][0] = f(7, 6**2) x[2][0] = f(7**2, 6**3) x[3][0] = f(7**3, 6**4) x[4][0] = f(7**4, 6**5) x[5][0] = f(7**5, 6**6) A = [[0]*6 for _ in range(6)] for i in range(6): A[0][i] = f(1, 6) for i in range(1, 6): for j in range(6): for k in range(i): A[i][j] += f(1, 6)*A[k][j] A[i][j] %= MOD if j>=i: A[i][j] += f(1, 6) A[i][j] %= MOD c = (N-1)//6 res = mat_mul(mat_pow(A, c), x) print(res[(N-1)%6][0])