結果

問題 No.2230 Good Omen of White Lotus
コンテスト
ユーザー otoshigo
提出日時 2023-02-24 23:10:04
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 364 ms / 2,000 ms
コード長 2,201 bytes
コンパイル時間 1,355 ms
コンパイル使用メモリ 100,644 KB
最終ジャッジ日時 2025-02-10 22:23:57
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 44
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<vector>
#include<algorithm>
#include<map>
#include<atcoder/modint>
using namespace std;
using ll=long long;
using mint=atcoder::modint998244353;
#define rep(i,n) for(int i=0;i<n;i++)
#define rrep(i,n) for(int i=(n)-1;i>=0;i--)
#define all(v) v.begin(),v.end()
#define rall(v) v.rbegin(),v.rend()
template<class T> bool chmax(T &a, T b){if (a < b){a = b;return true;} else return false;}
template<class T> bool chmin(T &a, T b){if (a > b){a = b;return true;} else return false;}

template <class S, S (*op)(S, S), S (*e)()>
struct segtree{
    int n;
    vector<S> dat;
    segtree(int N){
        int x = 1;
        while (x < N) x *= 2;
        n = x;
        dat.resize(2 * n - 1, e());
    }
    void set(int i, S x){
        assert(i >= 0 && i < n);
        i += n - 1;
        dat[i] = x;
        while (i > 0){
            i = (i - 1) / 2;
            dat[i] = op(dat[2 * i + 1], dat[2 * i + 2]);
        }
    }
    S get(int i){
        assert(i >= 0 && i < n);
        i += n - 1;
        return dat[i];
    }
    S prod(int a, int b, int k = 0, int l = 0, int r = -1){
        if (r < 0) r = n;
        if (r <= a || l >= b) return e();
        if (l >= a && r <= b) return dat[k];
        else{
            S vl = prod(a, b, 2 * k + 1, l, (l + r) / 2);
            S vr = prod(a, b, 2 * k + 2, (l + r) / 2, r);
            return op(vl, vr);
        }
    }
};

int op(int l,int r){
    return max(l,r);
}

int e(){
    return 0;
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    int H,W,N;
    ll P;
    cin>>H>>W>>N>>P;
    vector<pair<int,int>>vp;
    vector<int>Y;
    rep(i,N){
        int x,y;
        cin>>x>>y;
        Y.push_back(y);
        vp.push_back(make_pair(x,y));
    }
    sort(all(vp));
    sort(all(Y));
    Y.erase(unique(all(Y)),Y.end());
    int l=Y.size();
    map<int,int>mp;
    rep(i,l)mp[Y[i]]=i;
    segtree<int,op,e>seg(l);
    int M=0;
    for(auto[x,y]:vp){
        int i=mp[y],m=seg.prod(0,i+1);
        seg.set(i,max(seg.get(i),m+1));
        chmax(M,max(seg.get(i),m+1));
    }
    mint ans=1;
    rep(i,M)ans*=1-2/(mint)P;
    rep(i,H+W-3-M)ans*=1-1/(mint)P;
    ans=1-ans;
    cout<<ans.val()<<"\n";
}
0