結果

問題 No.2230 Good Omen of White Lotus
ユーザー 👑 KazunKazun
提出日時 2023-02-24 23:00:23
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 546 ms / 2,000 ms
コード長 4,945 bytes
コンパイル時間 348 ms
コンパイル使用メモリ 87,092 KB
実行使用メモリ 117,040 KB
最終ジャッジ日時 2023-10-11 06:47:03
合計ジャッジ時間 14,571 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 75 ms
71,360 KB
testcase_01 AC 74 ms
71,496 KB
testcase_02 AC 74 ms
71,400 KB
testcase_03 AC 76 ms
71,136 KB
testcase_04 AC 72 ms
71,408 KB
testcase_05 AC 73 ms
71,396 KB
testcase_06 AC 73 ms
71,344 KB
testcase_07 AC 72 ms
71,176 KB
testcase_08 AC 71 ms
71,344 KB
testcase_09 AC 73 ms
71,136 KB
testcase_10 AC 73 ms
71,276 KB
testcase_11 AC 74 ms
71,588 KB
testcase_12 AC 75 ms
71,332 KB
testcase_13 AC 74 ms
71,352 KB
testcase_14 AC 185 ms
79,240 KB
testcase_15 AC 87 ms
75,768 KB
testcase_16 AC 226 ms
79,780 KB
testcase_17 AC 227 ms
79,764 KB
testcase_18 AC 225 ms
79,788 KB
testcase_19 AC 229 ms
79,732 KB
testcase_20 AC 130 ms
78,112 KB
testcase_21 AC 107 ms
77,596 KB
testcase_22 AC 135 ms
78,200 KB
testcase_23 AC 130 ms
78,204 KB
testcase_24 AC 99 ms
95,044 KB
testcase_25 AC 98 ms
95,144 KB
testcase_26 AC 98 ms
95,160 KB
testcase_27 AC 284 ms
116,368 KB
testcase_28 AC 260 ms
116,644 KB
testcase_29 AC 296 ms
107,584 KB
testcase_30 AC 299 ms
107,740 KB
testcase_31 AC 325 ms
117,040 KB
testcase_32 AC 297 ms
117,040 KB
testcase_33 AC 524 ms
111,552 KB
testcase_34 AC 509 ms
111,432 KB
testcase_35 AC 544 ms
111,312 KB
testcase_36 AC 546 ms
111,216 KB
testcase_37 AC 529 ms
111,580 KB
testcase_38 AC 532 ms
111,372 KB
testcase_39 AC 532 ms
111,208 KB
testcase_40 AC 283 ms
103,912 KB
testcase_41 AC 224 ms
87,712 KB
testcase_42 AC 367 ms
94,560 KB
testcase_43 AC 148 ms
96,276 KB
testcase_44 AC 446 ms
101,048 KB
testcase_45 AC 251 ms
88,684 KB
testcase_46 AC 351 ms
103,780 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

class Segment_Tree():
    def __init__(self, L, calc, unit):
        """ calc を演算とするリスト L の Segment Tree を作成

        calc: 演算 (2変数関数, Monoid)
        unit: Monoid calc の単位元 (xe=ex=xを満たすe)
        """
        self.calc=calc
        self.unit=unit

        N=len(L); self.n=N
        d=max(1,(N-1).bit_length())
        k=1<<d

        self.data=data=[unit]*k+L+[unit]*(k-len(L))
        self.N=k
        self.depth=d

        for i in range(k-1,0,-1):
            data[i]=calc(data[i<<1], data[i<<1|1])

    def get(self, k):
        """ 第 k 要素を取得
        """
        assert 0<=k<self.N,"添字が範囲外"
        return self.data[k+self.N]

    def update(self, k, x):
        """第k要素をxに変え,更新を行う.

        k:数列の要素
        x:更新後の値
        """
        assert 0<=k<self.N,"添字が範囲外"
        m=k+self.N

        data=self.data; calc=self.calc
        data[m]=x

        while m>1:
            m>>=1
            data[m]=calc(data[m<<1], data[m<<1|1])

    def product(self, l, r, left_closed=True,right_closed=True):
        L=l+self.N+(not left_closed)
        R=r+self.N+(right_closed)

        vL=self.unit
        vR=self.unit

        data=self.data; calc=self.calc
        while L<R:
            if L&1:
                vL=calc(vL, data[L])
                L+=1

            if R&1:
                R-=1
                vR=calc(data[R], vR)

            L>>=1
            R>>=1

        return calc(vL,vR)

    def all_product(self):
        return self.data[1]

    def max_right(self, left, cond):
        """ 以下の2つをともに満たす x の1つを返す.\n
        (1) r=left or cond(data[left]*data[left+1]*...*data[r-1]): True
        (2) r=N or cond(data[left]*data[left+1]*...*data[r]): False
        ※ cond が単調減少の時, cond(data[left]*...*data[r-1]) を満たす最大の r となる.

        cond:関数(引数が同じならば結果も同じ)
        cond(unit): True
        0<=left<=N
        """

        assert 0<=left<=self.N,"添字が範囲外"
        assert cond(self.unit),"単位元が条件を満たさない."

        if left==self.N:
            return self.N

        left+=self.N
        sm=self.unit

        calc=self.calc; data=self.data
        first=True

        while first or (left & (-left))!=left:
            first=False
            while left%2==0:
                left>>=1
            if not cond(calc(sm, data[left])):
                while left<self.N:
                    left<<=1
                    if cond(calc(sm, data[left])):
                        sm=calc(sm, data[left])
                        left+=1
                return left-self.N
            sm=calc(sm, data[left])
            left+=1
        return self.N

    def min_left(self, right, cond):
        """ 以下の2つをともに満たす y の1つを返す.\n
        (1) l=right or cond(data[l]*data[l+1]*...*data[right-1]): True
        (2) l=0 or cond(data[l-1]*data[l]*...*data[right-1]): False
        ※ cond が単調増加の時, cond(data[l]*...*data[right-1]) を満たす最小の l となる.

        cond: 関数(引数が同じならば結果も同じ)
        cond(unit): True
        0<=right<=N
        """
        assert 0<=right<=self.N,"添字が範囲外"
        assert cond(self.unit),"単位元が条件を満たさない."

        if right==0:
            return 0

        right+=self.N
        sm=self.unit

        calc=self.calc; data=self.data
        first=1
        while first or (right & (-right))!=right:
            first=0
            right-=1
            while right>1 and right&1:
                right>>=1

            if not cond(calc(data[right], sm)):
                while right<self.N:
                    right=2*right+1
                    if cond(calc(data[right], sm)):
                        sm=calc(data[right], sm)
                        right-=1
                return right+1-self.N
            sm=calc(data[right], sm)
        return 0

    def __getitem__(self,k):
        return self.get(k)

    def __setitem__(self,k,x):
        return self.update(k,x)

    def __iter__(self):
        for i in range(self.n):
            yield self.get(i)

#==================================================
def solve():
    H,W,K,P=map(int,input().split())

    good=[[] for _ in range(H+1)]
    for k in range(K):
        x,y=map(int,input().split())
        good[x].append(y)

    S=Segment_Tree([0]*(W+1), max, 0)
    for x in range(H,-1,-1):
        good[x].sort(reverse=True)
        for y in good[x]:
            S[y]=S.product(y,W)+1
    M=S.all_product()

    Mod=998244353
    P_inv=pow(P,Mod-2,Mod)
    alpha=1-2*P_inv
    beta=1-P_inv
    return (1-pow(alpha,M,Mod)*pow(beta,H+W-3-M,Mod))%Mod

#==================================================
import sys
input=sys.stdin.readline
write=sys.stdout.write

print(solve())
0