結果

問題 No.1059 素敵な集合
ユーザー maspymaspy
提出日時 2020-05-22 21:30:32
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 527 ms / 2,000 ms
コード長 1,260 bytes
コンパイル時間 103 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 22,400 KB
最終ジャッジ日時 2024-07-23 09:24:49
合計ジャッジ時間 4,733 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
10,880 KB
testcase_01 AC 527 ms
22,400 KB
testcase_02 AC 119 ms
15,488 KB
testcase_03 AC 29 ms
10,752 KB
testcase_04 AC 30 ms
10,752 KB
testcase_05 AC 30 ms
10,752 KB
testcase_06 AC 94 ms
13,952 KB
testcase_07 AC 95 ms
13,952 KB
testcase_08 AC 133 ms
15,104 KB
testcase_09 AC 50 ms
12,032 KB
testcase_10 AC 165 ms
15,488 KB
testcase_11 AC 103 ms
13,952 KB
testcase_12 AC 66 ms
12,928 KB
testcase_13 AC 185 ms
17,152 KB
testcase_14 AC 34 ms
11,136 KB
testcase_15 AC 328 ms
20,992 KB
testcase_16 AC 114 ms
14,592 KB
testcase_17 AC 108 ms
14,976 KB
testcase_18 AC 114 ms
22,400 KB
testcase_19 AC 521 ms
22,272 KB
testcase_20 AC 516 ms
22,400 KB
testcase_21 AC 327 ms
21,504 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