CPMagical CP Library

Problem solution · Python

ABC241 D — Sequence Query

ABC241 D — Sequence Query: a Python solution using sorting and greedy selection. Learn the idea, check the complexity, and read the full code, with credit to KATO-Hiro AtCoder Solutions.

Technique
Sorting and greedy selection
Source
KATO-Hiro AtCoder Solutions
Length
175 lines
Start with the idea.

Try the problem first. If you get stuck, read the approach below, then write your own solution. The full code is at the bottom.

Approach

Sorting and greedy selection

For ABC241 D — Sequence Query, the implementation first exposes a useful order, then scans that order while making locally justified choices.

  1. Choose the key that reveals the greedy or grouping structure.
  2. Sort the relevant records by that key.
  3. Scan in order, maintaining the invariant that makes each local choice safe.

Code notes

  • 175 lines of Python from the credited upstream file abc241_d.py.
  • The implementation visibly relies on sequence storage, ordered lookup.
  • No explicit loop blocks detected.

Complexity

Sorting is typically the dominant term unless the subsequent scan uses a more expensive nested operation.

Check the problem constraints before deciding whether this complexity will pass.

Source

Code and credit

This code comes from KATO-Hiro AtCoder Solutions by KATO-Hiro and is used under the CC0-1.0 licence.

Open problem ↗Original source ↗
Full codeABC241 D — Sequence Query · PythonPython
Use this to learn the idea, then write your own version.
# -*- coding: utf-8 -*-  # https://github.com/tatyam-prime/SortedSet/blob/main/SortedMultiset.pyimport mathfrom bisect import bisect_left, bisect_right, insortfrom typing import Generic, Iterable, Iterator, TypeVar, Union, ListT = TypeVar('T') class SortedMultiset(Generic[T]):    BUCKET_RATIO = 50    REBUILD_RATIO = 170     def _build(self, a=None) -> None:        "Evenly divide `a` into buckets."        if a is None: a = list(self)        size = self.size = len(a)        bucket_size = int(math.ceil(math.sqrt(size / self.BUCKET_RATIO)))        self.a = [a[size * i // bucket_size : size * (i + 1) // bucket_size] for i in range(bucket_size)]        def __init__(self, a: Iterable[T] = []) -> None:        "Make a new SortedMultiset from iterable. / O(N) if sorted / O(N log N)"        a = list(a)        if not all(a[i] <= a[i + 1] for i in range(len(a) - 1)):            a = sorted(a)        self._build(a)     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 __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 _find_bucket(self, x: T) -> List[T]:        "Find the bucket which should contain x. self must not be empty."        for a in self.a:            if x <= a[-1]: return a        return a     def __contains__(self, x: T) -> bool:        if self.size == 0: return False        a = self._find_bucket(x)        i = bisect_left(a, 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 = self._find_bucket(x)        insort(a, x)        self.size += 1        if len(a) > len(self.a) * self.REBUILD_RATIO:            self._build()     def discard(self, x: T) -> bool:        "Remove an element and return True if removed. / O(√N)"        if self.size == 0: return False        a = self._find_bucket(x)        i = bisect_left(a, x)        if i == len(a) or a[i] != x: return False        a.pop(i)        self.size -= 1        if len(a) == 0: self._build()        return True     def lt(self, x: T) -> Union[T, None]:        "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) -> Union[T, None]:        "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) -> Union[T, None]:        "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) -> Union[T, None]:        "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, x: int) -> T:        "Return the x-th element, or IndexError if it doesn't exist."        if x < 0: x += self.size        if x < 0: raise IndexError        for a in self.a:            if x < len(a): return a[x]            x -= 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  def main():    import sys     input = sys.stdin.readline     q = int(input())    s = SortedMultiset()    ans = list()     for _ in range(q):        qi = list(map(int, input().split()))        x = qi[1]         if qi[0] == 1:            s.add(x)        elif qi[0] == 2:            k = qi[2]            count = s.index_right(x)             if count >= k:                ans.append(s[count - k])            else:                ans.append(-1)        else:            k = qi[2]            count = s.index(x)            remain = len(s) - count             if remain >= k:                ans.append(s[count + k - 1])            else:                ans.append(-1)        print('\n'.join(map(str, ans)))  if __name__ == "__main__":    main()

Keep practising

Did this explanation save you time? I'm a Grade 11 student building this free library to make difficult algorithms easier to understand.

Buy me a coffee ↗