結果

問題 No.2756 GCD Teleporter
ユーザー titiatitia
提出日時 2024-05-14 04:25:28
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 628 ms / 2,000 ms
コード長 1,903 bytes
コンパイル時間 170 ms
コンパイル使用メモリ 82,556 KB
実行使用メモリ 137,896 KB
最終ジャッジ日時 2024-05-14 04:25:47
合計ジャッジ時間 14,634 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 64 ms
67,484 KB
testcase_01 AC 63 ms
67,944 KB
testcase_02 AC 61 ms
67,596 KB
testcase_03 AC 61 ms
67,708 KB
testcase_04 AC 142 ms
83,068 KB
testcase_05 AC 295 ms
93,744 KB
testcase_06 AC 471 ms
118,520 KB
testcase_07 AC 393 ms
110,228 KB
testcase_08 AC 395 ms
108,048 KB
testcase_09 AC 346 ms
104,636 KB
testcase_10 AC 436 ms
114,032 KB
testcase_11 AC 554 ms
123,652 KB
testcase_12 AC 190 ms
86,748 KB
testcase_13 AC 488 ms
118,272 KB
testcase_14 AC 472 ms
118,400 KB
testcase_15 AC 207 ms
87,480 KB
testcase_16 AC 354 ms
102,324 KB
testcase_17 AC 369 ms
105,432 KB
testcase_18 AC 194 ms
86,672 KB
testcase_19 AC 518 ms
121,960 KB
testcase_20 AC 413 ms
111,168 KB
testcase_21 AC 276 ms
94,076 KB
testcase_22 AC 433 ms
112,184 KB
testcase_23 AC 147 ms
83,604 KB
testcase_24 AC 339 ms
103,068 KB
testcase_25 AC 381 ms
105,832 KB
testcase_26 AC 409 ms
111,512 KB
testcase_27 AC 457 ms
112,308 KB
testcase_28 AC 307 ms
100,456 KB
testcase_29 AC 386 ms
108,656 KB
testcase_30 AC 483 ms
118,772 KB
testcase_31 AC 279 ms
96,772 KB
testcase_32 AC 306 ms
109,292 KB
testcase_33 AC 165 ms
90,276 KB
testcase_34 AC 134 ms
83,456 KB
testcase_35 AC 332 ms
108,288 KB
testcase_36 AC 335 ms
107,984 KB
testcase_37 AC 37 ms
52,476 KB
testcase_38 AC 490 ms
113,996 KB
testcase_39 AC 628 ms
137,896 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