結果
問題 | No.1968 Distance |
ユーザー |
|
提出日時 | 2022-06-03 22:34:23 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 55 ms / 2,000 ms |
コード長 | 2,749 bytes |
コンパイル時間 | 215 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 62,720 KB |
最終ジャッジ日時 | 2024-09-21 03:01:09 |
合計ジャッジ時間 | 2,516 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 25 |
ソースコード
class UnionFindVerSize():def __init__(self, N):self._parent = [n for n in range(0, N)]self._size = [1] * Nself.group = Ndef find_root(self, x):if self._parent[x] == x: return xself._parent[x] = self.find_root(self._parent[x])stack = [x]while self._parent[stack[-1]]!=stack[-1]:stack.append(self._parent[stack[-1]])for v in stack:self._parent[v] = stack[-1]return self._parent[x]def unite(self, x, y):gx = self.find_root(x)gy = self.find_root(y)if gx == gy: returnself.group -= 1if self._size[gx] < self._size[gy]:self._parent[gx] = gyself._size[gy] += self._size[gx]return gyelse:self._parent[gy] = gxself._size[gx] += self._size[gy]return gxdef get_size(self, x):return self._size[self.find_root(x)]def is_same_group(self, x, y):return self.find_root(x) == self.find_root(y)import sys,random,bisectfrom collections import deque,defaultdictfrom heapq import heapify,heappop,heappushfrom itertools import permutationsfrom math import log,gcdinput = lambda :sys.stdin.readline().rstrip()mi = lambda :map(int,input().split())li = lambda :list(mi())def integral(l,r,f):#print(l,r,f)res = 0for d in range(len(f)):res += f[d] * (pow(r,d+1,mod)-pow(l,d+1,mod)) * pow(d+1,mod-2,mod)res %= modreturn resmod = 998244353N = int(input())res = 0for a in range(1,61):"""1/2^a <= x/(1-x) <= 1/2^(a-1)のとき(1-x)をa回する(完了したら1-x)→交互"""if a!=60:l = pow(pow(2,a,mod)+1,mod-2,mod)r = pow(pow(2,a-1,mod)+1,mod-2,mod)else:l = 0r = pow(pow(2,a-1,mod)+1,mod-2,mod)if 2**(a-1) > N:l = 0r = pow(pow(2,a-1,mod)+1,mod-2,mod)res += integral(l,r,[0,1])breaktmp_x = 1tmp_1_x = 2**(a-1)rest = 2**(a-1)-1p = 1while True:if p==1:if N <= rest+tmp_1_x:q = tmp_1_x * 2 % modiq = pow(q,mod-2,mod)res += integral(l,r,[iq,-iq])breakelse:rest += tmp_1_xtmp_1_x *= 2p = 1-pelse:if N <= rest+tmp_x:q = tmp_x * 2 % modiq = pow(q,mod-2,mod)res += integral(l,r,[0,iq])breakelse:rest += tmp_xtmp_x *= 2p = 1-pres = 2 * res % modprint(res)