結果

問題 No.1059 素敵な集合
ユーザー maspymaspy
提出日時 2020-05-22 21:30:32
言語 Python3
(3.11.6 + numpy 1.26.0 + scipy 1.11.3)
結果
AC  
実行時間 455 ms / 2,000 ms
コード長 1,260 bytes
コンパイル時間 310 ms
コンパイル使用メモリ 11,032 KB
実行使用メモリ 19,768 KB
最終ジャッジ日時 2023-09-30 15:24:54
合計ジャッジ時間 4,246 ms
ジャッジサーバーID
(参考情報)
judge14 / judge15
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 16 ms
8,196 KB
testcase_01 AC 455 ms
19,592 KB
testcase_02 AC 96 ms
12,644 KB
testcase_03 AC 15 ms
8,268 KB
testcase_04 AC 15 ms
8,160 KB
testcase_05 AC 15 ms
8,336 KB
testcase_06 AC 75 ms
11,316 KB
testcase_07 AC 73 ms
11,168 KB
testcase_08 AC 101 ms
12,636 KB
testcase_09 AC 34 ms
9,596 KB
testcase_10 AC 140 ms
12,684 KB
testcase_11 AC 80 ms
11,368 KB
testcase_12 AC 49 ms
10,440 KB
testcase_13 AC 155 ms
14,520 KB
testcase_14 AC 20 ms
8,688 KB
testcase_15 AC 277 ms
18,228 KB
testcase_16 AC 89 ms
12,032 KB
testcase_17 AC 87 ms
12,152 KB
testcase_18 AC 86 ms
19,684 KB
testcase_19 AC 442 ms
19,768 KB
testcase_20 AC 426 ms
19,628 KB
testcase_21 AC 264 ms
18,784 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys

read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines

class UnionFind:
    def __init__(self, N):
        self.root = list(range(N))
        self.size = [1] * (N)
        self.n_comp = N

    def find_root(self, x):
        root = self.root
        while root[x] != x:
            root[x] = root[root[x]]
            x = root[x]
        return x

    def merge(self, x, y):
        x = self.find_root(x)
        y = self.find_root(y)
        if x == y:
            return False
        sx, sy = self.size[x], self.size[y]
        if sx < sy:
            self.root[x] = y
            self.size[y] += sx
        else:
            self.root[y] = x
            self.size[x] += sy
        self.n_comp -= 1
        return True

L, R = map(int, read().split())

U = R + 10

uf = UnionFind(R + 1)

is_prime = [0] * U
is_prime[2] = 1
for p in range(3, R + 1, 2):
    is_prime[p] = 1
for p in range(3, R + 1, 2):
    if p * p > R:
        break
    for i in range(p * p, R + 1, p + p):
        is_prime[i] = 0
primes = [p for p, x in enumerate(is_prime) if x]

for p in primes:
    for i in range(L, R + 1):
        if p * i > R:
            break
        uf.merge(i, p * i)

x = uf.n_comp - L - 1
print(x)
0