結果
問題 | No.2873 Kendall's Tau |
ユーザー | titia |
提出日時 | 2024-09-09 02:55:51 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 2,157 ms / 4,500 ms |
コード長 | 5,499 bytes |
コンパイル時間 | 388 ms |
コンパイル使用メモリ | 82,248 KB |
実行使用メモリ | 189,800 KB |
最終ジャッジ日時 | 2024-09-09 02:56:28 |
合計ジャッジ時間 | 35,727 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 65 ms
68,508 KB |
testcase_01 | AC | 64 ms
69,856 KB |
testcase_02 | AC | 64 ms
68,188 KB |
testcase_03 | AC | 65 ms
68,056 KB |
testcase_04 | AC | 64 ms
68,228 KB |
testcase_05 | AC | 64 ms
68,656 KB |
testcase_06 | AC | 65 ms
68,536 KB |
testcase_07 | AC | 2,140 ms
169,252 KB |
testcase_08 | AC | 2,157 ms
174,724 KB |
testcase_09 | AC | 2,094 ms
168,900 KB |
testcase_10 | AC | 2,119 ms
173,032 KB |
testcase_11 | AC | 2,099 ms
165,808 KB |
testcase_12 | AC | 2,103 ms
189,800 KB |
testcase_13 | AC | 726 ms
107,704 KB |
testcase_14 | AC | 1,632 ms
164,176 KB |
testcase_15 | AC | 422 ms
93,000 KB |
testcase_16 | AC | 446 ms
92,704 KB |
testcase_17 | AC | 1,636 ms
153,396 KB |
testcase_18 | AC | 1,160 ms
136,168 KB |
testcase_19 | AC | 1,583 ms
149,924 KB |
testcase_20 | AC | 425 ms
91,540 KB |
testcase_21 | AC | 1,289 ms
138,580 KB |
testcase_22 | AC | 580 ms
100,892 KB |
testcase_23 | AC | 1,306 ms
136,108 KB |
testcase_24 | AC | 280 ms
83,028 KB |
testcase_25 | AC | 408 ms
91,112 KB |
testcase_26 | AC | 1,613 ms
151,448 KB |
testcase_27 | AC | 983 ms
121,876 KB |
testcase_28 | AC | 1,825 ms
168,632 KB |
testcase_29 | AC | 1,883 ms
175,692 KB |
testcase_30 | AC | 357 ms
87,752 KB |
testcase_31 | AC | 555 ms
95,980 KB |
testcase_32 | AC | 1,276 ms
134,704 KB |
ソースコード
import sys input = sys.stdin.readline from collections import Counter from operator import itemgetter # https://github.com/tatyam-prime/SortedSet/blob/main/SortedMultiset.py import math from bisect import bisect_left, bisect_right from typing import Generic, Iterable, Iterator, List, Tuple, TypeVar, Optional T = TypeVar('T') class SortedMultiset(Generic[T]): BUCKET_RATIO = 16 SPLIT_RATIO = 24 def __init__(self, a: Iterable[T] = []) -> None: "Make a new SortedMultiset from iterable. / O(N) if sorted / O(N log N)" a = list(a) n = self.size = len(a) if any(a[i] > a[i + 1] for i in range(n - 1)): a.sort() num_bucket = int(math.ceil(math.sqrt(n / self.BUCKET_RATIO))) self.a = [a[n * i // num_bucket : n * (i + 1) // num_bucket] for i in range(num_bucket)] def __iter__(self) -> Iterator[T]: for i in self.a: for j in i: yield j def __reversed__(self) -> Iterator[T]: for i in reversed(self.a): for j in reversed(i): yield j def __eq__(self, other) -> bool: return list(self) == list(other) def __len__(self) -> int: return self.size def __repr__(self) -> str: return "SortedMultiset" + str(self.a) def __str__(self) -> str: s = str(list(self)) return "{" + s[1 : len(s) - 1] + "}" def _position(self, x: T) -> Tuple[List[T], int, int]: "return the bucket, index of the bucket and position in which x should be. self must not be empty." for i, a in enumerate(self.a): if x <= a[-1]: break return (a, i, bisect_left(a, x)) def __contains__(self, x: T) -> bool: if self.size == 0: return False a, _, i = self._position(x) return i != len(a) and a[i] == x def count(self, x: T) -> int: "Count the number of x." return self.index_right(x) - self.index(x) def add(self, x: T) -> None: "Add an element. / O(√N)" if self.size == 0: self.a = [[x]] self.size = 1 return a, b, i = self._position(x) a.insert(i, x) self.size += 1 if len(a) > len(self.a) * self.SPLIT_RATIO: mid = len(a) >> 1 self.a[b:b+1] = [a[:mid], a[mid:]] def _pop(self, a: List[T], b: int, i: int) -> T: ans = a.pop(i) self.size -= 1 if not a: del self.a[b] return ans def discard(self, x: T) -> bool: "Remove an element and return True if removed. / O(√N)" if self.size == 0: return False a, b, i = self._position(x) if i == len(a) or a[i] != x: return False self._pop(a, b, i) return True def lt(self, x: T) -> Optional[T]: "Find the largest element < x, or None if it doesn't exist." for a in reversed(self.a): if a[0] < x: return a[bisect_left(a, x) - 1] def le(self, x: T) -> Optional[T]: "Find the largest element <= x, or None if it doesn't exist." for a in reversed(self.a): if a[0] <= x: return a[bisect_right(a, x) - 1] def gt(self, x: T) -> Optional[T]: "Find the smallest element > x, or None if it doesn't exist." for a in self.a: if a[-1] > x: return a[bisect_right(a, x)] def ge(self, x: T) -> Optional[T]: "Find the smallest element >= x, or None if it doesn't exist." for a in self.a: if a[-1] >= x: return a[bisect_left(a, x)] def __getitem__(self, i: int) -> T: "Return the i-th element." if i < 0: for a in reversed(self.a): i += len(a) if i >= 0: return a[i] else: for a in self.a: if i < len(a): return a[i] i -= len(a) raise IndexError def pop(self, i: int = -1) -> T: "Pop and return the i-th element." if i < 0: for b, a in enumerate(reversed(self.a)): i += len(a) if i >= 0: return self._pop(a, ~b, i) else: for b, a in enumerate(self.a): if i < len(a): return self._pop(a, b, i) i -= len(a) raise IndexError def index(self, x: T) -> int: "Count the number of elements < x." ans = 0 for a in self.a: if a[-1] >= x: return ans + bisect_left(a, x) ans += len(a) return ans def index_right(self, x: T) -> int: "Count the number of elements <= x." ans = 0 for a in self.a: if a[-1] > x: return ans + bisect_right(a, x) ans += len(a) return ans N=int(input()) X=[list(map(int,input().split())) for i in range(N)] P=0 Q=0 R=0 S=0 X.sort(key=itemgetter(1),reverse=True) X.sort(key=itemgetter(0)) PS=SortedMultiset() for x,y in X: #print(P,Q) P+=PS.index(y) #Q+=len(PS)-PS.index(y+1) PS.add(y) X.sort(key=itemgetter(1)) X.sort(key=itemgetter(0)) PS=SortedMultiset() for x,y in X: #print(P,Q) #P+=PS.index(y) Q+=len(PS)-PS.index(y+1) PS.add(y) #print(P,Q) C=Counter() C2=Counter() for x,y in X: C[x]+=1 C2[y]+=1 for c in C: R+=C[c]*(N-C[c]) for c2 in C2: S+=C2[c2]*(N-C2[c2]) #print(P,Q,R//2,S//2) print((P-Q)/(R*S//4)**(1/2))