結果

問題 No.2230 Good Omen of White Lotus
ユーザー もろこし
提出日時 2025-10-03 15:20:23
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 2,263 bytes
コンパイル時間 522 ms
コンパイル使用メモリ 82,692 KB
実行使用メモリ 114,404 KB
最終ジャッジ日時 2025-10-03 15:21:00
合計ジャッジ時間 35,639 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 40 TLE * 4
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
sys.setrecursionlimit(10**7)
def ii(): return int(input())
def mi(d=0): return map(lambda x:int(x)-d,input().split())
INF = float("inf")
MOD = 998244353
def answer(s):
  print(s)
  exit()
def cross(xy):
  x,y = xy
  return (x+1,y),(x-1,y),(x,y+1),(x,y-1)
################################################
class SegTree:
  def __init__(self,n,inf,order,arr):
    self.inf = inf
    self.order = order
    self.bitlen = 2**n.bit_length()
    self.len = n
    self.n = self.bitlen*2-1
    arr += [inf]*(self.bitlen - n)
    self.arr = [None]*(self.bitlen-1) + arr

    def PreProcessing(index):
      if self.arr[index] == None:
        self.arr[index] = self.order(PreProcessing(2*index+1),PreProcessing(2*index+2))
      return self.arr[index]
    PreProcessing(0)
  
  #データ更新:index番目をxに変更
  def Renew(self,index,x):
    i = self.n - self.bitlen + index
    self.arr[i] = x
    while i != 0:
      i = (i-1)//2
      result = self.order(self.arr[i*2+1],self.arr[i*2+2])
      if result != self.arr[i]:
        self.arr[i] = result
  
  #index番目の値に+add
  def Add(self,index,add):
    i = self.n - self.bitlen + index
    self.arr[i] += add
    while i != 0:
      i = (i-1)//2
      self.arr[i] += add

  #[x,y)の範囲の結果を取得
  def GetRange(self,x,y,k=0,l=0,r=-1):
    if r == -1:
      r = self.bitlen
    if y<=l or r<=x:
      return self.inf
    elif x<=l and r<=y:
      return self.arr[k]
    else:
      return self.order(
        self.GetRange(x,y,2*k+1,l,(l+r)//2),
        self.GetRange(x,y,2*k+2,(l+r)//2,r)
        )
  
h,w,n,p = mi()
xy = [list(mi()) for _ in range(n)]
xy.sort()
# print(xy)
dp = SegTree(w,0,max,[0]*w)
for x,y in xy:
  # print(x,y)
  dp.Renew(y-1,dp.GetRange(0,y)+1)
  # print(dp.arr[-dp.bitlen:])
t = dp.GetRange(0,w)
# print(t)

move = h+w-3
p1 = p-1
a = 1
for b in range((move-t).bit_length()):
  if move-t>>b & 1:
    a = (a*p1) % MOD
  p1 = p1**2 % MOD

p2 = p-2
c = 1
for b in range(t.bit_length()):
  if t>>b & 1:
    c = (c*p2)%MOD
  p2 = p2**2 % MOD

print((1-((a*c)%MOD)*pow(p,-move,MOD))%MOD)

# q = ((p-1)**(move-t)%MOD * (p-2)**t%MOD)*pow(p,-move,MOD) %MOD
# print((1-q)%MOD)
0