結果
問題 | No.1632 Sorting Integers (GCD of M) |
ユーザー |
![]() |
提出日時 | 2021-07-30 21:05:44 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 58 ms / 2,000 ms |
コード長 | 1,534 bytes |
コンパイル時間 | 540 ms |
コンパイル使用メモリ | 82,404 KB |
実行使用メモリ | 69,228 KB |
最終ジャッジ日時 | 2024-09-21 17:30:20 |
合計ジャッジ時間 | 4,568 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 59 |
ソースコード
import sysfrom math import gcdsys.setrecursionlimit(10**7)def I(): return int(sys.stdin.readline().rstrip())def MI(): return map(int,sys.stdin.readline().rstrip().split())def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))def LI2(): return list(map(int,sys.stdin.readline().rstrip()))def S(): return sys.stdin.readline().rstrip()def LS(): return list(sys.stdin.readline().rstrip().split())def LS2(): return list(sys.stdin.readline().rstrip())N = I()C = [0]+LI()mod = 10**9+7def matrix_multiplication_mod(A,B,p):l = len(A)m = len(B)n = len(B[0])C = [[0]*n for _ in range(l)]for i in range(l):for j in range(n):C[i][j] = sum((A[i][k]*B[k][j]) % p for k in range(m)) % preturn Cdef matrix_power_mod(A,n,p):l = len(A)C = [[0] * l for _ in range(l)]for i in range(l):C[i][i] = 1while n > 0:if n % 2 == 1:C = matrix_multiplication_mod(C,A,p)A = matrix_multiplication_mod(A,A,p)n >>= 1return CX = [[10,1],[0,1]]for i in range(1,10):if C[i] == N:ans = i*matrix_power_mod(X,N,mod)[0][1]ans %= modprint(ans)exit()a = 0for i in range(1,10):if not C[i]:continuefor j in range(i+1,10):if not C[j]:continuea = gcd(a,9*(j-i))s = 0r = 0for i in range(1,10):c = C[i]b = matrix_power_mod(X,c,a)[0][1]s += i*b*pow(10,r,a)s %= ar += cans = gcd(a,s)print(ans)