結果

問題 No.2756 GCD Teleporter
ユーザー titiatitia
提出日時 2024-05-14 04:24:24
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
実行時間 -
コード長 1,903 bytes
コンパイル時間 150 ms
コンパイル使用メモリ 12,928 KB
実行使用メモリ 50,228 KB
最終ジャッジ日時 2024-05-14 04:25:12
合計ジャッジ時間 42,610 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 449 ms
30,464 KB
testcase_01 AC 388 ms
30,464 KB
testcase_02 AC 420 ms
30,464 KB
testcase_03 AC 405 ms
30,464 KB
testcase_04 AC 482 ms
31,692 KB
testcase_05 AC 738 ms
36,164 KB
testcase_06 AC 1,335 ms
43,912 KB
testcase_07 AC 1,149 ms
38,376 KB
testcase_08 AC 1,117 ms
36,340 KB
testcase_09 AC 995 ms
37,784 KB
testcase_10 AC 1,226 ms
41,052 KB
testcase_11 AC 1,472 ms
47,856 KB
testcase_12 AC 584 ms
35,184 KB
testcase_13 AC 1,452 ms
49,384 KB
testcase_14 AC 1,321 ms
47,072 KB
testcase_15 AC 638 ms
34,996 KB
testcase_16 AC 910 ms
37,728 KB
testcase_17 AC 1,011 ms
38,484 KB
testcase_18 AC 595 ms
35,040 KB
testcase_19 AC 1,429 ms
50,228 KB
testcase_20 AC 1,278 ms
47,120 KB
testcase_21 AC 826 ms
37,208 KB
testcase_22 AC 1,389 ms
48,320 KB
testcase_23 AC 492 ms
32,076 KB
testcase_24 AC 1,100 ms
41,012 KB
testcase_25 AC 1,163 ms
42,976 KB
testcase_26 AC 1,262 ms
46,844 KB
testcase_27 AC 1,382 ms
48,420 KB
testcase_28 AC 921 ms
37,884 KB
testcase_29 AC 1,115 ms
41,012 KB
testcase_30 AC 1,353 ms
48,048 KB
testcase_31 AC 820 ms
37,092 KB
testcase_32 AC 1,090 ms
46,804 KB
testcase_33 AC 653 ms
34,056 KB
testcase_34 AC 528 ms
31,352 KB
testcase_35 AC 1,248 ms
43,908 KB
testcase_36 AC 1,173 ms
42,188 KB
testcase_37 AC 32 ms
11,008 KB
testcase_38 AC 1,932 ms
47,456 KB
testcase_39 TLE -
権限があれば一括ダウンロードができます

ソースコード

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