結果

問題 No.28 末尾最適化
コンテスト
ユーザー ayaoni
提出日時 2020-12-04 14:28:54
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 998 ms / 5,000 ms
コード長 1,338 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 345 ms
コンパイル使用メモリ 82,208 KB
実行使用メモリ 81,860 KB
最終ジャッジ日時 2024-09-15 00:46:52
合計ジャッジ時間 1,718 ms
ジャッジサーバーID
（参考情報） 		judge6 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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)
0