結果

問題 No.2230 Good Omen of White Lotus
ユーザー tree_splaytree_splay
提出日時 2023-02-24 22:51:17
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 72 ms / 2,000 ms
コード長 1,218 bytes
コンパイル時間 2,303 ms
コンパイル使用メモリ 209,336 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-13 05:54:29
合計ジャッジ時間 5,072 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 44
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

#define INF 1234567890
#define ll long long
#define MOD 998244353

ll pw(ll a, ll n)
{
	ll ret = 1;
	while(n)
	{
		if (n&1) ret=ret*a%MOD;
		a=a*a%MOD;
		n >>= 1;
	}
	return ret;
}

ll inv(ll a)
{
	return pw(a, MOD-2);
}

int N, M, K, P;
int fen[200001];

int Max(int idx)
{
	int ret = 0;
	while(idx)
	{
		ret = max(ret, fen[idx]);
		idx &= idx-1;
	}
	return ret;
}

void Add(int idx, int val)
{
	while(idx <= 200000)
	{
		fen[idx] = max(fen[idx], val);
		idx += idx & -idx;
	}
}

int main()
{
	ios::sync_with_stdio(0); cin.tie(0);
	cin.exceptions(ios::badbit | ios::failbit);
	cin >> N >> M >> K >> P;
	vector<pair<int, int> > v;
	for(int i=0; i<K; i++)
	{
		int y, x;
		cin >> y >> x;
		v.push_back({y, x});
	}
	v.push_back({N, M});
	sort(v.begin(), v.end());

	for(auto [y,x] : v)
	{
		int dp = Max(x) + 1;
		Add(x, dp);

		if (y == N && x == M)
		{
			int cnt = dp - 1;
			int dnt = N+M-1 - cnt - 2;

			ll cval = 2 * inv(P) % MOD;
			ll dval = inv(P);

			ll res = 1;
			while(cnt--) res *= (1-cval), res %= MOD;
			while(dnt--) res *= (1-dval), res %= MOD;
			res = 1-res;

			res %= MOD; res += MOD; res %= MOD;
			cout << res << "\n";
			break;
		}
	}
	return 0;
}
0