結果
問題 | No.2197 Same Dish |
ユーザー |
![]() |
提出日時 | 2023-02-26 23:45:56 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 124 ms / 2,000 ms |
コード長 | 3,375 bytes |
コンパイル時間 | 912 ms |
コンパイル使用メモリ | 59,988 KB |
実行使用メモリ | 7,808 KB |
最終ジャッジ日時 | 2024-09-14 12:17:14 |
合計ジャッジ時間 | 2,926 ms |
ジャッジサーバーID (参考情報) |
judge6 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
/* -*- coding: utf-8 -*-** 2197.cc: No.2197 Same Dish - yukicoder*/#include<cstdio>#include<vector>#include<algorithm>#include<utility>using namespace std;/* constant */const int MAX_N = 100000;const int MAX_X = 200000;const int MOD = 998244353;/* typedef */typedef pair<int,int> pii;template<const int MOD>struct MI {int v;MI(): v() {}MI(int _v): v(_v % MOD) {}MI(long long _v): v(_v % MOD) {}MI operator+(const MI m) const { return MI(v + m.v); }MI operator-(const MI m) const { return MI(v + MOD - m.v); }MI operator*(const MI m) const { return MI((long long)v * m.v); }MI &operator+=(const MI m) { return (*this = *this + m); }MI &operator-=(const MI m) { return (*this = *this - m); }MI &operator*=(const MI m) { return (*this = *this * m); }bool operator==(const MI m) const { return v == m.v; }bool operator!=(const MI m) const { return v != m.v; }MI pow(int n) const { // a^n % MODMI pm = 1, a = *this;while (n > 0) {if (n & 1) pm *= a;a *= a;n >>= 1;}return pm;}MI inv() const { return pow(MOD - 2); }MI operator/(const MI m) const { return *this * m.inv(); }MI &operator/=(const MI m) { return (*this = *this / m); }};typedef MI<MOD> mi;template <typename T>struct SegTreeMaxDelay {int e2;vector<T> nodes, dls;T defv;SegTreeMaxDelay() {}void init(int n, T _defv) {defv = _defv;for (e2 = 1; e2 < n; e2 <<= 1);nodes.assign(e2 * 2, defv);dls.assign(e2 * 2, 0);}T &geti(int i) { return nodes[e2 - 1 + i]; }void seti(int i, T v) { geti(i) = v; }void setall() {for (int j = e2 - 2; j >= 0; j--)nodes[j] = max(nodes[j * 2 + 1], nodes[j * 2 + 2]);}void __update(int k) {if (dls[k] != 0) {int k0 = k * 2 + 1, k1 = k0 + 1;nodes[k0] += dls[k];nodes[k1] += dls[k];dls[k0] += dls[k];dls[k1] += dls[k];dls[k] = 0;}}void add_range(int r0, int r1, T v, int k, int i0, int i1) {if (r1 <= i0 || i1 <= r0) return;if (r0 <= i0 && i1 <= r1) {nodes[k] += v;dls[k] += v;return;}__update(k);int im = (i0 + i1) / 2;int k0 = k * 2 + 1, k1 = k0 + 1;add_range(r0, r1, v, k0, i0, im);add_range(r0, r1, v, k1, im, i1);nodes[k] = max(nodes[k0], nodes[k1]);}void add_range(int r0, int r1, T v) { add_range(r0, r1, v, 0, 0, e2); }T max_range(int r0, int r1, int k, int i0, int i1) {if (r1 <= i0 || i1 <= r0) return defv;if (r0 <= i0 && i1 <= r1) return nodes[k];__update(k);int im = (i0 + i1) / 2;T v0 = max_range(r0, r1, k * 2 + 1, i0, im);T v1 = max_range(r0, r1, k * 2 + 2, im, i1);return max(v0, v1);}T max_range(int r0, int r1) { return max_range(r0, r1, 0, 0, e2); }};/* global variables */pii ps[MAX_N];SegTreeMaxDelay<int> st;/* subroutines *//* main */int main() {int n, k;scanf("%d%d", &n, &k);for (int i = 0; i < n; i++) {int li, ri;scanf("%d%d", &li, &ri);ps[i] = pii(li, ri);}sort(ps, ps + n);st.init(MAX_X, 0);mi p = 1;for (int i = 0; i < n; i++) {int ri = ps[i].second, li = ps[i].first;int v = st.max_range(li, ri);p *= max(0, k - v);st.add_range(li, ri, 1);}mi s = mi(k).pow(n) - p;printf("%d\n", s.v);return 0;}