結果

問題 No.14 最小公倍数ソート
コンテスト
ユーザー Navier_Boltzmann
提出日時 2023-07-02 17:53:41
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 262 ms / 5,000 ms
コード長 1,383 bytes
コンパイル時間 298 ms
コンパイル使用メモリ 82,008 KB
実行使用メモリ 88,028 KB
最終ジャッジ日時 2024-07-16 13:41:44
合計ジャッジ時間 5,187 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

from collections import *
from itertools import *
from functools import *
from heapq import *
import sys,math
input = sys.stdin.readline

N = int(input())
A = list(map(int,input().split()))
def md(n):
    lower_divisors , upper_divisors = [], []
    i = 1
    while i*i <= n:
        if n % i == 0:
            lower_divisors.append(i)
            if i != n // i:
                upper_divisors.append(n//i)
        i += 1
    return lower_divisors + upper_divisors[::-1]
    
X = defaultdict(list)
for i,a in enumerate(A):
    for p in md(a):
        X[p].append((a,i))

for k in X.keys():
    X[k].sort()
    X[k] = deque(X[k])

called = [False]*N
ans = [0]
called[0]=True
for _ in range(N-1):
    
    I = A[ans[-1]]
    
    idx = -1
    G = (1<<60)
    jval = (1<<60)
    
    for p in md(I):
        
        if len(X[p])==0:
            continue
        
        val,jdx = X[p].popleft()
        while called[jdx]:
            if len(X[p])==0:
                break
            val,jdx = X[p].popleft()
        if called[jdx]:
            continue
        X[p].appendleft((val,jdx))
        tmp = I*val//math.gcd(I,val)
        if tmp<G:
            G = tmp
            idx = jdx
            jval = val
        elif tmp==G:
            if val<jval:
                jval = val
                idx = jdx
    called[idx]=True

    ans.append(idx)
print(*[A[i] for i in ans])
0