結果
| 問題 |
No.826 連絡網
|
| コンテスト | |
| ユーザー |
👑 Kazun
|
| 提出日時 | 2021-02-10 21:16:06 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 502 ms / 2,000 ms |
| コード長 | 3,090 bytes |
| コンパイル時間 | 148 ms |
| コンパイル使用メモリ | 82,308 KB |
| 実行使用メモリ | 103,936 KB |
| 最終ジャッジ日時 | 2024-07-08 08:30:58 |
| 合計ジャッジ時間 | 6,907 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 30 |
ソースコード
class Union_Find():
def __init__(self,N):
"""0,1,...,n-1を要素として初期化する.
N:要素数
"""
self.n=N
self.parents=[-1]*N
self.rank=[0]*N
def find(self, x):
"""要素xの属している族を調べる.
x:要素
"""
V=[]
while self.parents[x]>=0:
V.append(x)
x=self.parents[x]
for v in V:
self.parents[v]=x
return x
def union(self, x, y):
"""要素x,yを同一視する.
x,y:要素
"""
x=self.find(x)
y=self.find(y)
if x==y:
return
if self.rank[x]<self.rank[y]:
x,y=y,x
self.parents[x]+=self.parents[y]
self.parents[y]=x
if self.rank[x]==self.rank[y]:
self.rank[x]+=1
def size(self, x):
"""要素xの属している要素の数.
x:要素
"""
return -self.parents[self.find(x)]
def same(self, x, y):
"""要素x,yは同一視されているか?
x,y:要素
"""
return self.find(x) == self.find(y)
def members(self, x):
"""要素xが属している族の要素.
※族の要素の個数が欲しいときはsizeを使うこと!!
x:要素
"""
root = self.find(x)
return [i for i in range(self.n) if self.find(i) == root]
def roots(self):
"""族の名前のリスト
"""
return [i for i, x in enumerate(self.parents) if x < 0]
def group_count(self):
"""族の個数
"""
return len(self.roots())
def all_group_members(self):
"""全ての族の出力
"""
X={r:[] for r in self.roots()}
for k in range(self.n):
X[self.find(k)].append(k)
return X
def refresh(self):
for i in range(self.n):
_=self.find(i)
def __str__(self):
return '\n'.join('{}: {}'.format(r, self.members(r)) for r in self.roots())
def __repr__(self):
return self.__str__()
def Sieve_of_Eratosthenes(N,mode=False):
"""Nまでのエラトステネスの篩を実行
N:自然数
mode:False->素数のリスト,True->素数かどうかのリスト
(False->[2,3,5,...],True->[False,False,True,True,False,True,...])
"""
if N==0:
return [None]
T=[True]*(N+1)
T[0]=None
T[1]=False
x=4
while x<=N:
T[x]=False
x+=2
x=9
while x<=N:
T[x]=False
x+=6
a=5
Flag=1
while a*a<=N:
if T[a]:
b=a*a
c=2*a
while b<=N:
T[b]=False
b+=c
a+=2 if Flag else 4
Flag^=1
if mode:
return T
else:
return [k for k in range(N+1) if T[k]]
#================================================
N,P=map(int,input().split())
U=Union_Find(N+1)
A=Sieve_of_Eratosthenes(N)
for p in A:
q=2*p
while q<=N:
U.union(p,q)
q+=p
print(U.size(P))
Kazun