#1008057 (C++17) No.2873 Kendall's Tau

• 問題: No.2873 Kendall's Tau
• ユーザー: eve__fuyuki
• 提出日時: 2024-09-07 01:33:30
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 213 ms / 4,500 ms
• コード長: 1,557 bytes
• コンパイル時間: 2,656 ms
• コンパイル使用メモリ: 218,428 KB
• 実行使用メモリ: 9,240 KB
• 最終ジャッジ日時: 2024-09-07 01:33:38
• 合計ジャッジ時間: 7,501 ms
• ジャッジサーバーID （参考情報）: judge5 / judge4

入力	結果	実行時間; 実行使用メモリ
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,944 KB
testcase_07	AC	213 ms; 8,324 KB
testcase_08	AC	196 ms; 9,240 KB
testcase_09	AC	189 ms; 8,188 KB
testcase_10	AC	193 ms; 9,092 KB
testcase_11	AC	187 ms; 8,320 KB
testcase_12	AC	188 ms; 8,820 KB
testcase_13	AC	75 ms; 6,940 KB
testcase_14	AC	152 ms; 8,180 KB
testcase_15	AC	32 ms; 6,940 KB
testcase_16	AC	34 ms; 6,940 KB
testcase_17	AC	136 ms; 7,028 KB
testcase_18	AC	100 ms; 6,944 KB
testcase_19	AC	130 ms; 6,976 KB
testcase_20	AC	33 ms; 6,944 KB
testcase_21	AC	114 ms; 6,944 KB
testcase_22	AC	48 ms; 6,944 KB
testcase_23	AC	112 ms; 6,940 KB
testcase_24	AC	15 ms; 6,940 KB
testcase_25	AC	30 ms; 6,940 KB
testcase_26	AC	163 ms; 7,028 KB
testcase_27	AC	85 ms; 6,940 KB
testcase_28	AC	157 ms; 7,932 KB
testcase_29	AC	169 ms; 8,468 KB
testcase_30	AC	24 ms; 6,940 KB
testcase_31	AC	45 ms; 6,940 KB
testcase_32	AC	111 ms; 6,940 KB

ソースコード

#include <bits/stdc++.h>
using namespace std;

void fast_io() {
	ios_base::sync_with_stdio(false);
	cin.tie(nullptr);
}

#include <atcoder/fenwicktree>
using namespace atcoder;
int main() {
	fast_io();
	int n;
	cin >> n;
	vector<pair<int, int>> xy(n);
	vector<int> vals;
	for (int i = 0; i < n; i++) {
		cin >> xy[i].first >> xy[i].second;
		vals.push_back(xy[i].first);
		vals.push_back(xy[i].second);
	}
	sort(vals.begin(), vals.end());
	vals.erase(unique(vals.begin(), vals.end()), vals.end());
	int m = vals.size();
	for (int i = 0; i < n; i++) {
		xy[i].first =
			lower_bound(vals.begin(), vals.end(), xy[i].first) - vals.begin();
		xy[i].second =
			lower_bound(vals.begin(), vals.end(), xy[i].second) - vals.begin();
	}
	long long r = (long long)n * (n - 1) / 2, s = r;
	{
		vector<int> cnt_x(m), cnt_y(m);
		for (auto [x, y] : xy) {
			r -= cnt_x[x];
			cnt_x[x]++;
			s -= cnt_y[y];
			cnt_y[y]++;
		}
	}
	long long p = 0;
	{
		sort(xy.begin(), xy.end(), [&](auto a, auto b) {
			if (a.first == b.first) return a.second > b.second;
			return a.first < b.first;
		});
		fenwick_tree<long long> fw(m);
		for (auto [_, y] : xy) {
			p += fw.sum(0, y);
			fw.add(y, 1);
		}
	}
	long long q = 0;
	{
		sort(xy.begin(), xy.end(), [&](auto a, auto b) {
			if (a.first == b.first) return a.second < b.second;
			return a.first < b.first;
		});
		fenwick_tree<long long> fw(m);
		for (auto [_, y] : xy) {
			q += fw.sum(y + 1, m);
			fw.add(y, 1);
		}
	}
	long double ans = (p - q) / sqrtl(r) / sqrtl(s);
	cout << fixed << setprecision(16) << ans << endl;
}