Défi de Clem // Percolation // un exemple de code dans un langage 'atypique'.

PROCEDURE MaFenêtre()
nlig, ncol est un entier
la0, ha0 est un entier      // utile uniquement pour la représentation dans le champ image. 
n_arbres_initial est un entier
n_arbres_brules est un entier

pct_clairiere est un réel
pct_p est un réel

tb est un tableau dynamique // Etat de chaque case ( ou encore de chaque arbre)

q_etat est une Enumération
    clairiere
    Arbre_sain
    arbre_en_feu
    arbre_en_cendres
FIN
liste_arbres_en_feu est un tableau de 0 par 2 entiers

i,j est un entier
r est un réel
i1,j1 , i2,j2  est un entier

nlig = 200 
ncol = 200
la0 = PartieEntière(Image1..Largeur  / ncol)
ha0 = PartieEntière(Image1..Hauteur  / nlig)
Image1..Largeur = la0 * ncol
Image1..Hauteur = ha0 * nlig

pct_p = 0.59
i = Saisie( " Donnez-le pourcentage de propagation du feu d'un arbe à l'autre",pct_p )
SI i = -1 ALORS RENVOYER Faux  // L'utilisateur a cliqué sur Annuler
SI pct_p < 0 _OU_ pct_p > 1 ALORS 
    Erreur ( " Le nombre saisi est uinvalide "," Abandon ") 
    RENVOYER Faux
FIN

pct_clairiere = 0.02
i = Saisie( " Donnez-le pourcentage de clairiere dans la grille ",pct_clairiere )
SI i = -1 ALORS RENVOYER Faux  // L'utilisateur a cliqué sur Annuler
SI pct_clairiere < 0 _OU_ pct_clairiere > 1 ALORS 
    Erreur ( " Le nombre saisi est uinvalide "," Abandon ") 
    RENVOYER Faux
FIN

tb = allouer un tableau de nlig par ncol q_etat
// Initialisation clairiere // Arbre
InitHasard( HeureVersEntier(HeureSys()))
dDébutDessin(Image1)

n_arbres_initial = 0
POUR i = 1 _A_ nlig
    POUR j = 1 _A_ ncol
        r = Hasard()
        SI r > pct_clairiere ALORS 
            tb[i,j]= Arbre_sain
            n_arbres_initial ++
        SINON
            tb[i,j] = clairiere
        FIN
        f_dessine(i,j )
    FIN
FIN
SI n_arbres_initial = 0 ALORS 
    Erreur ( " Aucun arbre, uniquement de la clairiere ! ", "fin du traitement")
    RENVOYER Faux
FIN

// Choisir un arbre  où le feu se déclenche
i1 = Hasard(nlig)
j1 = Hasard(ncol)
TANTQUE tb[i1,j1] <> Arbre_sain  // Boucle 
    i1 = Hasard(nlig)
    j1 = Hasard(ncol)
FIN
tb[i1,j1]= arbre_en_feu
TableauSupprimeTout(liste_arbres_en_feu)
TableauAjouteLigne(liste_arbres_en_feu,i1,j1)
f_dessine(i1,j1)
n_arbres_brules =1


// Et debut de la propagation
TANTQUE TableauOccurrence( liste_arbres_en_feu) > 0
    // Je prends le premier arbre, je prends ses 4 voisins, et je propage le feu, avec une probablilité p
    i1 = liste_arbres_en_feu[1,1]
    j1 = liste_arbres_en_feu[1,2]
    TableauSupprimeLigne(liste_arbres_en_feu,1)
    POUR direction  = 1 A 4
        SELON direction
            CAS 1 
                i2=i1+1  ;  j2= j1
            CAS 2 
                i2=i1-1  ;  j2= j1
            CAS 3 
                i2=i1  ;  j2= j1 -1
            CAS 4 
                i2=i1  ;  j2= j1+1
            AUTRE CAS   
        FIN
        // Ici, on pourrait aménager pour traiter le cas du tore.
        SI i2 > 0 _ET_ i2<= nlig _ET_ j2 > 0 _ET_ j2 <= ncol ALORS
            SI tb[i2,j2] = Arbre_sain ALORS 
                r = Hasard()
                SI r <= pct_p ALORS 
                    tb[i2,j2] = arbre_en_feu
                    f_dessine(i2,j2)
                    TableauAjouteLigne(liste_arbres_en_feu,i2,j2) 
                    n_arbres_brules++
                FIN
            FIN
        FIN
    FIN
    tb[i1,j1]= arbre_en_cendres
    f_dessine(i1,j1)
FIN
Info( n_arbres_brules, n_arbres_initial)

PROCEDURE f_dessine(i,j)
coul est un entier
qq est un q_etat
qq= tb[i,j]
SELON qq
    CAS clairiere
        coul = JauneClair
    CAS Arbre_sain
        coul = VertClair
    CAS arbre_en_feu
        coul = RougeClair
    CAS arbre_en_cendres
        coul = GrisFoncé    
FIN
dRectangle(  (j-1)* la0, ( i-1)*ha0, j*la0,i*ha0, coul,coul)
MultitâcheRepeint()
RENVOYER Vrai