🚲 Tanguy Ortolo a écrit 12197 commentaires

À un moment il faut donner un numéro de téléphone, mais ça ne change rien si c'est celui de quelqu'un d'autre.

Tiens, je n'avais pas pensé à Pypy. Pour info, ma solution termine la seconde partie – ligne par ligne, avec calcul d'unions d'intervalles – en 2,7 s sous Pypy.

J'y ai réfléchi cette nuit, on peut faire encore mieux. Premièrement parce qu'il n'est pas vraiment nécessaire de considérer les intersections de droites, mais seulement celles des segments qui délimitent les zones de chaque capteur (enfin, les segments juste un point plus loin en fait).

Mais surtout parce que le point cherché est soit :

• dans un coin du terrain et sur un segment ;
• sur un bord du terrain et à l'intersection de deux segments ;
• n'importe où ailleurs, mais à l'intersection de quatre segments.

Bref, in n'a plus qu'à tester si un coin ne se trouverait pas sur un segment, et sinon, si on construit la liste des intersections, on n'a plus qu'à chercher si on ne trouve pas dedans :

• un point qui serait sur un bord ;
• un point qui s'y trouverait quatre fois.

Je ne réponds pas à un numéro qui n'est pas dans mon annuaire perso (même si c'est un livreur, un commerçant ou un artisan qui m'appelle pour un truc important, figurez-vous qu'il va me laisser un mot sur le répondeur !!! dingue non ?)

Quand on a une messagerie vocale, oui. J'ai volontairement désactivé la mienne, et les fournisseurs divers risqueraient de moyennement apprécier le fait que je leur demande de laisser un SMS sans possibilité de laisser un message vocal…

Ça, c'est bien quand on supporte les messages vocaux. Personnellement, je déteste ça, et après des années d'hésitation, j'ai finalement désactivé ma messagerie vocale. J'ai juste un message d'absence qui invite à laisser un SMS. Le genre de truc qu'un démarcheur ne fera jamais depuis son poste fixe d'ailleurs.

Il me semble que c'est un brin plus nuancé. Je ne suis pas certain qu'il soit illégal d'enregistrer un appel qu'on reçoit dans le consentement de l'interlocuteur. Mais même si ça l'était, s'il ne l'apprend jamais il ne pourra pas s'en plaindre, en attendant c'est bien pratique pour vérifier le nom qu'à donné le démarcheur par exemple.

Rendre public cet enregistrement, ça pourrait peut-être bien être risqué en revanche. Et ça pourrait bien ne pas être recevable comme preuve en effet.

Farpaitement !

Et j'y ai aussi découvert les dataclasses. Très pratique ça, avec en plus l'option congélateur, ça va me simplifier la vie la prochaine fois que j'aurai affaire à des coordonnées.

Euh, ça dépend lesquels. J'ai affaire à une boîte particulièrement coriace, leur nom, en tout cas le nom qu'ils prétendent avoir, est en deux mots, ça commence par Habitat et ça finit par Fermeture. Reçu peut-être une dizaine d'appels, genre un par mois depuis un an. Ils sont ma prochaine cible.

Je viens de tomber sur une solution que je trouve assez géniale : https://github.com/JakubDotPy/aoc2022/blob/master/day15/part2.py

La solution sera forcément à une intersection de frontières de zones couvertes par les capteurs. Bon, pas forcément, elle pourrait être dans un coin, mais c'est plus que douteux. On pourrait quand même tester.

Bref, entre les intersections et les quatre coins de la grille, il n'est pas difficile de vérifier si l'un de ces points ne serait pas par hasard hors de portée de tous les capteurs…

Comment prouver qu'il y eu démarchage abusif ?

Plusieurs possibilités :

• Je suis sur Bloctel. Ce n'était pas un de mes fournisseurs, ni une association, ni un journal. L'appel était donc forcément abusif.
• C'est pour vendre une formation CPF : ce genre de démarchage est interdit tout court, donc forcément abusif.
• C'est pour vendre de la rénovation énergétique : ce genre de démarchage est interdit tout court, donc forcément abusif.

Pour ce qui est de prouver que j'ai été appelé, lorsque je fais un signalement, je laisse mes coordonnées, donc au besoin je peux fournir mon journal d'appels et un enregistrement de la conversation.

Et si la DGCCRF décide de lancer des poursuites, ils demanderont à la police ou à la gendarmerie d'enquêter, et mon opérateur pourra sans problème confirmer que j'ai bien reçu l'appel en question.

C'était quand? On en parlait récemment après un article; et après une recherche rapide, j'ai pu confirmer qu'il y a une loi qui interdit tout démarchage pour une formation CPF, adoptée ce 8 décembre, avec 375 000 euros d'amende pour une personne morale (donc l'organisme de formation, mais aussi le sous-traitant? Ou l'un d'eux paye pour les 2? Ou ils partagent l'amende?).

Peut-être les deux. Mais pour répondre à la question, c'était avant cette loi, dont je ne pouvais pas me prévaloir. Sauf que je suis sur Bloctel, donc sauf quelques domaines bien précis dont la formation professionnelle ne fait pas partie, il est de toute façon interdit de me démarcher pour me vendre des trucs.

Vu la réaction du type que j'avais en face, ils avaient l'air plutôt conscients d'être hors-la-loi, ce qui suggère qu'ils l'étaient bel et bien. D'ailleurs ça me rappelle un peu mieux la fin du dialogue. Je n'avais pas raccroché aussi sec :

• Vous n'étiez pas intéressé par la formation ?
• Non, désolé, j'avais juste besoin d'informations pour le signalement.
• Pour signaler quoi ?
• Eh bien, votre appel illégal.
• Mais ça n'a rien d'illégal Monsieur !
• Dites, vous avez entendu parler de Bloctel ?
• Oui, bien sûr, et nous vérifions les gens que nous appelons.
• Ben non. J'y suis, sur Bloctel justement, donc vous n'aviez pas le droit de m'appeler.
• Mais si notre appel vous dérangeait, vous n'aviez qu'à dire que vous n'étiez pas intéressé. Pourquoi avez-vous fait croire que ça vous intéressait, vous nous avez tendu un piège ?
• C'est ça. Ça vous a plu ?
• Mais c'est dégueulasse ce que vous faites !
• Ah non, moi j'aide à faire respecter la loi. Ce que vous faites en revanche, c'est illégal. Si j'ai besoin d'une formation, je suis assez grand pour aller la chercher tout seul. Allez, bonne journée.

Si ce n'est pas en vigueur, quelle est la base pour dire que c'est illicite? C'est parce que tu es sur Bloctel (le cas "automatiquement illicite" que tu cites)? Parce qu'hormis cela, il faudrait par exemple prouver qu'ils ont eu ton numéro de manière illicite (c'est à dire qu'un organisme leur a donné alors que tu avais explicité refusé, ou simplement rien accepté; mais comment prouver cela?).

Bloctel en effet. Mais pour aller plus loin, lorsqu'on fait l'objet d'un traitement de données, et le démarchage téléphonique inclut évidemment un traitement de données, au minimum les numéros de téléphone, on a le droit d'accéder à ces données, qui doivent inclure leur origine. Une boîte qui refuse de fournir les données qu'elle a sur celui qui le demande, ou qui refuse de préciser d'où ils les tiennent, est déjà en violation de la loi.

Sinon similairement à ǝpɐןƃu∀ nǝıɥʇʇɐW-ǝɹɹǝıԀ plus bas, il y a aussi la question de savoir s'ils auront vraiment une amende. Ils pourraient se retourner en disant qu'ils n'ont jamais appelé (c'est alors une parole contre une autre). Pour vérifier les dires, il faudrait qu'un juge acte pour avoir accès aux appels téléphoniques, mais le feraient-ils sur le signalement d'une personne seulement (je doute malheureusement que beaucoup signalent et au final, chaque organisme a probablement un ou 2 signalements à peine à leur actif)?

Alors pour la preuve, j'ai mon journal d'appel, et au besoin un enregistrement de la conversation. Ah oui, j'oubliais de le préciser, j'enregistre. Sans le dire, sinon ils raccrochent direct, donc c'est un peu limite comme preuve, mais la jurisprudence évolue un peu sur le sujet.

Pour les plaintes à propos de démarchage au CPF, ça peut aussi se faire directement sur cette plate-forme, et je suppose qu'ils en tiennent un peu compte quand même, vu le niveau d'arnaques qu'il y a eu. Mais autrement, c'est vrai, je n'ai aucune certitude que ce sera suivi d'effets. La seule certitude que j'ai, c'est que si personne ne signale les abus, ils ne seront jamais sanctionnés.

J'ai un dernier cas: des fois, on reçoit des appels et quand on répond, il n'y a rien. Personne parle. Ça m'arrive de temps en temps. Mais je ne saurais pas classer cela.

Ça vient de robots d'appels. Le principe, c'est que la boîte de démarchage a un vrai centre d'appel dédié à cela, et un robot qui appelle des gens à un rythme calculé selon la disponibilité des esclaves employés. Lorsque quelqu'un répond, il lui passe le premier troufion disponible, mais parfois il n'y en a pas, et l'appelé se retrouve en attente.

Il va sans dire que je ne suis pas pleinement satisfait de la recherche en semi-force brute pour la deuxième étape, intellectuellement parlant.

Mais bon, la flemme de concevoir et d'implémenter un algorithme d'union de parallélogrammes. Il y a des tas de cas, et autant avec des rectangles ça se fait bien en itérant sur chaque coordonnée, autant avec les parallélogrammes, c'est compliqué.

Voilà, avec une implémentation des unions d'intervalles sécants.

from __future__ import annotations

import re

from collections import namedtuple
from typing import Tuple, Iterable, List, Optional, Type

import aoc


Coords = Tuple[int, int]


class Point(namedtuple('Point', ['x', 'y'])):
    __slots__ = ()

    def __str__(self) -> str:
        return "({},{})".format(self.x, self.y)

    def dist(self, other: Point) -> int:
        return abs(self.x - other.x) + abs(self.y - other.y)


class Interval:
    def __init__(self, start: int, end: int) -> None:
        if start >= end:
            raise ValueError("unsupported empty or negative interval")
        self.start = start
        self.end = end

    def __str__(self) -> str:
        return "[{},{}[".format(self.start, self.end)

    def __len__(self) -> int:
        return self.end - self.start

    def __contains__(self, value: int) -> bool:
        return value >= self.start and value < self.end

    def intersects(self, other: Interval) -> bool:
        return max(self.start, other.start) < min(self.end, other.end)

    def union_update(self, other: Interval) -> bool:
        if not self.intersects(other):
            return False
        self.start = min(self.start, other.start)
        self.end = max(self.end, other.end)
        return True


class Sensor:
    def __init__(self, position: Point, beacon: Point):
        self.position = position
        self.beacon = beacon
        self.range = position.dist(beacon)

    def covers(self, y: int) -> bool:
        dist = abs(y - self.position.y)
        return dist <= self.range

    def coverage(self, y: int) -> Optional[Interval]:
        if not self.covers(y):
            return None
        else:
            dist = abs(y - self.position.y)
            return Interval(self.position.x - self.range + dist,
                            self.position.x + self.range + 1 - dist)

    re_line = re.compile(r'^Sensor at x=(-?\d+), y=(-?\d+): '
                         + r'closest beacon is at x=(-?\d+), y=(-?\d+)\n?$')

    @classmethod
    def import_line(class_, line: str) -> Sensor:
        match = class_.re_line.match(line)
        if match is None:
            raise ValueError("invalid sensor description")
        position = Point(int(match.group(1)), int(match.group(2)))
        beacon = Point(int(match.group(3)), int(match.group(4)))
        return class_(position, beacon)


class SensorArray:
    def __init__(self, sensors: Iterable[Sensor]):
        self.sensors = list(sensors)

    @classmethod
    def import_lines(class_, lines: Iterable[str]) -> SensorArray:
        return class_(Sensor.import_line(line) for line in lines)

    def coverage(self, y: int) -> List[Interval]:
        intervals: List[Interval] = []
        for sensor in self.sensors:
            if not sensor.covers(y):
                continue
            new_interval = sensor.coverage(y)
            if new_interval is None:
                continue
            new_intervals: List[Interval] = []
            for interval in intervals:
                if not new_interval.union_update(interval):
                    # new_interval did not absorb interval
                    new_intervals.append(interval)
            new_intervals.append(new_interval)
            intervals = new_intervals
        return intervals


def solve1(lines: Iterable[str]) -> int:
    """Solve part 1 of today's puzzle"""
    array = SensorArray.import_lines(lines)
    beacons = {sensor.beacon for sensor in array.sensors}
    y = 2000000
    coverage = array.coverage(y)
    total = 0
    for interval in coverage:
        total += len(interval)
        total -= sum((beacon.y == y and beacon.x in interval)
                     for beacon in beacons)
    return total


def solve2(lines: Iterable[str]) -> int:
    """Solve part 2 of today's puzzle"""
    array = SensorArray.import_lines(lines)
    position: Optional[Point] = None
    for y in range(4000000):
        if y % 100000 == 0:
            print(y)
        coverage = array.coverage(y)
        if len(coverage) < 2:
            continue
        return min(interval.end for interval in coverage) * 4000000 + y
    raise ValueError("no solution found")

Un détail tout de même, ça semble évident mais ça peut demander un peu de concentration pour ne pas l'oublier en court de route : tant que vous n'êtes pas arrivé à vos fins, c'est à dire que vous n'avez pas assez d'informations, il ne faut pas donner le moindre indice qui laisserait soupçonner que vous n'appréciez pas ce démarchage.

Éviter par exemple des phrases comme « Au fait, comment avez-vous eu mon numéro de téléphone ? » Si votre interlocuteur commence à soupçonner que son appel vous dérange, il risque de ne pas vous croire si vous vous prétendez intéressé, et l'appel ne durera pas longtemps. Pas assez longtemps pour obtenir les infos que vous souhaitez en tout cas.

Ces démarcheurs savent très bien que ce qu'ils font est illégal, même s'ils prétendent parfois le contraire. Et sachant qu'ils peuvent être poursuivis, ou avoir diverses emmerdes comme être radiés de la plate-forme CPF, ils commencent à être assez méfiants.

L'autre jour, j'ai oublié ce principe et demandé directement à une démarcheuse où elle travaillait. Elle m'a fait répéter la question, et après avoir compris ma demande, elle a raccroché sans y mettre la moindre forme. Ni refus, ni excuse, ni au-revoir, juste raccroché. Si ce n'est pas de la méfiance, ça…

À titre d'exemple, j'ai pu piéger une entreprise de formation professionnelle qui se livrait au démarchage téléphonique. Sans doute sous-traité, mais ce n'est pas mon problème, les agissements d'un sous-traitant engagent évidemment la responsabilité de son donneur d'ordre.

Je me suis montré intéressé par une formation en langue anglaise, très utile pour mon travail. Le démarcheur m'a alors fait rappeler par un responsable des inscriptions, et je m'attendais au pire, c'est à dire à ce qu'il me demande mon numéro de sécurité sociale et qu'il réinitialise avec mon aide mon mot de passe pour accéder lui-même à mon compte CPF et m'inscrire à sa formation. Ç'aurait été une arnaque selon les règles, mais non, c'était un peu plus honnête, il m'a juste demandé mon adresse électronique, celle inscrite sur le compte CPF, ce qui lui a suffit à me pré-inscrire à des formations.

Vérification faite, les formations apparaissaient bien comme des propositions à valider, avec un organisme de formation tout à fait identifiable, en tout cas par la plate-forme CPF. J'en avais assez, et la fin de la discussion téléphonique a été assez intéressante :

• C'est bon, je vois bien les formations. Je vous remercie, j'ai toutes les informations qu'il me faut.
• Attendez, je peux vous expliquer la façon dont les formations vont se passer.
• Ce ne sera pas nécessaire, vraiment, j'ai tout ce qu'il me faut.
• Non, je ne vous ai pas expliqué…
• Si si, j'en sais assez je vous dit.
• Assez pour quoi ?
• Assez pour signaler votre démarchage illicite.
• Vous n'étiez pas intéressé par la formation ?
• Non, désolé, j'avais juste besoin d'informations pour le signalement.
• Vous nous avez tendu un piège alors ?
• C'est ça. Ça vous a plu ?
• Mais c'est dégueulasse ce que vous faites !
• Ah non, moi j'aide à faire respecter la loi. Ce que vous faites en revanche, c'est illégal. Si j'ai besoin d'une formation, je suis assez grand pour aller la chercher tout seul. Allez, bonne journée.

Avec ladite optimisation :

from __future__ import annotations

import enum
import io
import itertools

from typing import Iterable, Iterator, List, Tuple

import numpy as np
import numpy.typing as npt


Coords = Tuple[int, int]


class Material(enum.Enum):
    AIR = enum.auto()
    ROCK = enum.auto()
    SAND = enum.auto()


class Slice:
    def __init__(self, array: npt.ArrayLike, leak: Coords) -> None:
        self.matrix: npt.NDArray = np.array(array)
        self.leak = leak

    def pour(self) -> Iterator[bool]:
        """Inject a unit of sand.

        Yields False if it was blocked, True if it fell into the void.
        Stops when the leak is clogged up by a mountain of sand."""
        # List of coordinates to pour sand from
        sources: List[Coords] = []
        # Start at the original leak
        y, x = self.leak
        while True:
            # Our unit of sand is injected at current coordinates (y, x)
            while y < self.matrix.shape[0] - 1:
                y_ = y + 1
                for x_ in (x, x - 1, x + 1):
                    if self.matrix[y_, x_] is Material.AIR:
                        # Sand can flow down:
                        # * current coordinate is a good place to pour next
                        #   sand unit from;
                        sources.append((y, x))
                        # * current sand unit falls down one level.
                        y = y_
                        x = x_
                        break
                else:
                    # Ways down exhausted, sand is blocked
                    self.matrix[y, x] = Material.SAND
                    # Get out of the descent loop
                    break
            else:
                # Descent ended, sand fell all the way down
                yield True
                continue
            # Descent did not end, sand got blocked
            yield False
            if not sources:
                # Nowhere to pour sand from, even the leak is clogged up
                break
            # Pour next sand unit from last good place
            y, x = sources.pop()

    def __str__(self) -> str:
        result = io.StringIO()
        for line in self.matrix:
            for point in line:
                if point is Material.AIR:
                    result.write(' ')
                elif point is Material.ROCK:
                    result.write('█')
                elif point is Material.SAND:
                    result.write('░')
                else:
                    assert False  # we covered all cases
            result.write('\n')
        return result.getvalue()

    def add_segment(self, start: Coords, end: Coords) -> None:
        y1, x1 = start
        y2, x2 = end
        if y1 == y2:
            x1, x2 = min(x1, x2), max(x1, x2)
            ys = (y1 for _ in range(x1, x2 + 1))
            xs = (x for x in range(x1, x2 + 1))
        elif x1 == x2:
            y1, y2 = min(y1, y2), max(y1, y2)
            ys = (y for y in range(y1, y2 + 1))
            xs = (x1 for _ in range(y1, y2 + 1))
        for y, x in zip(ys, xs):
            self.matrix[y, x] = Material.ROCK


def import_segments(
        lines: Iterable[str]
        ) -> Tuple[List[Coords], List[Tuple[Coords, Coords]]]:
    points: List[Coords] = []
    segments: List[Tuple[Coords, Coords]] = []

    def coords(word: str) -> Coords:
        parts = word.split(',')
        if len(parts) != 2:
            raise ValueError("unexpected number of coordinates")
        return int(parts[1]), int(parts[0])

    for line in lines:
        words = line.rstrip().split(' -> ')
        segment_points = [coords(word) for word in words]
        points.extend(segment_points)
        segments.extend(itertools.pairwise(segment_points))
    return points, segments


def import_slice1(lines: Iterable[str]) -> Slice:
    points, segments = import_segments(lines)
    y_leak, x_leak = 0, 500
    points.append((y_leak, x_leak))
    xmin = min(x for _, x in points) - 1
    xmax = max(x for _, x in points) + 1
    ymin = min(y for y, _ in points)
    ymax = max(y for y, _ in points)
    if ymin < 0:
        raise ValueError("unexpected negative coordinate")
    xshift = xmin
    matrix = np.full(((ymax + 1), xmax - xshift + 1), Material.AIR)
    result = Slice(matrix, (0, 500 - xshift))
    for (y1, x1), (y2, x2) in segments:
        result.add_segment((y1, x1 - xshift), (y2, x2 - xshift))
    return result


def import_slice2(lines: Iterable[str]) -> Slice:
    points, segments = import_segments(lines)
    y_leak, x_leak = 0, 500
    points.append((y_leak, x_leak))
    xmin = min(x for _, x in points)
    xmax = max(x for _, x in points)
    ymin = min(y for y, _ in points)
    ymax = max(y for y, _ in points) + 2  # including the floor
    # Update xmin and xmax to accomodate a pyramid of sand
    xmin = min(xmin, x_leak - (ymax - y_leak))
    xmax = max(xmax, x_leak + (ymax - y_leak))
    if ymin < 0:
        raise ValueError("unexpected negative coordinate")
    xshift = xmin
    matrix = np.full(((ymax + 1), xmax - xshift + 1), Material.AIR)
    result = Slice(matrix, (0, 500 - xshift))
    for (y1, x1), (y2, x2) in segments:
        result.add_segment((y1, x1 - xshift), (y2, x2 - xshift))
    result.add_segment((ymax, xmin - xshift), (ymax, xmax - xshift))
    return result


def solve1(lines: Iterable[str]) -> int:
    """Solve part 1 of today's puzzle"""
    slice_ = import_slice1(lines)
    for i, fallen in enumerate(slice_.pour()):
        if fallen:
            return i
    raise ValueError("Simulation ended with unexpected sand leak clogged")


def solve2(lines: Iterable[str]) -> int:
    """Solve part 2 of today's puzzle"""
    slice_ = import_slice2(lines)
    for i, fallen in enumerate(slice_.pour()):
        if fallen:
            raise ValueError("Simulation ended with unexpected sand under the floor")
    return i + 1

Bravo pour l'implémentation de mon idée !

Visiblement, il y a une optimisation possible : une si un grain de sable, lâché du point de départ, tombe à un endroit, on peut lâcher ce grain et les suivants de cet endroit, jusqu'à ce qu'il soit ensablé…

Je me demande s'il ne serait pas possible de résoudre la partie 2 bien plus vite, en calculant l'aire des zones protégées par des rochers.

En effet, le nombre d'unités de sables tombées, c'est l'aire d'une triangle plein, moins celle des rochers qu'il contient, moins l'aire des zones protégées par ces derniers.

Ceci dit, déterminer la zone protégée par l'ensemble des rochers, c'est loin d'être évident. La zone protégée par un segment horizontal, c'est facile, c'est un triangle dessous. La zone protégée par un segment vertical, c'est facile, c'est rien du tout. Mais quand on commence à avoir des segments horizontaux et verticaux qui se touchent, ça devient compliqué.

Allez, une autre piste, sans doute pas beaucoup plus simple. On dirait qu'on peut déterminer la zone protégée par des rochers, non pas seulement par calcul, mais aussi par une simulation bizarre : faire couler de l'air depuis le bas vers le haut et regarder où il s'accumule. Mais ça pose le problème d'introduire cet air : il faut qu'il s'accumule et qu'il traverse en même temps les segments…

from __future__ import annotations

import enum
import io
import itertools

from typing import Iterable, List, Tuple

import numpy as np
import numpy.typing as npt


Coords = Tuple[int, int]


class Material(enum.Enum):
    AIR = enum.auto()
    ROCK = enum.auto()
    SAND = enum.auto()


class Slice:
    def __init__(self, array: npt.ArrayLike, leak: Coords) -> None:
        self.matrix: npt.NDArray = np.array(array)
        self.leak = leak

    def pour(self) -> bool:
        """Inject a unit of sand.

        Return False if it was blocked, True if it fell into the void."""
        y, x = self.leak
        while y < self.matrix.shape[0] - 1:
            y_ = y + 1
            for x_ in (x, x - 1, x + 1):
                if self.matrix[y_, x_] is Material.AIR:
                    # Sand can flow down
                    y = y_
                    x = x_
                    break
            else:
                # Blocked
                self.matrix[y, x] = Material.SAND
                return False
        # Sand fell all the way down into the void
        return True

    def __str__(self) -> str:
        result = io.StringIO()
        for line in self.matrix:
            for point in line:
                if point is Material.AIR:
                    result.write(' ')
                elif point is Material.ROCK:
                    result.write('█')
                elif point is Material.SAND:
                    result.write('░')
                else:
                    assert False  # we covered all cases
            result.write('\n')
        return result.getvalue()

    def add_segment(self, start: Coords, end: Coords) -> None:
        y1, x1 = start
        y2, x2 = end
        if y1 == y2:
            x1, x2 = min(x1, x2), max(x1, x2)
            ys = (y1 for _ in range(x1, x2 + 1))
            xs = (x for x in range(x1, x2 + 1))
        elif x1 == x2:
            y1, y2 = min(y1, y2), max(y1, y2)
            ys = (y for y in range(y1, y2 + 1))
            xs = (x1 for _ in range(y1, y2 + 1))
        for y, x in zip(ys, xs):
            self.matrix[y, x] = Material.ROCK


def import_segments(
        lines: Iterable[str]
        ) -> Tuple[List[Coords], List[Tuple[Coords, Coords]]]:
    points: List[Coords] = []
    segments: List[Tuple[Coords, Coords]] = []

    def coords(word: str) -> Coords:
        parts = word.split(',')
        if len(parts) != 2:
            raise ValueError("unexpected number of coordinates")
        return int(parts[1]), int(parts[0])

    for line in lines:
        words = line.rstrip().split(' -> ')
        segment_points = [coords(word) for word in words]
        points.extend(segment_points)
        segments.extend(itertools.pairwise(segment_points))
    return points, segments


def import_slice1(lines: Iterable[str]) -> Slice:
    points, segments = import_segments(lines)
    y_leak, x_leak = 0, 500
    points.append((y_leak, x_leak))
    xmin = min(x for _, x in points) - 1
    xmax = max(x for _, x in points) + 1
    ymin = min(y for y, _ in points)
    ymax = max(y for y, _ in points)
    if ymin < 0:
        raise ValueError("unexpected negative coordinate")
    xshift = xmin
    matrix = np.full(((ymax + 1), xmax - xshift + 1), Material.AIR)
    result = Slice(matrix, (0, 500 - xshift))
    for (y1, x1), (y2, x2) in segments:
        result.add_segment((y1, x1 - xshift), (y2, x2 - xshift))
    return result


def import_slice2(lines: Iterable[str]) -> Slice:
    points, segments = import_segments(lines)
    y_leak, x_leak = 0, 500
    points.append((y_leak, x_leak))
    xmin = min(x for _, x in points)
    xmax = max(x for _, x in points)
    ymin = min(y for y, _ in points)
    ymax = max(y for y, _ in points) + 2  # including the floor
    # Update xmin and xmax to accomodate a pyramid of sand
    xmin = min(xmin, x_leak - (ymax - y_leak))
    xmax = max(xmax, x_leak + (ymax - y_leak))
    if ymin < 0:
        raise ValueError("unexpected negative coordinate")
    xshift = xmin
    matrix = np.full(((ymax + 1), xmax - xshift + 1), Material.AIR)
    result = Slice(matrix, (0, 500 - xshift))
    for (y1, x1), (y2, x2) in segments:
        result.add_segment((y1, x1 - xshift), (y2, x2 - xshift))
    result.add_segment((ymax, xmin - xshift), (ymax, xmax - xshift))
    return result


def solve1(lines: Iterable[str]) -> int:
    """Solve part 1 of today's puzzle"""
    slice_ = import_slice1(lines)
    for i in range(10000):
        if slice_.pour():
            return i
    raise ValueError("Simulantion did not end within expected limit")


def solve2(lines: Iterable[str]) -> int:
    """Solve part 2 of today's puzzle"""
    slice_ = import_slice2(lines)
    for i in range(100000):
        _ = slice_.pour()
        if slice_.matrix[slice_.leak] is Material.SAND:
            return i + 1
    raise ValueError("Simulantion did not end within expected limit")

On peut tricher en moins craignos. Un indice :

$ file 13.in
13.in: JSON data

Pour la définition de aoc.group_lines(lines: str), cf. hier

from __future__ import annotations

import itertools

from typing import Iterable, Iterator, List, Optional, Tuple, Union

import aoc


class Element:
    def __init__(self, value: Union[List[Element], int]):
        self.value = value

    def __lt__(self, other: Element) -> bool:
        if isinstance(self.value, int) and isinstance(other.value, int):
            return self.value < other.value
        if isinstance(self.value, list) and isinstance(other.value, list):
            for self_elt, other_elt in zip(self.value, other.value):
                if self_elt < other_elt:
                    return True
                if other_elt < self_elt:
                    return False
            # Ran out of elements on one of the lists
            return len(self.value) < len(other.value)
        if isinstance(self.value, int):
            return Element([self]) < other
        if isinstance(other.value, int):
            return self < Element([other])
        assert False  # should not happen: we covered all cases

    def __eq__(self, other: object) -> bool:
        if not isinstance(other, Element):
            return NotImplemented
        if isinstance(self.value, int) and isinstance(other.value, int):
            return self.value == other.value
        if isinstance(self.value, list) and isinstance(other.value, list):
            return (len(self.value) == len(other.value)
                    and all(self_elt == other_elt for self_elt, other_elt
                            in zip(self.value, other.value)))
        if isinstance(self.value, int):
            return Element([self]) == other
        if isinstance(other.value, int):
            return self == Element([other])
        assert False  # should not happen: we covered all cases

    def __str__(self) -> str:
        if isinstance(self.value, int):
            return str(self.value)
        elif isinstance(self.value, list):
            return '[{}]'.format(
                    ','.join(str(element) for element in self.value))
        assert False  # should not happen: we covered all cases


def import_element(chars: Iterable[str]) -> Element:
    element: Optional[Element] = None
    lists: List[List[Element]] = []  # List[List[Element]]
    current_list: Optional[List[Element]] = None
    current_int: Optional[int] = None
    for char in chars:
        if char == '[':
            if current_int is not None:
                raise ValueError("unexpected beginning of list")
            if current_list is not None:
                lists.append(current_list)
            current_list = []
        elif char == ']':
            if current_list is None:
                raise ValueError("unexpected end of list")
            # If we were parsing an int, add it before closing current list
            if current_int is not None:
                current_list.append(Element(current_int))
                current_int = None
            if lists:
                # We are closing a sub-list
                prev_list = lists.pop()
                prev_list.append(Element(current_list))
                current_list = prev_list
            else:
                # We are closing the top-level list
                element = Element(current_list)
                current_list = None
        elif char.isdecimal():
            value = int(char)
            if current_int is None:
                current_int = value
            else:
                current_int = 10 * current_int + value
                continue
        elif char == ',':
            if current_list is None:
                raise ValueError("unexpected separator")
            if current_int is not None:
                current_list.append(Element(current_int))
                current_int = None
        elif char == '\n':
            pass
        else:
            raise ValueError("unexpected character")
    if element is None:
        raise ValueError("nothing to parse")
    return element


def import_pairs(lines: Iterable[str]) -> Iterator[Tuple[Element, Element]]:
    for group in aoc.group_lines(lines):
        elements = [import_element(line) for line in group]
        if len(elements) != 2:
            raise ValueError("unexpected group length")
        yield (elements[0], elements[1])


def import_elements(lines: Iterable[str]) -> Iterator[Element]:
    for line in lines:
        if line and line != '\n':
            yield import_element(line)


def solve1(lines: Iterable[str]) -> int:
    """Solve part 1 of today's puzzle"""
    pairs = import_pairs(lines)
    return sum(i + 1 for i, (elt1, elt2) in enumerate(pairs) if elt1 < elt2)


def solve2(lines: Iterable[str]) -> int:
    """Solve part 2 of today's puzzle"""
    elements: Iterable[Element] = import_elements(lines)
    div1 = import_element("[[2]]")
    div2 = import_element("[[6]]")
    elements = sorted(itertools.chain(elements, (div1, div2)))
    key = 1
    for i, element in enumerate(elements):
        if element == div1 or element == div2:
            key *= i + 1
    return key

Au fait, vu les visualisations, vous avez évidemment remarqué que loin d'être une gorge de rivière, nous avions plutôt affaire à l'île de Numenor une montagne en forme d'étoile, n'est-ce pas ?

from __future__ import annotations

import itertools
import math
import re

from typing import Callable, Iterable, Iterator, List, Tuple


class Test:
    def __init__(self, div: int, rcpt1: int, rcpt2: int):
        self.div = div
        self.rcpt1 = rcpt1
        self.rcpt2 = rcpt2

    def __call__(self, n: int):
        if n % self.div == 0:
            return self.rcpt1
        else:
            return self.rcpt2


class Monkey:
    def __init__(self, number: int, items: List[int], op: Callable[[int], int], test: Test) -> None:
        self.number = number
        self.items = items
        self.op = op
        self.test = test
        self.activity = 0

    _re_monkey = re.compile(r'^Monkey (\d+):?\n?$')
    _re_items = re.compile(r'^\s*Starting items: (.*)\n?$')
    _re_op = re.compile(r'^\s*Operation: (.*)\n?$')
    _re_test = re.compile(r'^\s*Test: divisible by (\d+)\n?$')
    _re_true = re.compile(r'^\s*If true: throw to monkey (\d+)\n?$')
    _re_false = re.compile(r'^\s*If false: throw to monkey (\d+)\n?$')

    @classmethod
    def import_lines(class_, lines: Iterable[str]) -> Monkey:
        lines = iter(lines)
        number = int(class_._re_monkey.match(next(lines)).group(1))
        items = [int(word) for word in class_._re_items.match(next(lines)).group(1).split(', ')]
        op = import_op(class_._re_op.match(next(lines)).group(1))
        div = int(class_._re_test.match(next(lines)).group(1))
        rcpt1 = int(class_._re_true.match(next(lines)).group(1))
        rcpt2 = int(class_._re_false.match(next(lines)).group(1))
        test = Test(div, rcpt1, rcpt2)
        return class_(number, items, op, test)


def group_key() -> Callable[[str], int]:
    key = 0
    def aux(line: str):
        nonlocal key
        if line.rstrip() == '':
            key += 1
        return key
    return aux


def group_lines(lines: Iterable[str]) -> Iterable[Iterable[str]]:
    for _, group in itertools.groupby(lines, group_key()):
        yield (line for line in group if line.rstrip() != "")


_re_op_unary = re.compile("^new = old ([+*]) (\d+)$")
_re_op_binary = re.compile("^new = old ([+*]) old$")
def import_op(line: str) -> Callable[[int], int]:
    if (m := _re_op_unary.match(line)) is not None:
        arg = int(m.group(2))
        if m.group(1) == '+':
            return lambda old: old + arg
        if m.group(1) == '*':
            return lambda old: old * arg
        raise ValueError("unrecognized operator")
    if (m := _re_op_binary.match(line)) is not None:
        if m.group(1) == '+':
            return lambda old: 2 * old
        if m.group(1) == '*':
            return lambda old: old ** 2
        raise ValueError("unrecognized operator")
    raise ValueError("unrecognized operation arity")


class Game:
    def __init__(self, monkeys: dict[int, Monkey]) -> None:
        self.monkeys = monkeys

    def round(self) -> None:
        for monkey in self.monkeys.values():
            for item in monkey.items:
                item = monkey.op(item)
                item //= 3
                rcpt = monkey.test(item)
                monkey.activity += 1
                self.monkeys[rcpt].items.append(item)
            monkey.items = []

    @classmethod
    def import_lines(class_, lines: Iterable[str]) -> Game:
        monkeys = {}
        for group in group_lines(lines):
            monkey = Monkey.import_lines(group)
            monkeys[monkey.number] = monkey
        return class_(monkeys)

    def business(self) -> int:
        activity = sorted(monkey.activity for monkey in self.monkeys.values())
        return activity[-1] * activity[-2]


class Game2(Game):
    def __init__(self, *args, **kwargs) -> None:
        super().__init__(*args, **kwargs)
        self.div = math.lcm(*(monkey.test.div for monkey in self.monkeys.values()))

    def round(self) -> None:
        for monkey in self.monkeys.values():
            for item in monkey.items:
                item = monkey.op(item) % self.div
                rcpt = monkey.test(item)
                monkey.activity += 1
                self.monkeys[rcpt].items.append(item)
            monkey.items = []


def solve1(lines: Iterable[str]) -> int:
    """Solve part 1 of today's puzzle"""
    game = Game.import_lines(lines)
    for _ in range(20):
        game.round()
    return game.business()


def solve2(lines: Iterable[str]) -> int:
    """Solve part 2 of today's puzzle"""
    game = Game2.import_lines(lines)
    for _ in range(10000):
        game.round()
    return game.business()

Ma première idée fonctionnait sur l'exemple mais pas sur les données réelles : parcourir récursivement une matrice de proche en proche en évitant simplement les endroits où on est déjà passé, ça atteint vite la limites de longueur de la pile d'appels.

Du coup, j'ai fait du Dijkstra. Enfin quelque chose inspiré de son algorithme en tout cas. Ça pourrait être fait de façon entièrement itérative, mais ça ne valait pas la peine, ça reste raisonnable en récursif.

Une fois qu'on a ça, la deuxième partie du puzzle ne pose pas spécialement de problème. Il y a tout de même deux façons de l'implémenter, une bête et méchante (on prend la première solution et on l'applique plusieurs fois), et une un rien plus futée.

from typing import Iterable, Iterator, List, Optional, Set, Tuple

import numpy as np


Coords = Tuple[int, int]


class Map:
    def __init__(self, matrix: np.ndarray) -> None:
        self.matrix = matrix
        self.ly, self.lx = matrix.shape

    def neighs(self, y: int, x: int) -> Iterator[Coords]:
        """Yield neighbours that are reachable from the given coordinates,
        considering movement rules (no climbing)"""
        for (dx, dy) in ((-1, 0), (1, 0), (0, -1), (0, 1)):
            y_ = y + dy
            x_ = x + dx
            if y_ < 0 or y_ >= self.ly or x_ < 0 or x_ >= self.lx:
                continue
            if self.matrix[y_, x_] - self.matrix[y, x] <= 1:
                yield y_, x_

    def _distances(self, starts: Iterable[Coords], end: Coords,
                   distances: np.ndarray) -> None:
        """Update distances matrix by computing walking distance from possible
        starts"""
        if starts == []:
            # Nothing left to explore
            return
        nexts = []  # List[Coords]
        for start in starts:
            start_dist = distances[start]
            if start_dist < 0:
                raise ValueError(
                        'cannot compute distances from uncharted point')
            for neigh in self.neighs(*start):
                neigh_dist = distances[neigh]
                if neigh_dist < 0 or start_dist + 1 < neigh_dist:
                    # This point has either never been checked before, or has
                    # been but we have a shorter path to it
                    distances[neigh] = start_dist + 1
                    nexts.append(neigh)
        # Update distances from the points we have just updated
        self._distances(nexts, end, distances)

    def min_dist(self, starts: List[Coords], end: Coords) -> int:
        distances = np.full_like(self.matrix, -1)
        """Return the minimal distance to reach end, starting from starts"""
        for start in starts:
            distances[start] = 0
        self._distances(starts, end, distances)
        return distances[end]


def import_lines(lines: Iterable[str]) -> Tuple[Map, Coords, Coords]:
    matrix = []  # type: List[List[int]]
    start = None
    end = None
    for y, line in enumerate(lines):
        matrix.append([])
        for x, char in enumerate(line.rstrip()):
            if char == 'S':
                start = (y, x)
                height = 0
            elif char == 'E':
                end = (y, x)
                height = 25
            else:
                height = ord(char) - ord('a')
            matrix[-1].append(height)
    if start is None or end is None:
        raise ValueError("no start or end position found")
    return Map(np.array(matrix)), start, end


def solve_both(lines: Iterable[str]) -> Tuple[int, int]:
    """Solve part 1 of today's puzzle"""
    map_, start, end = import_lines(lines)
    min1 = map_.min_dist([start], end)
    min2 = map_.min_dist([coords for coords, height  # type: ignore
                          in np.ndenumerate(map_.matrix)
                          if height == 0], end)
    return min1, min2