結果
問題 |
No.2873 Kendall's Tau
|
ユーザー |
![]() |
提出日時 | 2024-07-06 12:18:22 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 337 ms / 4,500 ms |
コード長 | 1,638 bytes |
コンパイル時間 | 3,436 ms |
コンパイル使用メモリ | 264,760 KB |
実行使用メモリ | 18,144 KB |
最終ジャッジ日時 | 2024-07-06 12:19:31 |
合計ジャッジ時間 | 9,902 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/fenwicktree> using namespace std; using namespace atcoder; using ll = long long; int main() { int N; cin >> N; vector<pair<int, int>> p(N); for (int i = 0; i < N; i++) cin >> p[i].first >> p[i].second; // x座標を圧縮 vector<int> xs; for (int i = 0; i < N; i++) xs.push_back(p[i].first); sort(xs.begin(), xs.end()); xs.erase(unique(xs.begin(), xs.end()), xs.end()); for (int i = 0; i < N; i++) p[i].first = lower_bound(xs.begin(), xs.end(), p[i].first) - xs.begin(); // y座標を圧縮 vector<int> ys; for (int i = 0; i < N; i++) ys.push_back(p[i].second); sort(ys.begin(), ys.end()); ys.erase(unique(ys.begin(), ys.end()), ys.end()); for (int i = 0; i < N; i++) p[i].second = lower_bound(ys.begin(), ys.end(), p[i].second) - ys.begin(); // x座標でソート sort(p.begin(), p.end()); // x毎にy座標を管理 vector<vector<int>> v(xs.size()); for (int i = 0; i < N; i++) v[p[i].first].push_back(p[i].second); // BITでP,Q,R,Sを求める fenwick_tree<ll> fw(ys.size()); ll P, Q, R, S; P = Q = 0, R = S = (ll)N * (N - 1) / 2; for (int i = 0; i < xs.size(); i++) { for (auto y: v[i]) { P += fw.sum(0, y); Q += fw.sum(y + 1, ys.size()); } for (auto y: v[i]) fw.add(y, 1); R -= (ll)v[i].size() * (v[i].size() - 1) / 2; } for (int i = 0; i < ys.size(); i++) { ll cnt = fw.sum(i, i+1); S -= cnt * (cnt - 1) / 2; } cout << fixed << setprecision(16) << (double)(P-Q)/sqrt(R)/sqrt(S) << endl; }