結果

問題 No.1611 Minimum Multiple with Double Divisors
ユーザー shotoyooshotoyoo
提出日時 2021-07-21 23:00:04
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 2,558 bytes
コンパイル時間 172 ms
コンパイル使用メモリ 82,376 KB
実行使用メモリ 115,056 KB
最終ジャッジ日時 2024-07-17 20:10:04
合計ジャッジ時間 20,591 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 866 ms
114,644 KB
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 AC 648 ms
111,668 KB
testcase_11 AC 635 ms
114,312 KB
testcase_12 AC 625 ms
114,168 KB
testcase_13 AC 653 ms
113,988 KB
testcase_14 AC 694 ms
114,316 KB
testcase_15 AC 722 ms
114,004 KB
testcase_16 AC 671 ms
114,120 KB
testcase_17 AC 725 ms
114,668 KB
testcase_18 AC 651 ms
114,136 KB
testcase_19 AC 213 ms
103,364 KB
testcase_20 AC 196 ms
103,128 KB
testcase_21 AC 213 ms
103,436 KB
testcase_22 AC 199 ms
103,412 KB
testcase_23 AC 202 ms
103,388 KB
testcase_24 AC 208 ms
103,468 KB
testcase_25 AC 205 ms
103,240 KB
testcase_26 AC 211 ms
103,284 KB
testcase_27 AC 210 ms
102,972 KB
testcase_28 AC 149 ms
101,892 KB
testcase_29 AC 150 ms
101,636 KB
testcase_30 AC 145 ms
102,256 KB
testcase_31 AC 153 ms
102,412 KB
testcase_32 AC 159 ms
102,292 KB
testcase_33 AC 155 ms
102,200 KB
testcase_34 AC 155 ms
102,328 KB
testcase_35 AC 157 ms
102,656 KB
testcase_36 AC 157 ms
102,236 KB
testcase_37 AC 157 ms
102,284 KB
testcase_38 AC 148 ms
102,264 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = lambda : sys.stdin.readline().rstrip()

sys.setrecursionlimit(2*10**5+10)
write = lambda x: sys.stdout.write(x+"\n")
debug = lambda x: sys.stderr.write(x+"\n")
writef = lambda x: print("{:.12f}".format(x))

def hurui(n):
    """線形篩
    pl: 素数のリスト
    mpf: iを割り切る最小の素因数
    """
    pl = []
    mpf = [None]*(n+1)
    for d in range(2,n+1):
        if mpf[d] is None:
            mpf[d] = d
            pl.append(d)
        for p in pl:
            if p*d>n or p>mpf[d]:
                break
            mpf[p*d] = p
    return pl, mpf
from collections import defaultdict
def factor(num):
    d = defaultdict(int)
    if num==1:
        d.update({1:1})
        return d
    while num>1:
        d[mpf[num]] += 1
        num //= mpf[num]
    return d
def fs(num):
    f = factor(num)
    ans = [1]
    for k,v in f.items():
        tmp = []
        for i in range(len(ans)):
            val = 1
            for _ in range(v):
                val *= k
                ans.append(ans[i]*val)
    return ans
def sub(x,p):
    res = 0
    while x%p==0:
        x //= p
        res += 1
    return res
inf = 10**15
from heapq import heappop as hpp, heappush as hp, heapify
def dijkstra(start,d):
    n = len(d)
    vals = [inf] * n
    h = [(0, start)] # (距離, ノード番号)
    vals[start] = 0
    # order = []
    while h:
        val, u = hpp(h)
        if val>vals[u]:
            continue
        # order.append(u)
        for v in range(u+1,n):
            dd = d[u][v]
#         for d,v in ns[u]:
            if vals[v]>val+dd:
                vals[v] = val+dd
                hp(h, (vals[v], v))
    return vals #, order

pl, mpf = hurui(10**6)
t = int(input())
ans = []
for i in range(t):
    x = int(input())
    for p in pl:
        if x%p!=0:
            res = p
            break
    else:
        assert 0
    m = 10
    count = [-1]*m # 回数がiの素数の最小
    for i in range(50)[::-1]:
        v = sub(x,pl[i])
        if v<m:
            count[v] = pl[i]
    inf = 10**9
    d = [[inf]*m for _ in range(m)]
    for i in range(m):
        v = count[i]
        if v==-1:
            continue
        cost = v
        for j in range(i+1,m):
            d[i][j] = cost
            cost *= v
    n = m
    for k in range(n):
        for i in range(k):
            for j in range(k+1,n):
                d[i][j] = min(d[i][j], d[i][k]*d[k][j])
    for i in range(m):
        if 2*i+1<m:
            res = min(res, d[i][2*i+1])
    ans.append(res*x)
write("\n".join(map(str, ans)))
0