結果

問題 No.2873 Kendall's Tau
ユーザー katonyonkokatonyonko
提出日時 2024-09-11 13:53:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,541 ms / 4,500 ms
コード長 2,387 bytes
コンパイル時間 313 ms
コンパイル使用メモリ 82,128 KB
実行使用メモリ 230,920 KB
最終ジャッジ日時 2024-09-11 13:53:33
合計ジャッジ時間 26,655 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 43 ms
55,948 KB
testcase_01 AC 43 ms
55,252 KB
testcase_02 AC 44 ms
56,696 KB
testcase_03 AC 43 ms
56,192 KB
testcase_04 AC 44 ms
56,672 KB
testcase_05 AC 45 ms
57,192 KB
testcase_06 AC 43 ms
56,604 KB
testcase_07 AC 1,428 ms
188,964 KB
testcase_08 AC 1,532 ms
219,164 KB
testcase_09 AC 1,438 ms
189,308 KB
testcase_10 AC 1,541 ms
230,920 KB
testcase_11 AC 1,442 ms
188,844 KB
testcase_12 AC 1,482 ms
226,936 KB
testcase_13 AC 474 ms
110,768 KB
testcase_14 AC 1,197 ms
204,412 KB
testcase_15 AC 296 ms
104,840 KB
testcase_16 AC 337 ms
101,536 KB
testcase_17 AC 1,039 ms
166,344 KB
testcase_18 AC 875 ms
161,032 KB
testcase_19 AC 1,050 ms
162,008 KB
testcase_20 AC 311 ms
105,132 KB
testcase_21 AC 881 ms
149,212 KB
testcase_22 AC 416 ms
114,668 KB
testcase_23 AC 926 ms
148,436 KB
testcase_24 AC 184 ms
85,828 KB
testcase_25 AC 281 ms
101,892 KB
testcase_26 AC 1,045 ms
160,248 KB
testcase_27 AC 687 ms
144,364 KB
testcase_28 AC 1,251 ms
175,024 KB
testcase_29 AC 1,339 ms
209,172 KB
testcase_30 AC 240 ms
90,460 KB
testcase_31 AC 394 ms
108,696 KB
testcase_32 AC 858 ms
149,972 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import io
import sys
from collections import defaultdict, deque, Counter
from itertools import permutations, combinations, accumulate
from heapq import heappush, heappop
from bisect import bisect_right, bisect_left
from math import gcd
import math

_INPUT = """\
6
3
1 4
1 5
9 2
2
1 1
2 2
5
4 4
5 7
0 4
2 10
-6 4
"""

class BIT:
    def __init__(self, n):
        self._n = n
        self.data = [0] * n
    def add(self, p, x):
        assert 0 <= p < self._n
        p += 1
        while p <= self._n:
            self.data[p - 1] += x
            p += p & -p
    #合計にはrを含まない
    def sum(self, l, r):
        assert 0 <= l <= r <= self._n
        return self._sum(r) - self._sum(l)
    def _sum(self, r):
        s = 0
        while r > 0:
            s += self.data[r - 1]
            r -= r & -r
        return s
    #pの位置をxという値にセット
    def set(self, p, x):
        self.add(p, -self.sum(p, p+1) + x)

def input():
  return sys.stdin.readline()[:-1]

def solve(test):
  N=int(input())
  data=[list(map(int, input().split())) for _ in range(N)]
  d={x:i for i,x in enumerate(sorted(set([x for x,y in data]+[y for x,y in data])))}
  data=[(d[x],d[y]) for x,y in data]
  data.sort(key=lambda x:(x[0],x[1]))
  bit=BIT(len(d))
  Q=0
  for x,y in data:
    Q+=bit.sum(y+1,len(d))
    bit.add(y,1)
  bit=BIT(len(d))
  data.sort(key=lambda x:(x[0],-x[1]))
  # print(data)
  P=0
  for x,y in data:
    P+=bit.sum(0,y)
    bit.add(y,1)
  C1=Counter([x for x,y in data])
  C2=Counter([y for x,y in data])
  R=sum([(N-C1[x])*C1[x] for x in C1])//2
  S=sum([(N-C2[y])*C2[y] for y in C2])//2
  print((P-Q)/(R*S)**0.5)

def random_input():
  from random import randint,shuffle
  N=randint(1,10)
  M=randint(1,N)
  A=list(range(1,M+1))+[randint(1,M) for _ in range(N-M)]
  shuffle(A)
  return (" ".join(map(str, [N,M]))+"\n"+" ".join(map(str, A))+"\n")*3

def simple_solve():
  return []

def main(test):
  if test==0:
    solve(0)
  elif test==1:
    sys.stdin = io.StringIO(_INPUT)
    case_no=int(input())
    for _ in range(case_no):
      solve(0)
  else:
    for i in range(1000):
      sys.stdin = io.StringIO(random_input())
      x=solve(1)
      y=simple_solve()
      if x!=y:
        print(i,x,y)
        print(*[line for line in sys.stdin],sep='')
        break

#0:提出用、1:与えられたテスト用、2:ストレステスト用
main(0)
0