結果

問題 No.2230 Good Omen of White Lotus
ユーザー shobonvipshobonvip
提出日時 2023-02-24 22:32:07
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 726 ms / 2,000 ms
コード長 1,926 bytes
コンパイル時間 253 ms
コンパイル使用メモリ 82,396 KB
実行使用メモリ 180,532 KB
最終ジャッジ日時 2024-09-13 05:46:37
合計ジャッジ時間 14,102 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
54,168 KB
testcase_01 AC 39 ms
53,556 KB
testcase_02 AC 37 ms
53,668 KB
testcase_03 AC 37 ms
54,284 KB
testcase_04 AC 37 ms
53,288 KB
testcase_05 AC 37 ms
54,068 KB
testcase_06 AC 37 ms
52,592 KB
testcase_07 AC 38 ms
53,964 KB
testcase_08 AC 38 ms
52,956 KB
testcase_09 AC 37 ms
53,328 KB
testcase_10 AC 37 ms
54,376 KB
testcase_11 AC 37 ms
53,364 KB
testcase_12 AC 39 ms
54,148 KB
testcase_13 AC 42 ms
56,068 KB
testcase_14 AC 220 ms
82,336 KB
testcase_15 AC 60 ms
67,068 KB
testcase_16 AC 288 ms
85,036 KB
testcase_17 AC 286 ms
85,172 KB
testcase_18 AC 292 ms
84,908 KB
testcase_19 AC 286 ms
85,040 KB
testcase_20 AC 113 ms
77,376 KB
testcase_21 AC 94 ms
76,820 KB
testcase_22 AC 115 ms
77,568 KB
testcase_23 AC 117 ms
77,696 KB
testcase_24 AC 49 ms
65,724 KB
testcase_25 AC 50 ms
66,276 KB
testcase_26 AC 50 ms
66,472 KB
testcase_27 AC 290 ms
105,204 KB
testcase_28 AC 295 ms
104,692 KB
testcase_29 AC 398 ms
172,416 KB
testcase_30 AC 394 ms
138,764 KB
testcase_31 AC 436 ms
180,404 KB
testcase_32 AC 430 ms
180,532 KB
testcase_33 AC 681 ms
144,600 KB
testcase_34 AC 675 ms
144,992 KB
testcase_35 AC 676 ms
145,864 KB
testcase_36 AC 714 ms
144,596 KB
testcase_37 AC 726 ms
144,856 KB
testcase_38 AC 678 ms
145,248 KB
testcase_39 AC 683 ms
144,724 KB
testcase_40 AC 324 ms
99,472 KB
testcase_41 AC 290 ms
94,896 KB
testcase_42 AC 439 ms
105,872 KB
testcase_43 AC 141 ms
79,688 KB
testcase_44 AC 605 ms
128,020 KB
testcase_45 AC 263 ms
99,276 KB
testcase_46 AC 423 ms
118,292 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

class SegmentTree:
	def __init__(self, n):
		tmp = 0
		while n != 0:
			n >>= 1
			tmp += 1
		self.seglen = 1 << tmp
		self.inf = 0
		self.seg = [self.inf]*(self.seglen*2)
	
	def update(self, x, val):
		# update x to val;
		x += self.seglen
		self.seg[x] = val
		while x > 1:
			x >>= 1
			self.seg[x] = max(self.seg[2*x], self.seg[2*x+1])
	
	def allupdate(self, a):
		targ = self.seglen
		for i in range(len(a)):
			self.seg[targ + i] = a[i]
		for i in range(self.seglen-1, -1, -1):
			self.seg[i] = max(self.seg[2*i], self.seg[2*i+1])

	def query(self, l, r):
		l = max(l, 0) + self.seglen
		r = max(r, 0) + self.seglen
		lres, rres = self.inf, self.inf
		while l < r:
			if l & 1:
				lres = max(lres, self.seg[l])
				l += 1
			if r & 1:
				r -= 1
				rres = max(rres, self.seg[r])
			l >>= 1
			r >>= 1
		return max(lres, rres)

def compress(arr):
	*XS, = set(arr)
	XS.sort()
	return {cmp_e: cmp_i for cmp_i, cmp_e in enumerate(XS)}

h,w,n,p = map(int,input().split())
ys = []
bucket = [[] for i in range(h)]
for i in range(n):
	x, y = map(int,input().split())
	bucket[x-1].append(y-1)
	ys.append(y-1)
y_c = compress(ys)
m = len(y_c)

seg = SegmentTree(m)
for i in range(h):
	bucket[i].sort()
	for j in bucket[i]:
		targ = y_c[j]
		ast = seg.query(0, targ+1) + 1
		ar = seg.query(targ, targ + 1)
		seg.update(targ, max(ar, ast))

tmp = seg.query(0, m)

mod = 998244353
yes = (1 - 2 * pow(p, mod-2, mod)) % mod
no = (1 - pow(p, mod-2, mod)) % mod

def best_fraction(n:int, mod:int, strmode=0) -> tuple:
	w = 9 * 10 ** 18
	t = pow(n, mod-2, mod)
	ansx, ansy = 0, 0
	for b in range(1, mod):
		if b*b > w: break
		a = n*b % mod
		c = t*b % mod
		if a*b < w:
			w = a*b
			ansx, ansy = a, b
		if c*b < w:
			w = c*b
			ansx, ansy = b, c
	if strmode:
		return f"{n} = {ansx}/{ansy} (mod {mod})"
	return (ansx, ansy)


ans = (1 - pow(yes, tmp, mod) * pow(no, h+w-3-tmp, mod)) % mod
print(ans)
#print(best_fraction(ans, mod, 1))
0