結果

問題 No.1611 Minimum Multiple with Double Divisors
ユーザー persimmon-persimmonpersimmon-persimmon
提出日時 2022-04-01 14:21:30
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
実行時間 -
コード長 1,336 bytes
コンパイル時間 118 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 25,600 KB
最終ジャッジ日時 2024-04-29 22:08:35
合計ジャッジ時間 7,992 ms
ジャッジサーバーID
(参考情報)
judge2 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 TLE -
testcase_01 -- -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #


# 指定した整数以下の素数を列挙。エラトステネスの篩
def eratosthenes(n):
    l0=list(range(2,n+1))
    l1=[1]*(n+1)
    l1[0]=0
    l1[1]=0
    for li in l0:
        if li>int(n**0.5):
            break
        if l1[li]==1:
            k=2
        while k*li<=n:
            l1[k*li]=0
            k+=1
    ret=[i for i,li in enumerate(l1) if li==1]
    return ret

# 素因数分解
def factorization(n,prime):
    arr = []
    temp = n
    for p in prime:
        if temp%p==0:
            cnt=0
            while temp%p==0:
                cnt+=1
                temp //= p
            arr.append([p, cnt])
        if temp==1:break
    if temp!=1:
        arr.append([temp, 1])
    if arr==[]:
        arr.append([n, 1])
    ps={p:q for p,q in arr}
    return ps
prime=eratosthenes(31)
prime.sort()
df={i:factorization(i,prime) for i in range(2,32)}
t=int(input())
xx=[int(input()) for _ in range(t)]
inf=float('inf')

for x in xx:
    ps=factorization(x,prime)
    for i in range(2,32):
        if i not in ps:ps[i]=0
    nx=1
    ans=inf
    for v in ps.values():nx*=v+1
    for i in range(2,32):
        for k,v in df[i].items():ps[k]+=v
        ny=1
        for v in ps.values():ny*=v+1
        if 2*nx==ny:
            ans=x*i
            break
        for k,v in df[i].items():ps[k]-=v
    print(ans)
0