結果
問題 |
No.2230 Good Omen of White Lotus
|
ユーザー |
![]() |
提出日時 | 2022-12-08 20:49:35 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 462 ms / 2,000 ms |
コード長 | 1,657 bytes |
コンパイル時間 | 3,453 ms |
コンパイル使用メモリ | 167,880 KB |
最終ジャッジ日時 | 2025-02-09 06:29:21 |
ジャッジサーバーID (参考情報) |
judge4 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 44 |
ソースコード
#include <iostream> #include <vector> #include <map> #include <atcoder/all> using namespace std; using namespace atcoder; typedef long long ll; ll my_pow(ll x, ll n, ll mod){ ll ret; if (n == 0){ ret = 1; } else if (n % 2 == 1){ ret = (x * my_pow((x * x) % mod, n / 2, mod)) % mod; } else{ ret = my_pow((x * x) % mod, n / 2, mod); } return ret; } ll inv(ll x, ll mod){ return my_pow(x, mod - 2, mod); } ll op(ll a, ll b){ return max(a, b); } ll e(){ return -1e9; } ll mod2 = 998244353; int main(){ ll H,W,N,P; cin >> H >> W >> N >> P; vector<ll> x(N); vector<ll> y(N); for (ll i = 0; i < N; i++){ cin >> x[i] >> y[i]; } vector<pair<ll,ll>> p(N); for (ll i = 0; i < N; i++){ p[i].first = x[i]; p[i].second = y[i]; } sort(p.begin(), p.end()); map<ll, ll> mp; mp[1]++; for (ll i = 0; i < N; i++){ mp[y[i]]++; } ll now = 1; for (auto itr = mp.begin(); itr != mp.end(); itr++){ itr->second = now; now++; } segtree<ll, op, e> seg(now); seg.set(mp[1], 0); for (ll i = 0; i < N; i++){ ll prd = seg.prod(0, mp[p[i].second]); ll val = seg.get(mp[p[i].second]); seg.set(mp[p[i].second], max(val, prd) + 1); } ll max_omen = seg.prod(0, now); //cout << max_omen << endl; ll ans = (my_pow(P, H + W - 3, mod2)); ans = (ans - ((my_pow(P - 1, H + W - 3 - max_omen, mod2)) * (my_pow(P - 2, max_omen, mod2))) % mod2 + mod2) % mod2; ans = (ans * inv((my_pow(P, H + W - 3, mod2)), mod2)) % mod2; cout << ans << endl; }