結果

問題 No.2873 Kendall's Tau
ユーザー sasa8uyauya
提出日時 2024-09-06 22:15:29
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 2,200 ms / 4,500 ms
コード長 964 bytes
コンパイル時間 268 ms
コンパイル使用メモリ 82,396 KB
実行使用メモリ 244,308 KB
最終ジャッジ日時 2024-09-06 22:16:18
合計ジャッジ時間 34,274 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

n=int(input())
p=[tuple(map(int,input().split())) for i in range(n)]

P=0
B=448
st1=[0]*B*B
st2=[0]*B
z=sorted(set([y for x,y in p]))
d={v:i for i,v in enumerate(z)}
p.sort(key=lambda a:(-a[0],a[1]))
for x,y in p:
  l=d[y]
  for i in range(l+1,((l+1)//B)*B+B):
    P+=st1[i]
  for i in range((l+1)//B+1,B):
    P+=st2[i]
  st1[l]+=1
  st2[l//B]+=1

Q=0
B=448
st1=[0]*B*B
st2=[0]*B
z=sorted(set([y for x,y in p]))
d={v:i for i,v in enumerate(z)}
p.sort(key=lambda a:(-a[0],-a[1]))
for x,y in p:
  l=d[y]
  for i in range(((l-1)//B)*B,l-1+1):
    Q+=st1[i]
  for i in range((l-1)//B):
    Q+=st2[i]
  st1[l]+=1
  st2[l//B]+=1

R=0
z=sorted(set([x for x,y in p]))
d={v:i for i,v in enumerate(z)}
c=[0]*len(z)
for x,y in p:
  c[d[x]]+=1
for i in range(len(z)):
  R+=c[i]*(n-c[i])
R//=2

S=0
z=sorted(set([y for x,y in p]))
d={v:i for i,v in enumerate(z)}
c=[0]*len(z)
for x,y in p:
  c[d[y]]+=1
for i in range(len(z)):
  S+=c[i]*(n-c[i])
S//=2

print((P-Q)/(R*S)**0.5)
0