結果

問題 No.2756 GCD Teleporter
ユーザー titiatitia
提出日時 2024-05-14 04:25:28
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 632 ms / 2,000 ms
コード長 1,903 bytes
コンパイル時間 405 ms
コンパイル使用メモリ 82,120 KB
実行使用メモリ 137,948 KB
最終ジャッジ日時 2024-12-20 09:50:23
合計ジャッジ時間 15,729 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 62 ms
68,384 KB
testcase_01 AC 63 ms
68,300 KB
testcase_02 AC 61 ms
67,916 KB
testcase_03 AC 63 ms
68,356 KB
testcase_04 AC 165 ms
82,804 KB
testcase_05 AC 261 ms
93,608 KB
testcase_06 AC 487 ms
118,156 KB
testcase_07 AC 424 ms
109,896 KB
testcase_08 AC 396 ms
107,908 KB
testcase_09 AC 351 ms
104,644 KB
testcase_10 AC 458 ms
113,748 KB
testcase_11 AC 550 ms
123,364 KB
testcase_12 AC 204 ms
86,420 KB
testcase_13 AC 540 ms
118,000 KB
testcase_14 AC 500 ms
118,100 KB
testcase_15 AC 214 ms
87,332 KB
testcase_16 AC 343 ms
102,108 KB
testcase_17 AC 387 ms
105,212 KB
testcase_18 AC 200 ms
86,592 KB
testcase_19 AC 561 ms
121,264 KB
testcase_20 AC 441 ms
110,820 KB
testcase_21 AC 268 ms
93,936 KB
testcase_22 AC 457 ms
111,712 KB
testcase_23 AC 154 ms
83,400 KB
testcase_24 AC 373 ms
102,996 KB
testcase_25 AC 391 ms
105,776 KB
testcase_26 AC 437 ms
111,224 KB
testcase_27 AC 450 ms
112,176 KB
testcase_28 AC 320 ms
100,080 KB
testcase_29 AC 418 ms
108,372 KB
testcase_30 AC 520 ms
118,836 KB
testcase_31 AC 290 ms
96,452 KB
testcase_32 AC 303 ms
109,368 KB
testcase_33 AC 176 ms
89,972 KB
testcase_34 AC 137 ms
83,380 KB
testcase_35 AC 360 ms
108,084 KB
testcase_36 AC 361 ms
107,752 KB
testcase_37 AC 40 ms
53,352 KB
testcase_38 AC 530 ms
113,496 KB
testcase_39 AC 632 ms
137,948 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.readline

N=int(input())
A=list(map(int,input().split()))

if set(A)=={1}:
    print(2*N)
    exit()

mod=998244353

# エラトステネスの篩を用いた素因数分解・約数列挙
MAX=5*10**5+10 # 使いたい最大値を指定

# Sieve[i]で、iの最も小さい約数を返す。
Sieve=[i for i in range(MAX)]

for i in range(2,MAX):
    if Sieve[i]!=i:
        continue
    
    for j in range(i,MAX,i):
        if Sieve[j]==j:
            Sieve[j]=i

# 素因数分解
def fact(x):
    D=dict()
    while x!=1:
        k=Sieve[x]
        if k in D:
            D[k]+=1
        else:
            D[k]=1
        x//=k
    return D

# 約数列挙
def faclist(x):
    LIST=[1]
    while x!=1:
        k=Sieve[x]
        count=0
        while x%k==0:
            count+=1
            x//=k

        LIST2=[]
        for l in LIST:
            for i in range(1,count+1):
                LIST2.append(l*k**i)
        LIST+=LIST2

    return LIST

LIST=[]
for a in A:
    LIST+=fact(a)

LIST=sorted(set(LIST))
D={LIST[i]:i for i in range(len(LIST))}

E=[[] for i in range(N+len(LIST))]

n=N+len(LIST)

# UnionFind

Group = [i for i in range(n)] # グループ分け
Nodes = [1]*(n) # 各グループのノードの数

def find(x):
    while Group[x] != x:
        x=Group[x]
    return x

def Union(x,y):
    if find(x) != find(y):
        if Nodes[find(x)] < Nodes[find(y)]:
            
            Nodes[find(y)] += Nodes[find(x)]
            Nodes[find(x)] = 0
            Group[find(x)] = find(y)
            
        else:
            Nodes[find(x)] += Nodes[find(y)]
            Nodes[find(y)] = 0
            Group[find(y)] = find(x)

for i in range(N):
    a=A[i]
    for u in fact(a):
        Union(i,D[u]+N)

#print([find(i) for i in range(n)])
k=LIST[0]

c=0
for i in range(n):
    if find(i)==i:
        c+=1

ANS=k*(c-1)
ANS2=2*c

print(min(ANS,ANS2))
0