hanoy_binary.py

def restore_disks(size, action_index, action_count, from_axis, to_axis):
    pivot = action_count // 2
    free_axis = _select_axis(from_axis, to_axis)
    part_count = action_count - pivot - 1

    if action_index == pivot:
        axises[from_axis].append(size)

        for i in range(size-1, 0, -1):
            axises[free_axis].append(i)

    elif action_index < pivot:
        axises[from_axis].append(size)
        restore_disks(size-1, action_index, part_count, from_axis, free_axis)
    else:
        action_index -= pivot + 1
        axises[to_axis].append(size)
        restore_disks(size-1, action_index, part_count, free_axis, to_axis)


def restore_towers(index):
    global axises
    rings_count = 6
    axises = [[], [], []]

    restore_disks(rings_count, index, 2 ** 6 - 1, 0, 2)
    print(axises)


for i in range(0, 2 ** 6 - 1):
    restore_towers(i)