結果

問題 No.826 連絡網
コンテスト
ユーザー c-yan
提出日時 2021-01-16 00:03:35
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 1,328 ms / 2,000 ms
コード長 1,176 bytes
コンパイル時間 560 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 114,420 KB
最終ジャッジ日時 2024-11-26 18:51:05
合計ジャッジ時間 13,019 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

from sys import setrecursionlimit, stdin


def find(parent, i):
    t = parent[i]
    if t < 0:
        return i
    t = find(parent, t)
    parent[i] = t
    return t


def unite(parent, i, j):
    i = find(parent, i)
    j = find(parent, j)
    if i == j:
        return
    parent[j] += parent[i]
    parent[i] = j


def make_prime_table(n):
    sieve = list(range(n + 1))
    sieve[0] = -1
    sieve[1] = -1
    for i in range(4, n + 1, 2):
        sieve[i] = 2
    for i in range(3, int(n ** 0.5) + 1, 2):
        if sieve[i] != i:
            continue
        for j in range(i * i, n + 1, i * 2):
            if sieve[j] == j:
                sieve[j] = i
    return sieve


def main():
    readline = stdin.readline
    setrecursionlimit(10 ** 6)

    N, P = map(int, readline().split())

    parent = [-1] * (N + 1)
    sieve = make_prime_table(N)
    for j in range(2 * 2, N + 1, 2):
        unite(parent, 2, j)
    for i in range(3, N + 1, 2):
        if sieve[i] != i:
            continue
        if 2 * i <= N:
            unite(parent, 2, i)
        for j in range(i * i, N + 1, i * 2):
            unite(parent, i, j)
    print(-parent[find(parent, P)])
main()
0