結果

問題 No.28 末尾最適化
コンテスト
ユーザー ayaoni
提出日時 2020-12-04 14:24:18
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,451 bytes
コンパイル時間 526 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 182,912 KB
最終ジャッジ日時 2024-09-15 00:41:44
合計ジャッジ時間 7,126 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 1 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

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())
    
    d = {}
    for p in prime_factor:
        d[p] = []
    for p in prime_factor:
        for i in range(N+1):
            x = X[i]
            r = 0
            while x % p == 0:
                x //= p
                r += 1
            d[p].append((r,i))

    ans = 10**18
    for p in prime_factor:
        d[p].sort()
        ans = min(ans,sum(d[p][i][0] for i in range(K))//factorization[p])
    print(ans)
0