結果

問題 No.1611 Minimum Multiple with Double Divisors
コンテスト
ユーザー persimmon-persimmon
提出日時 2022-04-01 14:21:30
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
TLE  
実行時間 -
コード長 1,336 bytes
コンパイル時間 268 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 31,232 KB
最終ジャッジ日時 2024-11-19 06:35:26
合計ジャッジ時間 55,668 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 18 TLE * 19
権限があれば一括ダウンロードができます

ソースコード

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