結果
問題 | No.2336 Do you like typical problems? |
ユーザー | srjywrdnprkt |
提出日時 | 2024-10-20 22:43:22 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 197 ms / 2,000 ms |
コード長 | 1,997 bytes |
コンパイル時間 | 2,375 ms |
コンパイル使用メモリ | 213,748 KB |
実行使用メモリ | 14,128 KB |
最終ジャッジ日時 | 2024-10-20 22:43:28 |
合計ジャッジ時間 | 6,085 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,816 KB |
testcase_01 | AC | 2 ms
6,820 KB |
testcase_02 | AC | 2 ms
6,820 KB |
testcase_03 | AC | 2 ms
6,816 KB |
testcase_04 | AC | 2 ms
6,816 KB |
testcase_05 | AC | 2 ms
6,816 KB |
testcase_06 | AC | 2 ms
6,820 KB |
testcase_07 | AC | 2 ms
6,816 KB |
testcase_08 | AC | 3 ms
6,816 KB |
testcase_09 | AC | 3 ms
6,816 KB |
testcase_10 | AC | 3 ms
6,820 KB |
testcase_11 | AC | 3 ms
6,816 KB |
testcase_12 | AC | 3 ms
6,816 KB |
testcase_13 | AC | 194 ms
13,468 KB |
testcase_14 | AC | 197 ms
14,128 KB |
testcase_15 | AC | 194 ms
13,404 KB |
testcase_16 | AC | 196 ms
13,468 KB |
testcase_17 | AC | 195 ms
13,328 KB |
testcase_18 | AC | 67 ms
12,848 KB |
testcase_19 | AC | 77 ms
12,080 KB |
testcase_20 | AC | 176 ms
13,704 KB |
ソースコード
#include <bits/stdc++.h> #include <atcoder/modint> using namespace std; using ll = long long; using namespace atcoder; using mint = modint998244353; int main(){ cin.tie(nullptr); ios_base::sync_with_stdio(false); /* N!通りの並べ方のうち、a_i != a_j となるものについて a_i > a_j となるパターンは N!/2 通り よって、全パターンの和は1つのパターンについて a_i != a_j となる(i<j)の数の期待値のN!/2倍 L_i = B_i, R_i = C_i と固定する。 a_i != a_j となる(i<j)の期待値は 全ての(i<j)の個数N(N-1)/2からa_i = a_jとなる(i<j)の個数の期待値を 引いたもの a_i = xとなる確率をW_iとすると、a_i = a_j = xとなる確率は sum_{i<j} W_iW_j = ((sum_{i} W_i)^2-sum_{i} W_i^2)/2 これをすべてのxについて足しあわせる。 座標圧縮してimos法 */ ll N, M; cin >> N; vector<ll> B(N), C(N), X; vector<mint> W(N); for (int i=0; i<N; i++){ cin >> B[i] >> C[i]; C[i]++; X.push_back(B[i]); X.push_back(C[i]); W[i] = mint(C[i]-B[i]).inv(); } sort(X.begin(), X.end()); X.erase(unique(X.begin(), X.end()), X.end()); for (int i=0; i<N; i++){ B[i] = lower_bound(X.begin(), X.end(), B[i]) - X.begin(); C[i] = lower_bound(X.begin(), X.end(), C[i]) - X.begin(); } M = X.size(); vector<mint> V1(M+1), V2(M+1); for (int i=0; i<N; i++){ V1[B[i]] += W[i]; V1[C[i]] -= W[i]; V2[B[i]] += W[i]*W[i]; V2[C[i]] -= W[i]*W[i]; } for (int i=1; i<=M; i++){ V1[i] += V1[i-1]; V2[i] += V2[i-1]; } mint ans=N*(N-1)/2, iv=mint(2).inv(); for (int i=0; i<M; i++){ ans -= (V1[i]*V1[i]-V2[i]) * iv * (X[i+1]-X[i]); } for (int i=1; i<=N; i++){ ans *= i; } ans *= iv; cout << ans.val() << endl; return 0; }