結果

問題 No.826 連絡網
ユーザー brthyyjp
提出日時 2020-09-01 04:32:17
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 379 ms / 2,000 ms
コード長 1,044 bytes
コンパイル時間 433 ms
コンパイル使用メモリ 82,504 KB
実行使用メモリ 99,452 KB
最終ジャッジ日時 2024-11-17 08:37:03
合計ジャッジ時間 7,367 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #

n, p = map(int, input().split())

import math

N = 10**6+50
spf = [-1]*(N+1)
for i in range(N+1):
    spf[i] = i
for i in range(2, int(math.sqrt(N))+1):
    if spf[i] == i:
        for j in range(i*2, N+1, i):
            if spf[j] == j:
                spf[j] = i

def Find(x, par):
    if par[x] < 0:
        return x
    else:
        par[x] = Find(par[x], par)
        return par[x]

def Unite(x, y, par, rank):
    x = Find(x, par)
    y = Find(y, par)

    if x != y:
        if rank[x] < rank[y]:
            par[y] += par[x]
            par[x] = y
        else:
            par[x] += par[y]
            par[y] = x
            if rank[x] == rank[y]:
                rank[x] += 1

def Same(x, y, par):
    return Find(x, par) == Find(y, par)

def Size(x, par):
    return -par[Find(x, par)]

par = [-1]*(n+1)
rank = [0]*(n+1)
for i in range(1, n+1):
    j = i
    s = set()
    while j != 1:
        q = spf[j]
        if q not in s:
            Unite(i, q, par, rank)
            s.add(q)
        j //= q
#print(par)
print(Size(p, par))
0