結果

問題 No.2230 Good Omen of White Lotus
ユーザー wolgnikwolgnik
提出日時 2023-02-24 22:55:30
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 548 ms / 2,000 ms
コード長 1,785 bytes
コンパイル時間 1,449 ms
コンパイル使用メモリ 85,908 KB
実行使用メモリ 107,656 KB
最終ジャッジ日時 2023-10-11 06:45:36
合計ジャッジ時間 14,064 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 71 ms
70,560 KB
testcase_01 AC 73 ms
70,848 KB
testcase_02 AC 73 ms
70,832 KB
testcase_03 AC 73 ms
70,844 KB
testcase_04 AC 71 ms
70,796 KB
testcase_05 AC 72 ms
70,848 KB
testcase_06 AC 72 ms
70,800 KB
testcase_07 AC 72 ms
70,808 KB
testcase_08 AC 73 ms
70,848 KB
testcase_09 AC 73 ms
70,772 KB
testcase_10 AC 73 ms
70,896 KB
testcase_11 AC 73 ms
70,708 KB
testcase_12 AC 72 ms
70,712 KB
testcase_13 AC 84 ms
75,164 KB
testcase_14 AC 212 ms
78,412 KB
testcase_15 AC 90 ms
75,748 KB
testcase_16 AC 249 ms
79,480 KB
testcase_17 AC 252 ms
79,808 KB
testcase_18 AC 254 ms
79,528 KB
testcase_19 AC 255 ms
79,664 KB
testcase_20 AC 144 ms
77,860 KB
testcase_21 AC 114 ms
77,344 KB
testcase_22 AC 150 ms
78,324 KB
testcase_23 AC 143 ms
77,764 KB
testcase_24 AC 93 ms
82,988 KB
testcase_25 AC 93 ms
83,056 KB
testcase_26 AC 94 ms
83,412 KB
testcase_27 AC 264 ms
107,656 KB
testcase_28 AC 253 ms
104,204 KB
testcase_29 AC 304 ms
94,996 KB
testcase_30 AC 313 ms
95,164 KB
testcase_31 AC 335 ms
105,228 KB
testcase_32 AC 293 ms
105,196 KB
testcase_33 AC 548 ms
100,260 KB
testcase_34 AC 491 ms
99,368 KB
testcase_35 AC 507 ms
99,160 KB
testcase_36 AC 517 ms
99,536 KB
testcase_37 AC 494 ms
99,228 KB
testcase_38 AC 545 ms
100,180 KB
testcase_39 AC 544 ms
100,140 KB
testcase_40 AC 259 ms
93,648 KB
testcase_41 AC 266 ms
83,120 KB
testcase_42 AC 396 ms
95,512 KB
testcase_43 AC 161 ms
84,068 KB
testcase_44 AC 427 ms
95,048 KB
testcase_45 AC 266 ms
83,148 KB
testcase_46 AC 360 ms
92,452 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