結果

問題 No.2230 Good Omen of White Lotus
ユーザー flygonflygon
提出日時 2023-02-25 00:23:29
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 564 ms / 2,000 ms
コード長 2,312 bytes
コンパイル時間 162 ms
コンパイル使用メモリ 82,200 KB
実行使用メモリ 129,032 KB
最終ジャッジ日時 2024-09-13 06:23:06
合計ジャッジ時間 12,734 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
53,196 KB
testcase_01 AC 39 ms
53,220 KB
testcase_02 AC 40 ms
53,712 KB
testcase_03 AC 40 ms
53,360 KB
testcase_04 AC 39 ms
52,960 KB
testcase_05 AC 40 ms
52,888 KB
testcase_06 AC 39 ms
53,272 KB
testcase_07 AC 38 ms
53,120 KB
testcase_08 AC 39 ms
54,012 KB
testcase_09 AC 39 ms
53,244 KB
testcase_10 AC 39 ms
52,696 KB
testcase_11 AC 39 ms
52,908 KB
testcase_12 AC 40 ms
52,984 KB
testcase_13 AC 44 ms
54,064 KB
testcase_14 AC 234 ms
91,068 KB
testcase_15 AC 58 ms
65,964 KB
testcase_16 AC 316 ms
99,168 KB
testcase_17 AC 316 ms
99,152 KB
testcase_18 AC 310 ms
98,936 KB
testcase_19 AC 317 ms
98,872 KB
testcase_20 AC 114 ms
77,748 KB
testcase_21 AC 90 ms
76,800 KB
testcase_22 AC 118 ms
78,644 KB
testcase_23 AC 120 ms
79,228 KB
testcase_24 AC 61 ms
73,632 KB
testcase_25 AC 61 ms
73,464 KB
testcase_26 AC 62 ms
74,520 KB
testcase_27 AC 349 ms
129,032 KB
testcase_28 AC 362 ms
124,712 KB
testcase_29 AC 380 ms
123,812 KB
testcase_30 AC 385 ms
123,984 KB
testcase_31 AC 398 ms
124,920 KB
testcase_32 AC 400 ms
125,428 KB
testcase_33 AC 560 ms
119,828 KB
testcase_34 AC 543 ms
119,744 KB
testcase_35 AC 564 ms
119,892 KB
testcase_36 AC 553 ms
120,068 KB
testcase_37 AC 548 ms
119,468 KB
testcase_38 AC 543 ms
119,608 KB
testcase_39 AC 533 ms
119,352 KB
testcase_40 AC 253 ms
98,748 KB
testcase_41 AC 213 ms
88,212 KB
testcase_42 AC 370 ms
105,736 KB
testcase_43 AC 127 ms
84,920 KB
testcase_44 AC 453 ms
113,032 KB
testcase_45 AC 223 ms
88,356 KB
testcase_46 AC 345 ms
103,500 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def add(x,y):
  return x + y

def e(a):
  if a == min:
    return 10**18
  if a == max:
    return -10**18
  if a == add:
    return 0

class SegTree:
  def __init__(self, segf, init_val):
    n = len(init_val)
    self.segf = segf
    self.e = e(segf)
    self.seg_len = 1 << n.bit_length()
    self.seg = [self.e] * 2*self.seg_len
    
    for i in range(n):
      self.seg[i+ self.seg_len] = init_val[i]
    for i in range(self.seg_len)[::-1]:
      self.seg[i] = segf(self.seg[i<<1], self.seg[i<<1|1])

  def point_add(self,pos,x):
    pos += self.seg_len
    self.seg[pos] += x 
    while True: 
      pos >>= 1 
      if not pos: break 
      self.seg[pos] = self.segf(self.seg[pos<<1],  self.seg[pos<<1|1])

  def point_update(self, pos, x):
    pos += self.seg_len
    self.seg[pos] = x 
    while True: 
      pos >>= 1 
      if not pos: break 
      self.seg[pos] = self.segf(self.seg[pos<<1],  self.seg[pos<<1|1])

  def get_range(self,l,r):
    l += self.seg_len; r += self.seg_len
    res = self.e
    while l < r: 
      if l & 1: 
        res = self.segf(res, self.seg[l])
        l += 1 
      if r & 1: 
        r -= 1 
        res = self.segf(res, self.seg[r])
      l >>= 1; r >>= 1
    return res
  
  # ------ range_add & get_point -------
  def range_add(self, l,r,x):
    l += self.seg_len; r += self.seg_len
    self.seg[l] += x
    self.seg[r] += x
    while l < r: 
      if l & 1: 
        self.seg[l] = self.segf(x,self.seg[l])
        l += 1 
      if r & 1: 
        r -= 1 
        self.seg[r] = self.segf(x,self.seg[r])
      l >>= 1; r >>= 1

  def get_point(self,pos):
    pos += self.seg_len
    res = self.seg[pos]
    while True: 
      pos >>=1
      if not pos: break
      res = self.segf(res, self.seg[pos])
    return res
  
mod = 998244353
h,w,n,p = map(int,input().split())
xy = [list(map(int,input().split())) for i in range(n)]
xx = [[] for i in range(h+1)]
for x,y in xy:
  xx[x].append(y)

for i in range(1, h+1):
  xx[i].sort()

st = SegTree(max, [0]*(w+1))

for i in range(1, h+1):
  xxi = xx[i]
  for y in xxi:
    tmp = st.get_range(0, y+1)
    tmp += 1
    st.point_update(y, tmp)
  
cnt = st.get_range(0, w+1)

ex1 = (p-1)*pow(p, mod-2, mod)
ex2 = (p-2)*pow(p, mod-2, mod)

ans = 1 - pow(ex1, h+w-1-cnt-2, mod) * pow(ex2, cnt, mod)%mod
ans %= mod
print(ans)
0