結果

問題 No.28 末尾最適化
ユーザー Navier_BoltzmannNavier_Boltzmann
提出日時 2023-06-08 21:12:12
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,038 ms / 5,000 ms
コード長 812 bytes
コンパイル時間 372 ms
コンパイル使用メモリ 82,020 KB
実行使用メモリ 84,576 KB
最終ジャッジ日時 2024-12-31 05:34:22
合計ジャッジ時間 2,113 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import *
from itertools import *
from functools import *
from heapq import *
import sys,math
input = sys.stdin.readline
mod = 10**8 + 9
P = [2,3,5,7,11,13,17,19,23,29,31]
def answer():
    seed,N,K,B = map(int,input().split())
    X = [seed]
    for _ in range(N):
        seed = 1 + (seed**2 + seed*12345)%mod
        X.append(seed)
    
    Z = []
    for p in P:
        while B%p==0:
            B //= p
            Z.append(p)
    C = Counter(Z)
    ans = 1<<60
    for k,v in C.items():
        Y = []
        for x in X:
            y = x
            tmp = 0
            while y%k==0:
                tmp += 1
                y //= k
            Y.append(tmp)
        Y.sort()
        ans = min(ans,sum(Y[:K])//v)

    print(ans)
    
for _ in range(int(input())):
    answer()
0