結果

問題 No.2873 Kendall's Tau
ユーザー hongrockhongrock
提出日時 2024-09-06 22:15:11
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 511 ms / 4,500 ms
コード長 1,276 bytes
コンパイル時間 2,393 ms
コンパイル使用メモリ 217,704 KB
実行使用メモリ 40,392 KB
最終ジャッジ日時 2024-09-06 22:15:36
合計ジャッジ時間 9,763 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 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,944 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 375 ms
27,384 KB
testcase_08 AC 438 ms
35,652 KB
testcase_09 AC 377 ms
26,916 KB
testcase_10 AC 511 ms
40,392 KB
testcase_11 AC 393 ms
26,016 KB
testcase_12 AC 446 ms
33,528 KB
testcase_13 AC 86 ms
11,008 KB
testcase_14 AC 341 ms
31,184 KB
testcase_15 AC 55 ms
9,728 KB
testcase_16 AC 51 ms
8,832 KB
testcase_17 AC 269 ms
22,388 KB
testcase_18 AC 220 ms
22,032 KB
testcase_19 AC 320 ms
23,792 KB
testcase_20 AC 52 ms
8,960 KB
testcase_21 AC 214 ms
18,348 KB
testcase_22 AC 75 ms
10,516 KB
testcase_23 AC 205 ms
19,232 KB
testcase_24 AC 22 ms
6,940 KB
testcase_25 AC 49 ms
9,084 KB
testcase_26 AC 254 ms
19,448 KB
testcase_27 AC 165 ms
18,180 KB
testcase_28 AC 362 ms
30,876 KB
testcase_29 AC 387 ms
27,388 KB
testcase_30 AC 34 ms
6,940 KB
testcase_31 AC 67 ms
10,092 KB
testcase_32 AC 236 ms
20,960 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

#define pb emplace_back
#define mp make_pair
 
using ll = long long;
using pii = pair<int,int>;
 
constexpr int mod = 998244353;
constexpr int inf = 0x3f3f3f3f;
constexpr int N = 2e5 + 10;

int n, m, b[N];
int s[N];

int lowbit(int x){
  return x & (-x);
}
void add(int x, int v){
  for(; x<=m; x+=lowbit(x)) s[x] += v;
}
int sum(int x){
  int ret = 0;
  for(; x>0; x-=lowbit(x))  ret += s[x];
  return ret;
}

int query(int L, int R){
  if(L > R) return 0;
  return sum(R) - sum(L - 1);
}

void _main(){
  unordered_map<int, int> cntx, cnty;
  map<int, vector<int>> g;
  cin >> n;
  ll P = 0, Q = 0, R = 0, S = 0;
  int x, y;
  m = 0;
  for(int i=1; i<=n; ++i){
    cin >> x >> y;
    R += (i - ++cntx[x]);
    S += (i - ++cnty[y]);
    g[x].pb(y);
    b[m++] = y;
  }
  sort(b, b+m);
  m = unique(b, b+m) - b;
  for(auto &[x, vec] : g){
    for(int y : vec){
      int j = lower_bound(b, b+m, y) - b + 1;
      P += query(1, j - 1);
      Q += query(j + 1, m);
    }
    for(int y : vec){
      int j = lower_bound(b, b+m, y) - b + 1;
      add(j, 1);
    }
  }
  cout << setprecision(8) << (P - Q) * 1.0 / sqrt(1.0 * R * S) << '\n';
}

int main(){
  ios::sync_with_stdio(0);
  cin.tie(0); cout.tie(0);
  _main();
  return 0;
}
0