結果

問題 No.2230 Good Omen of White Lotus
ユーザー wolgnikwolgnik
提出日時 2023-02-24 22:55:30
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 595 ms / 2,000 ms
コード長 1,785 bytes
コンパイル時間 175 ms
コンパイル使用メモリ 82,140 KB
実行使用メモリ 107,116 KB
最終ジャッジ日時 2024-09-13 05:58:11
合計ジャッジ時間 11,549 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
53,220 KB
testcase_01 AC 39 ms
54,424 KB
testcase_02 AC 38 ms
52,712 KB
testcase_03 AC 38 ms
52,296 KB
testcase_04 AC 39 ms
53,692 KB
testcase_05 AC 38 ms
53,256 KB
testcase_06 AC 38 ms
53,604 KB
testcase_07 AC 38 ms
53,268 KB
testcase_08 AC 39 ms
54,124 KB
testcase_09 AC 38 ms
53,656 KB
testcase_10 AC 38 ms
53,168 KB
testcase_11 AC 38 ms
53,584 KB
testcase_12 AC 40 ms
54,048 KB
testcase_13 AC 45 ms
60,300 KB
testcase_14 AC 167 ms
77,492 KB
testcase_15 AC 57 ms
68,228 KB
testcase_16 AC 216 ms
78,412 KB
testcase_17 AC 217 ms
78,824 KB
testcase_18 AC 214 ms
78,944 KB
testcase_19 AC 215 ms
78,384 KB
testcase_20 AC 111 ms
76,936 KB
testcase_21 AC 82 ms
76,144 KB
testcase_22 AC 114 ms
76,900 KB
testcase_23 AC 111 ms
76,724 KB
testcase_24 AC 62 ms
77,852 KB
testcase_25 AC 63 ms
78,100 KB
testcase_26 AC 65 ms
78,516 KB
testcase_27 AC 240 ms
107,116 KB
testcase_28 AC 225 ms
104,512 KB
testcase_29 AC 279 ms
94,884 KB
testcase_30 AC 289 ms
94,572 KB
testcase_31 AC 313 ms
104,772 KB
testcase_32 AC 270 ms
104,304 KB
testcase_33 AC 547 ms
98,512 KB
testcase_34 AC 499 ms
98,744 KB
testcase_35 AC 490 ms
98,872 KB
testcase_36 AC 508 ms
98,988 KB
testcase_37 AC 493 ms
98,496 KB
testcase_38 AC 521 ms
98,404 KB
testcase_39 AC 595 ms
98,748 KB
testcase_40 AC 239 ms
91,280 KB
testcase_41 AC 238 ms
81,976 KB
testcase_42 AC 372 ms
94,904 KB
testcase_43 AC 129 ms
83,292 KB
testcase_44 AC 413 ms
94,196 KB
testcase_45 AC 227 ms
81,840 KB
testcase_46 AC 325 ms
91,460 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.readline
H, W, N, P = map(int, input().split())
mod = 998244353
a = [[] for _ in range(H + 1)]

for _ in range(N):
  x, y = map(int, input().split())
  a[x].append(y)

class SegTree:
  def segfunc(self, x, y):
    return max(x, y)
  def __init__(self, n, ide_ele, init_val):
    #####単位元######
    self.ide_ele = ide_ele
    #num:n以上の最小の2のべき乗
    self.num = 2 ** (n - 1).bit_length()
    self.seg = [self.ide_ele] * 2 * self.num
    #set_val
    for i in range(n):
      self.seg[i + self.num - 1] = init_val[i]    
    #built
    for i in range(self.num - 2, -1, -1) :
      self.seg[i] = self.segfunc(self.seg[2 * i + 1], self.seg[2 * i + 2]) 
  def update(self, k, x):
    k += self.num - 1
    self.seg[k] = x
    while k + 1:
      k = (k - 1) // 2
      self.seg[k] = self.segfunc(self.seg[k * 2 + 1], self.seg[k * 2 + 2]) 
  def query(self, p, q):
    if q <= p:
      return self.ide_ele
    p += self.num - 1
    q += self.num - 2
    res = self.ide_ele
    while q - p > 1:
      if p & 1 == 0:
        res = self.segfunc(res, self.seg[p])
      if q & 1 == 1:
        res = self.segfunc(res, self.seg[q])
        q -= 1
      p = p // 2
      q = (q - 1) // 2
    if p == q:
      res = self.segfunc(res, self.seg[p])
    else:
      res = self.segfunc(self.segfunc(res, self.seg[p]), self.seg[q])
    return res

seg = SegTree(W + 1, 0, [0] * (W + 1))

for x in range(H, -1, -1):
  s = []
  a[x].sort(reverse = True)
  for y in a[x]:
    q = seg.query(y, W + 1)
    seg.update(y, q + 1)

  #print(s)
q = seg.query(0, W + 1)
res = pow(P - 2, q, mod) * pow(P, (mod - 2) * q, mod)
res %= mod
res *= pow(P - 1, (H + W - 3 - q), mod) * pow(P, (mod - 2) * (H + W - 3 - q), mod)
res %= mod
res = 1 - res
res %= mod
print(res)
0