結果
問題 | No.2872 Depth of the Parentheses |
ユーザー |
![]() |
提出日時 | 2024-09-07 15:16:03 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 163 ms / 2,000 ms |
コード長 | 23,151 bytes |
コンパイル時間 | 563 ms |
コンパイル使用メモリ | 82,420 KB |
実行使用メモリ | 208,164 KB |
最終ジャッジ日時 | 2024-09-07 15:16:22 |
合計ジャッジ時間 | 19,259 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 TLE * 5 |
ソースコード
def main():x,K=MI()global ansans=0prob1 = x*pow(100, -1, MOD)%MODprob0 = (100-x)*pow(100,-1,MOD)%MODrat=pow(prob1,K,MOD)*pow(prob0,K,MOD)%MODdef rec(pos, bit, now):global ansif pos==2*K:if now:returndep=0best=0for i in range(2*K):if bit>>i&1:dep+=1best=max(best,dep)else:dep-=1ans+=bestreturn# 1 ( を追加.rec(pos+1, bit+(1<<pos), now+1)# 0 ) を追加.if now>=1:rec(pos+1, bit, now-1)rec(0, 0, 0)print(ans*rat%MOD)# user config############DEBUG_MODE=1############# importimport sysimport itertoolsimport bisectimport mathfrom collections import *from functools import cachefrom heapq import *# aliasDD = defaultdictBSL = bisect.bisect_leftBSR = bisect.bisect_right# configinput = sys.stdin.readlinesys.setrecursionlimit(10**7)# inputdef II(): return int(input())def IS(): return input()[:-1]def MI(): return map(int,input().split())def LM(): return list(MI())def LL(n): return [LM() for _ in range(n)]def MI_1(): return map(lambda x:int(x)-1,input().split())def LM_1(): return list(MI_1())def LL_1(n): return [LM_1() for _ in range(n)]def ALPHABET_TO_NUM(string, upper=False): return list(map(lambda elm:ord(elm)-ord("A") if upper else ord(elm)-ord("a"), string))# functionsdef DB(*args,**kwargs):global DEBUG_MODEif not DEBUG_MODE:returnif args:print(*args)returnfor name, value in kwargs.items():print(f"{name} : {value}")def bit_count(num):length = num.bit_length()res = 0for i in range(length):if num >> i & 1:res += 1return resdef popcount64(n):# 63桁まで高速に動く.64桁まで正常に動く.c=(n&0x5555555555555555)+((n>>1)&0x5555555555555555)c=(c&0x3333333333333333)+((c>>2)&0x3333333333333333)c=(c&0x0f0f0f0f0f0f0f0f)+((c>>4)&0x0f0f0f0f0f0f0f0f)c=(c&0x00ff00ff00ff00ff)+((c>>8)&0x00ff00ff00ff00ff)c=(c&0x0000ffff0000ffff)+((c>>16)&0x0000ffff0000ffff)c=(c&0x00000000ffffffff)+((c>>32)&0x00000000ffffffff)return cdef argmax(*args):if len(args) == 1 and hasattr(args[0], '__iter__'):lst = args[0]else:lst = argsreturn lst.index(max(lst))def argmin(*args):if len(args) == 1 and hasattr(args[0], '__iter__'):lst = args[0]else:lst = argsreturn lst.index(min(lst))def prefix_op(lst, op=lambda x,y:x+y, e=0):N = len(lst)res = [e]*(N+1)for i in range(N):res[i+1] = op(res[i], lst[i])return resdef suffix_op(lst, op=lambda x,y:x+y, e=0):N = len(lst)res = [e]*(N+1)for i in range(N):res[N-1-i] = op(lst[N-1-i], res[N-i])return resdef sigma_LinearFunc(coeff1, coeff0, left, right, MOD=None):"""coeff1*x + coeff0の x = [left, right] の和を求める.MODで計算したい場合、区間の引数をMOD取った上で代入しても良い.そのとき、left > right となってもよい."""if MOD:# MODが素数でない場合にも対応するように、和公式を適応後に剰余を計算.return ((coeff0%MOD*((right-left+1)%MOD)%MOD) + (coeff1%MOD*((left+right)*(right-left+1)//2%MOD)%MOD))%MODreturn coeff0*(right-left+1) + coeff1*(left+right)*(right-left+1)//2def find_divisors(n):divs_small, divs_big = [], []i = 1while i * i <= n:if n % i == 0:divs_small.append(i)divs_big.append(n // i)i += 1if divs_big[-1] == divs_small[-1]:divs_big.pop()for e in reversed(divs_big):divs_small.append(e)return divs_smalldef prime_factorization(n):n_intact = nret = []i = 2while i * i <= n_intact:if n % i == 0:cnt = 0while n % i == 0:n //= icnt += 1ret.append((i,cnt))i += 1if n != 1: ret.append((n,1))return ret""" 矩形の二次元配列を扱う諸関数 """def copy_table(table):H,W = len(table), len(table[0])res = []for i in range(H):res.append([])for j in range(W):res[-1].append(table[i][j])return resdef sum_table(table, MOD=None):H,W = len(table), len(table[0])res = 0for i in range(H):for j in range(W):res += table[i][j]if MOD:res %= MODreturn resdef transpose_table(table):H,W = len(table), len(table[0])res = [[None]*H for _ in range(W)]for i in range(H):for j in range(W):res[j][i] = table[i][j]return resdef convert_table_to_bit(table, letter1="#", rev=False):H,W = len(table), len(table[0])res = []for h in range(H):rowBit = 0for w in range(W):if rev:if table[h][w] == letter1:rowBit += 2**welse:if table[h][W-w-1] == letter1:rowBit += 2**wres.append(rowBit)return resdef rotate_table_cc(S): return list(zip(*S))[::-1]def rotate_table_c(S): return list(zip(*S[::-1]))def mul_matrix(A, B, mod):N = len(A)K = len(A[0])M = len(B[0])res = [[0 for _ in range(M)] for _ in range(N)]for i in range(N) :for j in range(K) :for k in range(M) :res[i][k] += A[i][j] * B[j][k]res[i][k] %= modreturn resdef pow_matrix(mat, exp, mod):N = len(mat)res = [[1 if i == j else 0 for i in range(N)] for j in range(N)]while exp > 0 :if exp%2 == 1 :res = mul_matrix(res, mat, mod)mat = mul_matrix(mat, mat, mod)exp //= 2return resdef compress(lst):D = {e:i for i,e in enumerate(sorted(set(lst)))}return [D[e] for e in lst]def highDimCompress(lst):#(x,y)の配列や,(x,y,z)の配列が与えられたとき,軸ごとに座圧する.return list(zip(*list(map(compress,list(zip(*lst))))))#classes# https://github.com/tatyam-prime/SortedSet/blob/main/SortedSet.pyfrom bisect import bisect_left, bisect_rightfrom typing import Generic, Iterable, Iterator, List, Tuple, TypeVar, OptionalT = TypeVar('T')class SortedSet(Generic[T]):BUCKET_RATIO = 16SPLIT_RATIO = 24def __init__(self, a: Iterable[T] = []) -> None:"Make a new SortedSet from iterable. / O(N) if sorted and unique / O(N log N)"a = list(a)n = self.size = len(a)if any(a[i] > a[i + 1] for i in range(n - 1)):a.sort()if any(a[i] >= a[i + 1] for i in range(n - 1)):a, b = [], afor x in b:if not a or a[-1] != x:a.append(x)bucket_size = int(math.ceil(math.sqrt(n / self.BUCKET_RATIO)))self.a = [a[n * i // bucket_size : n * (i + 1) // bucket_size] for i in range(bucket_size)]def __iter__(self) -> Iterator[T]:for i in self.a:for j in i: yield jdef __reversed__(self) -> Iterator[T]:for i in reversed(self.a):for j in reversed(i): yield jdef __eq__(self, other) -> bool:return list(self) == list(other)def __len__(self) -> int:return self.sizedef __repr__(self) -> str:return "SortedSet" + str(self.a)def __str__(self) -> str:s = str(list(self))return "{" + s[1 : len(s) - 1] + "}"def _position(self, x: T) -> Tuple[List[T], int, int]:"return the bucket, index of the bucket and position in which x should be. self must not be empty."for i, a in enumerate(self.a):if x <= a[-1]: breakreturn (a, i, bisect_left(a, x))def __contains__(self, x: T) -> bool:if self.size == 0: return Falsea, _, i = self._position(x)return i != len(a) and a[i] == xdef add(self, x: T) -> bool:"Add an element and return True if added. / O(√N)"if self.size == 0:self.a = [[x]]self.size = 1return Truea, b, i = self._position(x)if i != len(a) and a[i] == x: return Falsea.insert(i, x)self.size += 1if len(a) > len(self.a) * self.SPLIT_RATIO:mid = len(a) >> 1self.a[b:b+1] = [a[:mid], a[mid:]]return Truedef _pop(self, a: List[T], b: int, i: int) -> T:ans = a.pop(i)self.size -= 1if not a: del self.a[b]return ansdef discard(self, x: T) -> bool:"Remove an element and return True if removed. / O(√N)"if self.size == 0: return Falsea, b, i = self._position(x)if i == len(a) or a[i] != x: return Falseself._pop(a, b, i)return Truedef lt(self, x: T) -> Optional[T]:"Find the largest element < x, or None if it doesn't exist."for a in reversed(self.a):if a[0] < x:return a[bisect_left(a, x) - 1]def le(self, x: T) -> Optional[T]:"Find the largest element <= x, or None if it doesn't exist."for a in reversed(self.a):if a[0] <= x:return a[bisect_right(a, x) - 1]def gt(self, x: T) -> Optional[T]:"Find the smallest element > x, or None if it doesn't exist."for a in self.a:if a[-1] > x:return a[bisect_right(a, x)]def ge(self, x: T) -> Optional[T]:"Find the smallest element >= x, or None if it doesn't exist."for a in self.a:if a[-1] >= x:return a[bisect_left(a, x)]def __getitem__(self, i: int) -> T:"Return the i-th element."if i < 0:for a in reversed(self.a):i += len(a)if i >= 0: return a[i]else:for a in self.a:if i < len(a): return a[i]i -= len(a)raise IndexErrordef pop(self, i: int = -1) -> T:"Pop and return the i-th element."if i < 0:for b, a in enumerate(reversed(self.a)):i += len(a)if i >= 0: return self._pop(a, ~b, i)else:for b, a in enumerate(self.a):if i < len(a): return self._pop(a, b, i)i -= len(a)raise IndexErrordef index(self, x: T) -> int:"Count the number of elements < x."ans = 0for a in self.a:if a[-1] >= x:return ans + bisect_left(a, x)ans += len(a)return ansdef index_right(self, x: T) -> int:"Count the number of elements <= x."ans = 0for a in self.a:if a[-1] > x:return ans + bisect_right(a, x)ans += len(a)return ans"""(num, cnt)を要素としたSSを管理してMultiset化したいとき用."""def exist(self, x):ret = self.gt((x,0))if ret is None:return Falseelif ret[0] == x:return Trueelse:return Falsedef increment(self, x):if not self.exist(x):self.add((x,1))else:num, cnt = self.gt((x,0))self.discard((x,cnt))self.add((x,cnt+1))def decrement(self, x):if not self.exist(x):returnnum, cnt = self.gt((x,0))if cnt == 1:self.discard((x,cnt))else:self.discard((x,cnt))self.add((x,cnt-1))def multi_add(self, x, y):if not self.exist(x):self.add((x,y))else:num, cnt = self.gt((x,0))self.discard((x,cnt))self.add((x,cnt+y))def multi_sub(self, x, y):if not self.exist(x):returnnum, cnt = self.gt((x,0))if cnt <= y:self.discard((x,cnt))else:self.discard((x,cnt))self.add((x,cnt-y))# https://github.com/tatyam-prime/SortedSet/blob/main/SortedMultiset.pyT = TypeVar('T')class SortedMultiset(Generic[T]):BUCKET_RATIO = 16SPLIT_RATIO = 24def __init__(self, a: Iterable[T] = []) -> None:"Make a new SortedMultiset from iterable. / O(N) if sorted / O(N log N)"a = list(a)n = self.size = len(a)if any(a[i] > a[i + 1] for i in range(n - 1)):a.sort()num_bucket = int(math.ceil(math.sqrt(n / self.BUCKET_RATIO)))self.a = [a[n * i // num_bucket : n * (i + 1) // num_bucket] for i in range(num_bucket)]def __iter__(self) -> Iterator[T]:for i in self.a:for j in i: yield jdef __reversed__(self) -> Iterator[T]:for i in reversed(self.a):for j in reversed(i): yield jdef __eq__(self, other) -> bool:return list(self) == list(other)def __len__(self) -> int:return self.sizedef __repr__(self) -> str:return "SortedMultiset" + str(self.a)def __str__(self) -> str:s = str(list(self))return "{" + s[1 : len(s) - 1] + "}"def _position(self, x: T) -> Tuple[List[T], int, int]:"return the bucket, index of the bucket and position in which x should be. self must not be empty."for i, a in enumerate(self.a):if x <= a[-1]: breakreturn (a, i, bisect_left(a, x))def __contains__(self, x: T) -> bool:if self.size == 0: return Falsea, _, i = self._position(x)return i != len(a) and a[i] == xdef count(self, x: T) -> int:"Count the number of x."return self.index_right(x) - self.index(x)def add(self, x: T) -> None:"Add an element. / O(√N)"if self.size == 0:self.a = [[x]]self.size = 1returna, b, i = self._position(x)a.insert(i, x)self.size += 1if len(a) > len(self.a) * self.SPLIT_RATIO:mid = len(a) >> 1self.a[b:b+1] = [a[:mid], a[mid:]]def _pop(self, a: List[T], b: int, i: int) -> T:ans = a.pop(i)self.size -= 1if not a: del self.a[b]return ansdef discard(self, x: T) -> bool:"Remove an element and return True if removed. / O(√N)"if self.size == 0: return Falsea, b, i = self._position(x)if i == len(a) or a[i] != x: return Falseself._pop(a, b, i)return Truedef lt(self, x: T) -> Optional[T]:"Find the largest element < x, or None if it doesn't exist."for a in reversed(self.a):if a[0] < x:return a[bisect_left(a, x) - 1]def le(self, x: T) -> Optional[T]:"Find the largest element <= x, or None if it doesn't exist."for a in reversed(self.a):if a[0] <= x:return a[bisect_right(a, x) - 1]def gt(self, x: T) -> Optional[T]:"Find the smallest element > x, or None if it doesn't exist."for a in self.a:if a[-1] > x:return a[bisect_right(a, x)]def ge(self, x: T) -> Optional[T]:"Find the smallest element >= x, or None if it doesn't exist."for a in self.a:if a[-1] >= x:return a[bisect_left(a, x)]def __getitem__(self, i: int) -> T:"Return the i-th element."if i < 0:for a in reversed(self.a):i += len(a)if i >= 0: return a[i]else:for a in self.a:if i < len(a): return a[i]i -= len(a)raise IndexErrordef pop(self, i: int = -1) -> T:"Pop and return the i-th element."if i < 0:for b, a in enumerate(reversed(self.a)):i += len(a)if i >= 0: return self._pop(a, ~b, i)else:for b, a in enumerate(self.a):if i < len(a): return self._pop(a, b, i)i -= len(a)raise IndexErrordef index(self, x: T) -> int:"Count the number of elements < x."ans = 0for a in self.a:if a[-1] >= x:return ans + bisect_left(a, x)ans += len(a)return ansdef index_right(self, x: T) -> int:"Count the number of elements <= x."ans = 0for a in self.a:if a[-1] > x:return ans + bisect_right(a, x)ans += len(a)return ansclass Comb:def __init__(self,table_len,mod):"""modが素数の場合しか使用できない.table_len に指定した数まで法modでのコンビネーションの計算が可能になる."""self.fac = [1,1]self.inv = [1,1]self.finv = [1,1]self.mod = modfor i in range(2,table_len+1):self.fac.append(self.fac[i-1]*i%mod)self.inv.append(-self.inv[mod%i]*(mod//i)%mod)self.finv.append(self.finv[i-1]*self.inv[i]%mod)def comb(self,a,b):return self.fac[a]*self.finv[b]*self.finv[a-b]%self.modclass GridBFS:def __init__(self, table):#二次元配列や文字列の配列を受け取る.self.table = tableself.H = len(table)self.W = len(table[0])self.wall = "#"def find(self, c):#table から引数の文字を探しインデックスを返す. 無い時、None.for h in range(self.H):for w in range(self.W):if self.table[h][w] == c:return (h,w)return Nonedef set_wall_string(self, string):#壁として扱う文字を "#!^" の様に文字列リテラルで格納. 初期値は、"#"self.wall = stringdef island(self, transition = [[-1,0],[0,1],[1,0],[0,-1]]):H, W = self.H, self.Wself.island_id = [[-1]*W for _ in range(H)]self.island_size = [[-1]*W for _ in range(W)]crr_id = 0id2size = dict()for sh in range(H):for sw in range(W):if self.table[sh][sw] in self.wall:continueif self.island_id[sh][sw] != -1:continuedeq = deque()deq.append((sh,sw))crr_size = 1self.island_id[sh][sw] = crr_idwhile deq:h,w = deq.popleft()for dh, dw in transition:nh, nw = h+dh, w+dwif (not 0<=nh<H) or (not 0<=nw<W):continueif self.table[nh][nw] in self.wall:continueif self.island_id[nh][nw] == -1:self.island_id[nh][nw] = crr_iddeq.append((nh, nw))crr_size += 1id2size[crr_id] = crr_sizecrr_id += 1for h in range(H):for w in range(W):if self.table[h][w] in self.wall:continueself.island_size[h][w] = id2size[self.island_id[h][w]]return self.island_id, self.island_sizedef distance(self, start, goal=None, transition = [[-1,0],[0,1],[1,0],[0,-1]]):#goal指定したら、goalまでの最短距離を、指定しなければdist配列を返す. 到達不可能は -1.#二次元配列上での遷移方法を transition で指定できる. 初期値は上下左右.H, W, tab, wall = self.H, self.W, self.table, self.wallINF = 1<<60deq = deque()deq.append(start)dist = [[INF]*W for _ in range(H)]dist[start[0]][start[1]] = 0if start == goal:return 0while deq:h,w = deq.popleft()for dh, dw in transition:nh = h+dhnw = w+dw# gridの範囲外.if (not 0<=nh<H) or (not 0<=nw<W):continue# wallに設定されている文字なら.if tab[nh][nw] in wall:continuenew_dist = dist[h][w] + 1#goalが引数で与えられていてgoalに達したら終了.if goal and (nh,nw)==goal:return new_distif dist[nh][nw] > new_dist:dist[nh][nw] = new_distdeq.append((nh,nw))# goalが設定されていていまだreturnされていないなら,# goalに到達できなかったということ.if goal:return -1return distclass nth_root:def __init__(self):self.ngs = [-1, -1, 4294967296, 2642246, 65536, 7132, 1626, 566, 256, 139, 85, 57, 41, 31, 24, 20, 16, 14, 12, 11, 10, 9, 8, 7, 7, 6, 6, 6, 5, 5, 5, 5, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3]""" pre-calculated by this code """# self.ngs = [-1,-1]# for n in range(2,64):# ok,ng=2**33,0# while abs(ok-ng)>1:# mid=ok+ng>>1# if mid**n>=2**64:# ok=mid# else:# ng=mid# self.ngs.append(ok)""" xが64bitに収まる時高速. そうでないときも動作するが速度の保証無し. """def calc(self, x, n, is_x_within_64bit=True):if x<=1 or n==1: return xif is_x_within_64bit:if n>=64: return 1ng = self.ngs[n]else:ng = xok = 0while abs(ok-ng)>1:mid = ok+ng>>1if mid**n <= x:ok = midelse:ng = midreturn ok# well-used constglobal DIRECTION_4, DIRECTION_8, DIRECTION_DIAGONAL, DIRECTION_URDL_TABLE, DIRECTION_URDL_COORD_PLANE, MOD, INF, LOWER_ALP, UPPER_ALP, ALL_ALP# clockwise from top.DIRECTION_4 = [[-1,0],[0,1],[1,0],[0,-1]]DIRECTION_8 = [[-1,0],[-1,1],[0,1],[1,1],[1,0],[1,-1],[0,-1],[-1,-1]]DIRECTION_DIAGONAL = [[-1,1],[1,1],[1,-1],[-1,-1]]DIRECTION_URDL_TABLE = {'U':(-1,0), 'R':(0,1), 'D':(1,0), 'L':(0,-1)}DIRECTION_URDL_COORD_PLANE = {'U':(0,1), 'R':(1,0), 'D':(0,-1), 'L':(-1,0)}MOD = 998244353INF = 1<<60LOWER_ALP = "abcdefghijklmnopqrstuvwxyz"UPPER_ALP = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"ALL_ALP = "abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZ"main()