結果
問題 | No.826 連絡網 |
ユーザー |
👑 ![]() |
提出日時 | 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))