結果
問題 | No.2833 Count Taiko Results |
ユーザー |
|
提出日時 | 2024-05-21 11:53:27 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 390 ms / 2,000 ms |
コード長 | 2,241 bytes |
コンパイル時間 | 187 ms |
コンパイル使用メモリ | 82,496 KB |
実行使用メモリ | 139,144 KB |
最終ジャッジ日時 | 2024-06-21 10:08:49 |
合計ジャッジ時間 | 10,225 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 58 |
ソースコード
class SegTree: def __init__(self, op, e, n, v=None): self._n = n self._op = op self._e = e self._log = (n - 1).bit_length() self._size = 1 << self._log self._d = [self._e] * (2 * self._size) if v is not None: for i in range(self._n): self._d[self._size + i] = v[i] for i in range(self._size - 1, 0, -1): self._update(i) def set(self, p, x): assert 0 <= p < self._n p += self._size self._d[p] = x for i in range(1, self._log + 1): self._update(p >> i) def get(self, p): assert 0 <= p < self._n return self._d[p + self._size] def prod(self, l, r): assert 0 <= l <= r <= self._n sml, smr = self._e, self._e l += self._size r += self._size while l < r: if l & 1: sml = self._op(sml, self._d[l]) l += 1 if r & 1: r -= 1 smr = self._op(self._d[r], smr) l >>= 1 r >>= 1 return self._op(sml, smr) def all_prod(self): return self._d[1] def _update(self, k): self._d[k] = self._op(self._d[2 * k], self._d[2 * k + 1]) #https://qiita.com/AkariLuminous/items/32cbf5bc3ffb2f84a898 N,K=list(map(int,input().split())) A=list(map(int,input().split())) B=list(map(int,input().split())) A.append(0) B.append(1) def MUL(a,b): return (a*b)%998244353 STST=SegTree(MUL,1,N,A) MOD=998244353 dp=[0]*(N+2)#Kコンボ未満 dp2=[0]*(N+2)#最大コンボがK rui=1 dp[0]=1 for i in range(1,N+2): if i>=K+1: rui=(rui-STST.prod(i-K-1,i-1)*dp[i-K-1])%MOD #K+1ノーツ目でミスするときはフルコンボ(0ノーツ目でミス)はつなげてはいけない。 dp[i]=(rui*B[i-1])%MOD #iノーツ目でミス rui=(rui*A[i-1])%MOD #iノーツ目をつないだ(コンボ数+1) rui=(rui+dp[i])%MOD #0コンボ rui=0 for i in range(N+2): if i>=K+2: rui=(rui-STST.prod(i-K-2,i-1)*dp2[i-K-2])%MOD dp2[i]=(rui*B[i-1])%MOD if i>=K+1: dp2[i]=(dp2[i]+B[i-1]*STST.prod(i-K-1,i-1)*dp[i-K-1])%MOD rui=(rui*A[i-1])%MOD rui=(rui+dp2[i])%MOD print(dp2[N+1])