結果

問題 No.2873 Kendall's Tau
ユーザー n0ma_run0ma_ru
提出日時 2024-09-07 01:40:15
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 325 ms / 4,500 ms
コード長 2,153 bytes
コンパイル時間 2,658 ms
コンパイル使用メモリ 214,712 KB
実行使用メモリ 18,720 KB
最終ジャッジ日時 2024-09-07 01:40:24
合計ジャッジ時間 9,216 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 2 ms
6,940 KB
testcase_06 AC 2 ms
6,948 KB
testcase_07 AC 325 ms
13,952 KB
testcase_08 AC 316 ms
17,388 KB
testcase_09 AC 302 ms
13,696 KB
testcase_10 AC 324 ms
18,720 KB
testcase_11 AC 325 ms
13,440 KB
testcase_12 AC 301 ms
16,384 KB
testcase_13 AC 93 ms
6,944 KB
testcase_14 AC 244 ms
15,104 KB
testcase_15 AC 51 ms
6,944 KB
testcase_16 AC 53 ms
6,940 KB
testcase_17 AC 219 ms
11,648 KB
testcase_18 AC 165 ms
10,752 KB
testcase_19 AC 215 ms
12,416 KB
testcase_20 AC 53 ms
6,940 KB
testcase_21 AC 176 ms
9,856 KB
testcase_22 AC 75 ms
6,944 KB
testcase_23 AC 181 ms
9,984 KB
testcase_24 AC 25 ms
6,944 KB
testcase_25 AC 49 ms
6,940 KB
testcase_26 AC 249 ms
10,752 KB
testcase_27 AC 137 ms
9,856 KB
testcase_28 AC 253 ms
15,228 KB
testcase_29 AC 264 ms
13,440 KB
testcase_30 AC 40 ms
6,944 KB
testcase_31 AC 74 ms
6,944 KB
testcase_32 AC 180 ms
10,880 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define ALL(v) v.begin(),v.end()
#define dbg(x) cerr << #x << ": " << (x) << endl;
template<class F, class S>
ostream& operator<<(ostream& os, pair<F,S>& p) {
    os << '(' << p.first << ',' << p.second << ')';
    return os;
}
template<class Iter>
void print(Iter beg, Iter end) {
    for (Iter itr = beg; itr != end; ++itr) {
        cerr << *itr << ' ';
    }
    cerr << '\n';
}
struct BIT {
    int n;
    vector<ll> d;
    BIT (int n) : n(n), d(n+1) {}

    void add(int idx, ll val) {
        ++idx;
        for (int x = idx; x <= n; x += x&-x) {
            d[x] += val;
        }
    }
    ll sum(int idx) {
        ll res = 0;
        for (int x = idx; x > 0; x -= x&-x) {
            res += d[x];
        }
        return res;
    }
    ll sum(int l, int r) {
        return sum(r) - sum(l);
    }
};
int n;
vector<int> x,y;
int main() {
    cin >> n;
    x.resize(n); y.resize(n);
    for (int i = 0; i < n; ++i) cin >> x[i] >> y[i];

    auto tmp = x;
    sort(ALL(tmp));
    tmp.erase(unique(ALL(tmp)), tmp.end());
    for (int i = 0; i < n; ++i) x[i] = lower_bound(ALL(tmp), x[i]) - tmp.begin();
    int xmax = tmp.size();
    tmp = y;
    sort(ALL(tmp));
    tmp.erase(unique(ALL(tmp)), tmp.end());
    for (int i = 0; i < n; ++i) y[i] = lower_bound(ALL(tmp), y[i]) - tmp.begin();
    int ymax = tmp.size();

    ll p=0, q=0, r=0, s=0;
    BIT bit(ymax);
    vector<vector<int>> ys(xmax);
    for (int i = 0; i < n; ++i) ys[x[i]].push_back(y[i]);
    for (int i = 0; i < xmax; ++i) {
        for (int j : ys[i]) {
            p += bit.sum(j);
            q += bit.sum(j+1, ymax);
        }
        for (int j : ys[i]) {
            bit.add(j, 1);
        }
    }
    vector<int> mx(xmax),my(ymax);
    for (int i = 0; i < n; ++i) {
        mx[x[i]]++;
        my[y[i]]++;
    }
    for (int i = 0; i < xmax; ++i) {
        r += (ll)mx[i] * (n - mx[i]);
    }
    for (int i = 0; i < ymax; ++i) {
        s += (ll)my[i] * (n - my[i]);
    }
    r >>= 1;
    s >>= 1;
    double ans = (p - q) / (sqrt(r) * sqrt(s));
    cout << fixed << setprecision(20) << ans << '\n';
}
0