Problem solution · Python

ABC425 D — Ulam-Warburton Automaton

ABC425 D — Ulam-Warburton Automaton: a Python solution using breadth-first search. Learn the idea, check the complexity, and read the full code, with credit to KATO-Hiro AtCoder Solutions.

Technique
Breadth-first search
Source
KATO-Hiro AtCoder Solutions
Length
70 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

Breadth-first search

For ABC425 D — Ulam-Warburton Automaton, the implementation explores reachable states in layers, which is the standard shape for unweighted shortest paths and minimum-step transitions.

  1. Model each valid configuration as a state and each legal move as an edge.
  2. Seed the queue with the starting state and mark it immediately.
  3. Expand each state once, recording distance or reachability for unseen neighbours.

Code notes

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

Complexity

Verify that each state and transition is processed only a bounded number of times; that determines the traversal cost.

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.

Full codeABC425 D — Ulam-Warburton Automaton · PythonPython
Use this to learn the idea, then write your own version.
# -*- coding: utf-8 -*-  def main():    import sys    from collections import deque     input = sys.stdin.readline     h, w = map(int, input().split())    s = [list(input().rstrip()) for _ in range(h)]     q = deque()    inf = 10**18    dist = [[inf for _ in range(w)] for _ in range(h)]     for i in range(h):        for j in range(w):            if s[i][j] == ".":                continue             q.append((i, j))            dist[i][j] = 0     dxy = [(1, 0), (-1, 0), (0, 1), (0, -1)]     while q:        y, x = q.popleft()        di = dist[y][x]         for dx, dy in dxy:            ny, nx = y + dy, x + dx             if not (0 <= ny < h):                continue            if not (0 <= nx < w):                continue            if dist[ny][nx] != inf:                continue             count = 0             for dx2, dy2 in dxy:                ny2, nx2 = ny + dy2, nx + dx2                 if not (0 <= ny2 < h):                    continue                if not (0 <= nx2 < w):                    continue                if dist[ny2][nx2] <= di:                    count += 1                    continue             if count == 1:                q.append((ny, nx))                dist[ny][nx] = dist[y][x] + 1     ans = 0     for i in range(h):        for j in range(w):            if dist[i][j] != inf:                ans += 1     print(ans)  if __name__ == "__main__":    main() 

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 ↗