結果
問題 | No.28 末尾最適化 |
ユーザー | ayaoni |
提出日時 | 2020-12-04 14:28:54 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 998 ms / 5,000 ms |
コード長 | 1,338 bytes |
コンパイル時間 | 345 ms |
コンパイル使用メモリ | 82,208 KB |
実行使用メモリ | 81,860 KB |
最終ジャッジ日時 | 2024-09-15 00:46:52 |
合計ジャッジ時間 | 1,718 ms |
ジャッジサーバーID (参考情報) |
judge6 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 44 ms
61,508 KB |
testcase_01 | AC | 998 ms
81,860 KB |
ソースコード
import sys sys.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()) def prime_factorization(n): # nの素因数分解 res = {} for i in range(2,int(n**.5)+1): if n % i == 0: r = 0 # nがiで何回割り切れるか while n % i == 0: n //= i r += 1 res[i] = r if n != 1: res[n] = 1 return res Q = I() mod = 10**8+9 for _ in range(Q): s,N,K,B = MI() X = [s] for _ in range(N): x = X[-1] X.append(1+(x**2+x*12345) % mod) factorization = prime_factorization(B) prime_factor = list(factorization.keys()) ans = 10**18 for p in prime_factor: A = [] for i in range(N+1): x = X[i] r = 0 while x % p == 0: x //= p r += 1 A.append(r) A.sort() ans = min(ans,sum(A[:K])//factorization[p]) print(ans)