結果

問題 No.2230 Good Omen of White Lotus
ユーザー KazunKazun
提出日時 2023-02-24 23:00:23
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 494 ms / 2,000 ms
コード長 4,945 bytes
コンパイル時間 185 ms
コンパイル使用メモリ 82,108 KB
実行使用メモリ 115,808 KB
最終ジャッジ日時 2024-09-13 05:59:28
合計ジャッジ時間 10,789 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
53,968 KB
testcase_01 AC 43 ms
53,108 KB
testcase_02 AC 39 ms
53,200 KB
testcase_03 AC 38 ms
53,148 KB
testcase_04 AC 38 ms
52,828 KB
testcase_05 AC 38 ms
53,368 KB
testcase_06 AC 40 ms
54,160 KB
testcase_07 AC 38 ms
53,972 KB
testcase_08 AC 39 ms
54,136 KB
testcase_09 AC 38 ms
53,488 KB
testcase_10 AC 39 ms
54,000 KB
testcase_11 AC 41 ms
53,940 KB
testcase_12 AC 40 ms
54,984 KB
testcase_13 AC 41 ms
55,116 KB
testcase_14 AC 154 ms
77,676 KB
testcase_15 AC 54 ms
64,856 KB
testcase_16 AC 194 ms
77,720 KB
testcase_17 AC 189 ms
77,944 KB
testcase_18 AC 192 ms
77,856 KB
testcase_19 AC 192 ms
78,596 KB
testcase_20 AC 97 ms
76,604 KB
testcase_21 AC 69 ms
73,912 KB
testcase_22 AC 100 ms
76,540 KB
testcase_23 AC 95 ms
76,420 KB
testcase_24 AC 63 ms
84,252 KB
testcase_25 AC 64 ms
84,152 KB
testcase_26 AC 63 ms
84,176 KB
testcase_27 AC 253 ms
115,164 KB
testcase_28 AC 225 ms
115,012 KB
testcase_29 AC 260 ms
107,520 KB
testcase_30 AC 265 ms
106,724 KB
testcase_31 AC 295 ms
115,668 KB
testcase_32 AC 267 ms
115,808 KB
testcase_33 AC 482 ms
110,224 KB
testcase_34 AC 467 ms
110,432 KB
testcase_35 AC 494 ms
110,132 KB
testcase_36 AC 492 ms
110,616 KB
testcase_37 AC 474 ms
110,424 KB
testcase_38 AC 488 ms
110,288 KB
testcase_39 AC 472 ms
110,468 KB
testcase_40 AC 248 ms
103,168 KB
testcase_41 AC 186 ms
86,924 KB
testcase_42 AC 317 ms
93,920 KB
testcase_43 AC 115 ms
95,128 KB
testcase_44 AC 399 ms
100,712 KB
testcase_45 AC 220 ms
87,208 KB
testcase_46 AC 324 ms
102,980 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