結果

問題 No.2873 Kendall's Tau
ユーザー ooaiu
提出日時 2026-11-12 17:27:35
言語 cLay
(20250308-1 + boost 1.92.0 + ACL)
コンパイル:
clayc _filename_
実行:
./a.out
結果
WA  
実行時間 -
コード長 571 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 6,660 ms
コンパイル使用メモリ 257,204 KB
実行使用メモリ 30,504 KB
最終ジャッジ日時 2026-11-12 17:27:50
合計ジャッジ時間 10,345 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge4_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14 WA * 18
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

ll@N,@(X,Y)[N];
coordcomp(N,X);
coordcomp(N,Y);
vector<ll>Z[N];
fenwick<int>fw;
fenwick<int>gw;
ll P = 0,Q = 0;
rep(i,N)Z[X[i]].push_back(Y[i]);
@[fw,gw].malloc(N,1);
rep(i,N){
for(int y:Z[i]){
P+=fw.range(0,y-1);
Q+=fw.range(y+1,N);
}
for(int y:Z[i]){
fw.add(y,1);
}
}
rrep(i,N){
for(int y:Z[i]){
P+=gw.range(y+1,N);
Q+=gw.range(0,y-1);
}
for(int y:Z[i]){
gw.add(y,1);
}
}
ll R=N*(N-1)/2,S=N*(N-1)/2;
ll cx[N]{},cy[N]{};
rep(i,N)cx[X[i]]++,cy[Y[i]]++;
rep(i,N)R-=cx[i]*(cx[i]-1)/2,S-=cy[i]*(cy[i]-1)/2;
double numer=P-Q;
double denom=std::sqrt(R*S);
wt(numer/denom/2.);
0