結果
問題 | No.2230 Good Omen of White Lotus |
ユーザー | jiangly |
提出日時 | 2023-02-24 21:54:14 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 69 ms / 2,000 ms |
コード長 | 3,985 bytes |
コンパイル時間 | 2,230 ms |
コンパイル使用メモリ | 208,308 KB |
実行使用メモリ | 5,504 KB |
最終ジャッジ日時 | 2023-10-11 06:10:25 |
合計ジャッジ時間 | 6,355 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,348 KB |
testcase_01 | AC | 2 ms
4,348 KB |
testcase_02 | AC | 2 ms
4,348 KB |
testcase_03 | AC | 1 ms
4,348 KB |
testcase_04 | AC | 2 ms
4,348 KB |
testcase_05 | AC | 1 ms
4,348 KB |
testcase_06 | AC | 2 ms
4,352 KB |
testcase_07 | AC | 2 ms
4,352 KB |
testcase_08 | AC | 2 ms
4,348 KB |
testcase_09 | AC | 2 ms
4,352 KB |
testcase_10 | AC | 1 ms
4,352 KB |
testcase_11 | AC | 2 ms
4,348 KB |
testcase_12 | AC | 1 ms
4,348 KB |
testcase_13 | AC | 2 ms
4,348 KB |
testcase_14 | AC | 25 ms
4,348 KB |
testcase_15 | AC | 2 ms
4,352 KB |
testcase_16 | AC | 36 ms
4,672 KB |
testcase_17 | AC | 36 ms
4,536 KB |
testcase_18 | AC | 36 ms
4,800 KB |
testcase_19 | AC | 36 ms
4,692 KB |
testcase_20 | AC | 4 ms
4,348 KB |
testcase_21 | AC | 2 ms
4,348 KB |
testcase_22 | AC | 4 ms
4,352 KB |
testcase_23 | AC | 5 ms
4,348 KB |
testcase_24 | AC | 2 ms
4,348 KB |
testcase_25 | AC | 2 ms
4,352 KB |
testcase_26 | AC | 2 ms
4,352 KB |
testcase_27 | AC | 41 ms
5,404 KB |
testcase_28 | AC | 43 ms
5,376 KB |
testcase_29 | AC | 43 ms
5,328 KB |
testcase_30 | AC | 49 ms
5,428 KB |
testcase_31 | AC | 45 ms
5,320 KB |
testcase_32 | AC | 43 ms
5,328 KB |
testcase_33 | AC | 68 ms
5,480 KB |
testcase_34 | AC | 68 ms
5,448 KB |
testcase_35 | AC | 69 ms
5,316 KB |
testcase_36 | AC | 68 ms
5,504 KB |
testcase_37 | AC | 68 ms
5,424 KB |
testcase_38 | AC | 67 ms
5,376 KB |
testcase_39 | AC | 68 ms
5,464 KB |
testcase_40 | AC | 21 ms
4,348 KB |
testcase_41 | AC | 21 ms
4,352 KB |
testcase_42 | AC | 44 ms
4,356 KB |
testcase_43 | AC | 5 ms
4,348 KB |
testcase_44 | AC | 57 ms
4,764 KB |
testcase_45 | AC | 23 ms
4,352 KB |
testcase_46 | AC | 36 ms
4,400 KB |
ソースコード
#include <bits/stdc++.h> using i64 = long long; template<class T> constexpr T power(T a, i64 b) { T res = 1; for (; b; b /= 2, a *= a) { if (b % 2) { res *= a; } } return res; } template<int P> struct MInt { int x; constexpr MInt() : x{} {} constexpr MInt(i64 x) : x{norm(x % P)} {} constexpr int norm(int x) const { if (x < 0) { x += P; } if (x >= P) { x -= P; } return x; } constexpr int val() const { return x; } explicit constexpr operator int() const { return x; } constexpr MInt operator-() const { MInt res; res.x = norm(P - x); return res; } constexpr MInt inv() const { assert(x != 0); return power(*this, P - 2); } constexpr MInt &operator*=(MInt rhs) { x = 1LL * x * rhs.x % P; return *this; } constexpr MInt &operator+=(MInt rhs) { x = norm(x + rhs.x); return *this; } constexpr MInt &operator-=(MInt rhs) { x = norm(x - rhs.x); return *this; } constexpr MInt &operator/=(MInt rhs) { return *this *= rhs.inv(); } friend constexpr MInt operator*(MInt lhs, MInt rhs) { MInt res = lhs; res *= rhs; return res; } friend constexpr MInt operator+(MInt lhs, MInt rhs) { MInt res = lhs; res += rhs; return res; } friend constexpr MInt operator-(MInt lhs, MInt rhs) { MInt res = lhs; res -= rhs; return res; } friend constexpr MInt operator/(MInt lhs, MInt rhs) { MInt res = lhs; res /= rhs; return res; } friend constexpr std::istream &operator>>(std::istream &is, MInt &a) { i64 v; is >> v; a = MInt(v); return is; } friend constexpr std::ostream &operator<<(std::ostream &os, const MInt &a) { return os << a.val(); } friend constexpr bool operator==(MInt lhs, MInt rhs) { return lhs.val() == rhs.val(); } friend constexpr bool operator!=(MInt lhs, MInt rhs) { return lhs.val() != rhs.val(); } }; template<int V, int P> constexpr MInt<P> CInv = MInt<P>(V).inv(); constexpr int P = 998244353; using Z = MInt<P>; template <typename T> struct Fenwick { int n; std::vector<T> a; Fenwick(int n = 0) { init(n); } void init(int n) { this->n = n; a.assign(n, T()); } void add(int x, T v) { for (int i = x + 1; i <= n; i += i & -i) { a[i - 1] += v; } } T sum(int x) { auto ans = T(); for (int i = x; i > 0; i -= i & -i) { ans += a[i - 1]; } return ans; } T rangeSum(int l, int r) { return sum(r) - sum(l); } int kth(T k) { int x = 0; for (int i = 1 << std::__lg(n); i; i /= 2) { if (x + i <= n && k >= a[x + i - 1]) { x += i; k -= a[x - 1]; } } return x; } }; struct Max { int x = 0; Max &operator+=(Max a) & { if (a.x > x) { x = a.x; } return *this; } }; int main() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); int H, W, N; Z p; std::cin >> H >> W >> N >> p; p = 1 / p; std::vector<std::pair<int, int>> a(N); for (int i = 0; i < N; i++) { int x, y; std::cin >> x >> y; x--, y--; a[i] = {x, y}; } std::sort(a.begin(), a.end()); Fenwick<Max> fen(W); int cnt = 0; for (auto [x, y] : a) { int f = 1 + fen.sum(y + 1).x; fen.add(y, {f}); cnt = std::max(cnt, f); } Z ans = 1 - power(1 - 2 * p, cnt) * power(1 - p, H + W - 3 - cnt); std::cout << ans << "\n"; return 0; }