結果

問題 No.826 連絡網
ユーザー stngstng
提出日時 2022-07-09 17:28:07
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,719 bytes
コンパイル時間 292 ms
コンパイル使用メモリ 82,364 KB
実行使用メモリ 117,248 KB
最終ジャッジ日時 2024-06-10 00:28:02
合計ジャッジ時間 3,695 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 35 ms
53,896 KB
testcase_01 AC 34 ms
54,048 KB
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import math

def sieve_of_eratosthenes(n):
    prime = [True for i in range(n+1)]
    prime[0] = False
    prime[1] = False

    sqrt_n = math.ceil(math.sqrt(n))
    for i in range(2, sqrt_n):
        if prime[i]:
            for j in range(2*i, n+1, i):
                prime[j] = False

    return prime

class UnionFind:
    def __init__(self, n):
        self.par = [i for i in range(n+1)]
        self.rank = [0] * (n+1)
        self.size = [1 for _ in range(n+1)]
    # 検索
    def find(self, x):
        if self.par[x] == x:
            return x
        else:
            self.par[x] = self.find(self.par[x])
            return self.par[x]
    # 併合
    def union(self, x, y):
        x = self.find(x)
        y = self.find(y)
        if self.rank[x] < self.rank[y]:
            self.par[x] = y
            self.size[y] += self.size[x]
        else:
            self.par[y] = x
            self.size[x] += self.size[y]
        if self.rank[x] == self.rank[y]:
            self.rank[x] += 1
    # 同じ集合に属するか判定
    def same_check(self, x, y):
        return self.find(x) == self.find(y)

n,p = map(int,input().split())
ki = UnionFind(n+1)
chk = [0]*(n+1)

prm = sieve_of_eratosthenes(n)
pm = []
for i in range(n+1):
    if prm[i] == True:
        pm.append(i)

tmp = 2
while tmp*2 <= n:
    #print(tmp*2,tmp,tmp)
    ki.union(2,tmp*2)
    tmp += 1

for i in range(1,len(pm)):
    now = pm[i]
    idx = i-1
    while now*pm[idx] <= n:
        #print(now,now,prm[idx])
        ki.union(now,now*pm[idx])
        idx += 1
        if idx >= len(pm):
            break

ans = 0
idx = ki.par[p]

for i in range(n):
    #print(ki.par[i+1])
    if ki.par[i+1] == idx:
        ans += 1

print(ans)
0