Fractales

Sur cette page de Wolfram, on peut voir différentes fractales :

Et en bas de cette même page les algorithmes utilisés pour transformer les cases :

Ci-dessous un programme Python qui permet d’appliquer le motif voulu, il suffit d’indiquer les coordonnées des cases blanches. Par exemple pour le (a), la seule case blanche est au milieu (coordonnées (1,1))

La case blanche en (1,1)
from PIL import Image, ImageDraw
import numpy as np

# Dimensions de l'image finale (à modifier comme vous le voulez)
source = Image.new("RGB", (600, 600), color="white")
draw = ImageDraw.Draw(source)

# Recherche s'il y a une case blanche
def has_common(str1, str2, pattern):
    min_len = min(len(str1), len(str2))
    for i in range(1, min_len + 1):
        for p in pattern:
            # Ne pas afficher si on trouve une case blanche
            if str1[-i] == p[0] and str2[-i] == p[1]:
                return False
    # Sinon case noire
    return True

# Construction de la fractale
def fractal(pattern):
    # On parcourt les colonnes
    for c in range(600):
        # Convertir 'c' en base 3 (ici à l'aide de numpy)
        t1 = np.base_repr(c, base=3)
        # On parcourt les lignes
        for l in range(600):
            t2 = np.base_repr(l, base=3)
            # Si pas de case blanche, afficher le point
            if has_common(t1, t2, pattern):
                draw.point((c, l), fill=(0, 0, 0))

# Exemple avec le tapis de Sierpiński
fractal(["11"])
source.show()
Motif (a) = case blanche en (1,1)

Quelques motifs

(b) = fractal([’00’,’02’,’20’,’22’])
(d) = fractal([’10’,’20’,’21’])
(e) = fractal([’00’,’02’,’11’,’20’,’22’])
(g) = fractal([’11’,’20’])
(h) = fractal([’02’,’11’,’20’])

Et on peut en inventer de nouveaux :

fractal([’02’,’10’,’20’])

Amusons-nous avec des mots du dictionnaire

Exercice 1

Écrire un programme qui cherche tous les mots de 5 lettres pouvant être formés à partir de mots de 9 lettres en prenant 1 lettre sur 2, en commençant par la première lettre. Exemples :

Une lettre sur deux

Vous pouvez utiliser ces 2 dictionnaires (mots de 9 et 5 lettres) : https://uabox.univ-angers.fr/s/b3GotFzpTYkL4MT

Corrigé en python

with open('dictionary_9.txt', 'r') as f:
    dictionary_9 = [line.strip() for line in f]

with open('dictionary_5.txt', 'r') as f:
    dictionary_5 = set(line.strip() for line in f)

compte = 0                      # Nb de mots trouvés

for word9 in dictionary_9:
  word5 = word9[::2]            # une lettre sur deux
  if word5 in dictionary_5:     # Si le mot existe
    print(word9, word5)         # on l'affiche
    compte += 1                 # et le compteur augmente de +1
print(f"Total = {compte}")  

Résultat :

ABLATIONS ALTOS
ABLUTIONS ALTOS
ABREGEAIS ARGAS
ABROGEAIS ARGAS
ABSTIENNE ASINE
ACCONIERS ACNES
...
VIELLEUSE VELUE
VIELLIONS VELOS
VIENDRAIS VEDAS
VOILETTES VIETS
Total = 913

corrigé en javascript

Cliquez sur ce lien puis bouton droit – Inspecter – Console. Copiez-collez le code suivant :

 dictionary_9.reduce((a, m) => {
     var mot5 = m.slice(0,2)+m[4]+m.slice(-2);
     return dictionary_5.includes(mot5) ? [...a, [m, mot5]] : a 
  }, [])

Exercice 2

Identifiez tous les mots de 10 lettres qui peuvent être composés de 2 mots de 5 lettres en respectant l’ordre des lettres. Exemple : RACHIDIENS s’écrit à partir de CHIEN et RADIS

Comme il y a beaucoup de solutions, on peut créer une fonction qui admet en paramètre un mot de 5 lettres et qui renvoie tous les mots de 10 lettres contenant ce mot ainsi que l’autre mot de 5 lettres pour compléter. Exemples :

>>> trouve("CHIEN")
BANCHAIENT CHIEN BANAT
CHARMAIENT CHIEN ARMAT
CHERRAIENT CHIEN ERRAT
CHIADERENT CHIEN ADRET
CHIENDENTS CHIEN DENTS
CHTHONIENS CHIEN THONS
CHTONIENNE CHIEN TONNE
CRASHAIENT CHIEN RASAT
DOUCHAIENT CHIEN DOUAT
LOUCHAIENT CHIEN LOUAT
MATCHAIENT CHIEN MATAT
RACHIDIENS CHIEN RADIS
TOUCHAIENT CHIEN TOUAT
TRICHAIENT CHIEN TRIAT

Corrigé en python

with open('dictionary_10.txt', 'r') as f:
    dictionary_10 = [line.strip() for line in f]

with open('dictionary_5.txt', 'r') as f:
    dictionary_5 = set(line.strip() for line in f)

def is_included(word1, word2):
    positions = []
    index = 0
    for letter in word2:
        index = word1.find(letter, index)
        if index == -1: return False
        positions.append(index)
        index += 1      
    return positions

def trouve(word5):
 for word10 in dictionary_10:
  positions = is_included(word10, word5)  
  if positions:
    reste = ''.join(word10[i] for i in range(10) if i not in positions)
    if reste in dictionary_5:
       print(word10, word5, reste)

exercice 3

Ci-dessous un jeu très simple trouvé dans le journal TV Télé Z. Dans notre cas nous allons travailler avec des mots de 9 et 5 lettres.

But : créer un programme qui, à partir d’une liste de mots de 5 lettres, va rechercher autant de mots de 9 lettres les contenant.

>>> jeu(["AIDER","CAPTER","GAINE","LOGER","FILES"])
PYRAMIDER AIDER PYR-M----
CAPTERAIT CAPTER ------AIT
GRAINIERS GAINE -R---I-RS
PLONGEOIR LOGER P--N--OI-
FILONIENS FILES ---ONI-N-

corrigé en python

from random import choice

with open('dictionary_9.txt', 'r') as f:
    dictionary_9 = [line.strip() for line in f]

with open('dictionary_5.txt', 'r') as f:
    dictionary_5 = set(line.strip() for line in f)

def is_included(word1, word2):
    positions = []
    index = 0
    for letter in word2:
        index = word1.find(letter, index)
        if index == -1: return False
        positions.append(index)
        index += 1      
    return positions

def trouve(word5):
 res = []         # On cherche toutes les solutions
 for word9 in dictionary_9:
  positions = is_included(word9, word5)  
  if positions: res.append((word9, positions))
 if len(res) > 0: return choice(res)      # on renvoie une solution au hasard
 return False

def jeu(arr):
 for mot in arr:
  r = trouve(mot)
  if r:     # Si un mot de 9 lettres a été trouvé
   cache = ''.join(r[0][i] if i not in r[1] else '-' for i in range(9))
   print(r[0], mot, cache)
  else: print(f"Rien trouvé pour {mot}")    

Quelques programmes pour les CASIO fx-180P

Présentation générale

Quelques modèles de ma collection

fx-180P, 180Pv, 180PA et 180P Plus

Autre version du programme Syracuse

Pas	Touches
01	2
02	Kin ÷ 1		# n / 2
03	Kout 1
04	-
05	1
06	Kin + 3 	# TOF = TOF + 1
07	=
08	1/x		# Si n = 1 arrêt sur une erreur
09	Kout 1
10	-
11	RND		# Partie entière quand on est en FIX 0
12	+
13	.
14	5
15	=		# Si n est pair on obtient 1 et 0 sinon
16	x > 0		# Retour pas n°1 si n pair
17	6		# Sinon...
18	Kin × 1
19	1
20	Kin + 1		# 3n+1
21	Kout 1
22	x ≤ M		# Si n <= max retour au pas n°1
23	Min		# sinon mise à jour du max
24	x > 0		

Exemple d’utilisation pour N = 27 puis N = 17 :

MODE FIX 0
KAC
27 Min Kin 1 P1

Le programme s'arrête sur -E-
AC puis MR pour voir le maximum (9232) et Kout 3 pour le temps de vol (111)

Pour un autre calcul faire : 
KAC
17 Min Kin 1 P1
Résultats : max = 52 et TOF = 12

Décomposition d’un entier en facteurs premiers

Programme P1 :

001	Min
002	1
003	Kin 1
004	Mode 7 0    # Fix 0

Programme P2 :

001	MR
002	÷
003	1
004	Kin + 1
005	Kout 1
006	-
007	RND
008	=
009	x > 0
010	+/-
011	x > 0       # Cela revient à tester x ≠ 0
012	Kout 1
013	HLT
014	1/x
015	×
016	MR
017	=
018	Min
019	1
020	Kin - 1
021	Kout 1
022	x ≤ M
023	AC

Utilisation : Décomposition de 2024 en facteurs premiers

2024 P1
P2
Affichage : 2
RUN
Affichage : 2
RUN
Affichage : 2
RUN
Affichage : 11
RUN
Affichage : 23
RUN
0

Donc 2024 = 2 × 2 × 2 × 11 × 23

Lorsque le nombre est premier :

13 P1
P2
Affichage : 13
RUN
0

Représentation du temps – Informatique et temps chronologique

Conflit de ressources

import threading
import time

# Nombre de tickets disponibles (ressource partagée)
nb_tickets  = 10

def achat_ticket():
  global nb_tickets 
  if nb_tickets  > 0:
    # Simuler une lecture de la base de données
    current_tickets = nb_tickets 
    # Simuler un délai
    time.sleep(0.1)
    # Simuler une écriture dans la base de données
    nb_tickets  = current_tickets - 1
    print(f"Ticket acheté, tickets restants : {nb_tickets }\n")
  else: print("Plus de tickets disponibles.\n")

# Plusieurs utilisateurs essayant d'acheter des tickets en même temps

for _ in range(20):
  threading.Thread(target=achat_ticket).start()

Version avec verrou

import threading
import time

# Nombre de tickets disponibles (ressource partagée)
nb_tickets  = 10
lock = threading.Lock()

def achat_ticket():
  global nb_tickets
  with lock:
   if nb_tickets  > 0:
    # Simuler une lecture de la base de données
    current_tickets = nb_tickets 
    # Simuler un délai
    time.sleep(0.1)
    # Simuler une écriture dans la base de données
    nb_tickets  = current_tickets - 1
    print(f"Ticket acheté, tickets restants : {nb_tickets }\n")
   else: print("Plus de tickets disponibles.\n")

# Plusieurs utilisateurs essayant d'acheter des tickets en même temps

for _ in range(20):
  threading.Thread(target=achat_ticket).start()

Version SQL

import sqlite3

def buy_ticket(event_name):
    try:
        # Connexion à la base de données SQLite (ou création de la base si elle n'existe pas)
        conn = sqlite3.connect('example.db')
        cursor = conn.cursor()

        # Démarrer une transaction
        cursor.execute('BEGIN TRANSACTION')

        # Étape 1 : Vérifier le nombre de tickets disponibles
        cursor.execute('SELECT available_tickets FROM tickets WHERE event_name = ?', (event_name,))
        result = cursor.fetchone()

        if result is None:
            print(f"Événement '{event_name}' introuvable.")
            conn.rollback()
            return

        available_tickets = result[0]
        print(f"Tickets disponibles avant l'achat : {available_tickets}")

        if available_tickets > 0:
            # Étape 2 : Réserver un ticket
            cursor.execute('UPDATE tickets SET available_tickets = available_tickets - 1 WHERE event_name = ?', (event_name,))
            print("Achat réussi !")

            # Valider la transaction
            conn.commit()
        else:
            print("Échec de l'achat : plus de tickets disponibles.")
            conn.rollback()

    except sqlite3.Error as e:
        print(f"Erreur SQL : {e}")
        conn.rollback()
    finally:
        conn.close()

# Exemple d'utilisation avec plusieurs "acheteurs"
if __name__ == "__main__":
    # Création de la table et insertion initiale (exécuté une seule fois)
    conn = sqlite3.connect('example.db')
    cursor = conn.cursor()
    cursor.execute('''
    CREATE TABLE IF NOT EXISTS tickets (
        id INTEGER PRIMARY KEY AUTOINCREMENT,
        event_name TEXT NOT NULL,
        available_tickets INTEGER NOT NULL
    )
    ''')
    cursor.execute('INSERT INTO tickets (event_name, available_tickets) VALUES (?, ?)', ('concert', 100))
    conn.commit()
    conn.close()

    # Simulation d'achats concurrents
    for _ in range(55): buy_ticket('concert')

Fuseaux horaires

from datetime import datetime
import pytz

# Obtenir l'heure actuelle en UTC
utc_now = datetime.now(pytz.utc)
print(f"Heure actuelle en UTC : {utc_now}")

for endroit in 'Europe/Paris', 'America/New_York', 'Africa/Cairo', 'Australia/Sydney', 'Asia/Tokyo', 'Antarctica/Troll':
 local_tz = pytz.timezone(endroit)
 local_time = utc_now.astimezone(local_tz)
 print(f"Heure locale {endroit} : {local_time}")

Représentation du temps – Dans les systèmes informatiques

Remarque sur les calculatrices TI : comme précisé dans la vidéo, la fonction dbd des TI-82/83/84 permet de trouver le nombre de jours entre 2 dates entre 1950 et 2049. Cependant, sur les TI-84 Plus CE cette plage a été décalée à l’intervalle 1980 – 2079.

Temps unix

import datetime
from datetime import timezone

dt = datetime.datetime(2024,8,5,14,33,17, tzinfo=timezone.utc)
>>> dt.timestamp()
1722868397.0

import time
>>> time.gmtime(1722868397.0)
time.struct_time(tm_year=2024, tm_mon=8, tm_mday=5, tm_hour=14, tm_min=33, tm_sec=17, tm_wday=0, tm_yday=218, tm_isdst=0)

>>> time.gmtime(0)
time.struct_time(tm_year=1970, tm_mon=1, tm_mday=1, tm_hour=0, tm_min=0, tm_sec=0, tm_wday=3, tm_yday=1, tm_isdst=0)

>>> time.gmtime(-2000 * 365 * 86400)
time.struct_time(tm_year=-29, tm_mon=5, tm_mday=1, tm_hour=0, tm_min=0, tm_sec=0, tm_wday=5, tm_yday=121, tm_isdst=0)

>>> time.time()
1722867076.132498

>>> datetime.datetime.now(timezone.utc).timestamp()

Date de pâques (Formule approximative !)

def paques(annee):
    cycle_lune = 29.53058853
    j_sem = 6  # Samedi 5/1/1901
    delta = nb_jours(5, 1, 1901, 21, 3, annee)
    nb_cycles = 1 + int(delta / cycle_lune)
    ajout = int(nb_cycles * cycle_lune)
    while (j_sem + ajout) % 7 != 0: ajout +=1
    nb = ajout - delta
    if nb <= 10: j, m = 21 + nb, 3
    else: j, m = nb - 10, 4
    return (m, j)

def formule(j,m,a):
    if m < 3:
        a -= 1
        m += 12
    return int(365.25 * a) + int(30.6 * (m + 1)) + j - 3

def nb_jours(j1, m1, a1, j2, m2, a2):
    return formule(j2, m2, a2) - formule(j1, m1, a1)

def dow(j,m,a): return 1 + (formule(j,m,a) % 7)

Date de pâques (formule exacte)

def easter(annee):
    a = annee % 19
    b = annee // 100
    c = annee % 100
    d = b // 4
    e = b % 4
    f = (b + 8) // 25
    g = (b - f + 1) // 3
    h = (19 * a + b - d - g + 15) % 30
    i = c // 4
    k = c % 4
    l = (32 + 2 * e + 2 * i - h - k) % 7
    m = (a + 11 * h + 22 * l) // 451
    mois = (h + l - 7 * m + 114) // 31
    jour = ((h + l - 7 * m + 114) % 31) + 1
    return (annee, mois, jour)

Représentation du temps – Calculs sur les dates et heures

Pourquoi 60 ?

import math

def nb_diviseurs(n):
    count = 0
    for i in range(1, int(math.sqrt(n)) + 1):
        if n % i == 0:
            if n == i * i: count += 1
            else: count += 2
    return count

maxi = 0
for n in range(1, 100):
    t = nb_diviseurs(n)
    if t > maxi:
        print(f"{n} a {t} diviseurs")
        maxi = t

# Résultats :

1 a 1 diviseur
2 a 2 diviseurs
4 a 3 diviseurs
6 a 4 diviseurs
12 a 6 diviseurs
24 a 8 diviseurs
36 a 9 diviseurs
48 a 10 diviseurs
60 a 12 diviseurs

Nombre de jours entre 2 dates

mois_jours = (31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31)

def est_bissextile(annee):
    return (annee % 4 == 0 and annee % 100 != 0) or annee % 400 == 0

def jours_dans_annee(jour, mois, annee):
    jours = jour
    for m in range(1, mois):
        jours += mois_jours[m - 1]
    if mois > 2 and est_bissextile(annee): jours += 1
    return jours

def jours_entre_dates(j1, m1, a1, j2, m2, a2):
    if (a1, m1, j1) > (a2, m2, j2):
        j1, m1, a1, j2, m2, a2 = j2, m2, a2, j1, m1, a1
    jours_total = 0
    if a1 == a2:
        return jours_dans_annee(j2, m2, a2) - jours_dans_annee(j1, m1, a1)  
    jours_total += (366 if est_bissextile(a1) else 365) - jours_dans_annee(j1, m1, a1)
    jours_total += jours_dans_annee(j2, m2, a2)  
    for annee in range(a1 + 1, a2):
        jours_total += 366 if est_bissextile(annee) else 365  
    return jours_total

# Exemple

>>> jours_entre_dates(15,10,1582,4,8,2024)
161366

Programme CULMINATION pour les ti-82/83/84

Disp "JJMM.AA ?"
Input D
Disp "LONGITUDE ?"
Input L
dbd(101+fPart(D),D→N
2π(N-81)/365→B
7.678sin(B+1.374)-9.87sin(2B)→E
Disp "CULMINATION"
12+E/60-L/15▸DMS

Programme LONGITUDE pour les ti-82/83/84

Disp "CULMINATION"
Disp "HH°MM'SS''?"
Input H
getDate→L₁
100L₁(3)+L₁(2)+remainder(L₁(1),100)/100→D
dbd(101+fPart(D),D→N
2π(N-81)/365→B
7.678sin(B+1.374)-9.87sin(2B)→E
Disp "LONGITUDE :"
15*(12+E/60-H)
Longitude trouvée pour le 2 août 2024

Représentation du temps : Calendriers

Phases de la lune

cycle_lunaire = 29.53058853
delta = cycle_lunaire / 8
phases = "Nouvelle lune","Premier croissant","Premier quartier", "Gibbeuse croissante","Pleine lune","Gibbeuse décroissante","Dernier quartier","Dernier croissant"
visuel = ['🌑','🌒','🌓','🌔','🌕','🌖','🌗','🌘']

def phase_lune(nb_jours):
 nb_cycles = int(nb_jours / cycle_lunaire)
 reste = nb_jours - nb_cycles * cycle_lunaire
 tranche = int(round(reste / delta, 0) % 8)
 return phases[tranche], visuel[tranche]

>>> phase_lune(365)
('Gibbeuse croissante', '🌔')

Calendrier chinois – signes du zodiaque

def zodiaque_chinois(annee):
 elements = "Bois", "Feu", "Terre", "Métal", "Eau"
 aspects = "Yang", "Yin"
 signes_zodiaque = (
        "Rat", "Buffle", "Tigre", "Lapin", "Dragon", "Serpent",
        "Cheval", "Chèvre", "Singe", "Coq", "Chien", "Cochon"
    )
 annee_ref = 1924 # Année du Rat, début d'un cycle de 60 ans
 diff_annees = annee - annee_ref
 ind_signe_zodiaque = diff_annees % 12
 ind_aspect = diff_annees % 2
 ind_element = (diff_annees % 10) // 2
 signe_zodiaque = signes_zodiaque[ind_signe_zodiaque]
 aspect = aspects[ind_aspect]
 element = elements[ind_element]
 annee_chinoise = annee_ref + (diff_annees % 60)
 return f"{annee_chinoise} : {signe_zodiaque} ({element} {aspect})"
 
>>> zodiaque_chinois(2024)
1964 : Dragon (Bois Yang)

Comparaison JULIEn | GREGORIEN | Personnel

def bissextile_julien(annee):
 return annee % 4 == 0
 
def bissextile_gregorien(annee):
 return (annee % 4 == 0 and annee % 100 != 0) or (annee % 400 == 0)
 
def bissextile_perso(annee):
 return annee % 4 == 0 and annee % 128 != 0
 
def compare(ans, methode):
 total = 0
 for an in range(ans):
  total += 365
  if methode(an): total += 1
 return total - 365.2421904482 * ans

>>> compare(10000, bissextile_julien)
78.09551799995825

>>> compare(10000, bissextile_gregorien)
3.095517999958247

>>> compare(10000, bissextile_perso)
-0.904482000041753

One-liners sur HP Prime

Réponses possibles aux 2 questions finales

ZIPSORT

cat(EXECON("CHAR(SORT({&1,&2}))",ASC(A),ASC(B)))

Exemple de déroulement avec zipsort(“oui”,”non”) :

  • On transforme les 2 chaines A et B en listes de codes ASCII : {111,117,105} et {110,111,110}
  • On trie les couples de lettres SORT({&1,&2}) : SORT({111,110}) -> {110,111}
  • Transformation en chaine : CHAR({110,111}) -> “no”
  • On obtient la liste : {“no”,”ou”,”in”}
  • Concaténation finale : “noouin”

Lettres manquantes

cat(EXECON("CAS((&2≠1+&1)*char(1+&1))",ASC(T)))

Exemple de déroulement avec manque(“abcdfg”) :

  • On transforme la chaine en liste ASCII : {97,98,99,100,102,103}
  • Est-ce que l’élément suivant est différent de l’élément courant + 1 ?
  • EXECON(“&2≠1+&1”,{97,98,99,100,102,103}) -> {0,0,0,1,0}
  • La lettre manquante a comme code ASCII la lettre courante + 1 : char(1+&1)
  • En multipliant le caractère par 0 on a une chaine vide, sinon on récupère le caractère manquant
  • Concaténation des chaines vides + lettres manquantes
  • manque(“abdefgijklmnpqrstv”) -> “chou”

Quelques challenges (proposition de solutions plus bas)

Autant de “x” que de “o”

Afficher 1 s’il y a autant de “x” que de “o” dans une chaine de caractère, sinon afficher 0.

xox("xxabcoo") -> 1
xox("xxooxoox") -> 1
xox("xaxoboxocccoox") -> 0

Au milieu

Ecrire une fonction qui prend en paramètres un caractère X et le place au milieu de Y répété N fois. Lorsque ce n’est pas possible, renvoyer X.

Exemples :

middle(10,"A","*") --> "*****A*****"  ("A" est au milieu de 10 "*")
middle(9,"A","*") --> "A"             ("A" ne peut pas être au milieu)
middle(2,"X","+") --> "+X+"

mot pur

Un mot pur est un mot dont la somme des positions dans l’alphabet de chaque lettre est divisible par la longueur totale du mot.

Par exemple, “abcb” est un mot pur car 1 + 2 + 3 + 2 = 8 et 8/4 = 2.

pure("ccc") -> 1
pure("bed") -> 0

Palindrome

Renvoyer 1 si une chaine est un palindrome, 0 sinon. La chaine pourra être écrite en minuscules et/ou majuscules.

palind("Laval") -> 1
palind("Angers") -> 0

Somme des N plus grands

Ecrire une fonction qui à partir d’une liste L et d’un entier N renvoie la somme des N entiers les plus grands de L.

sumgrand({4,9,2,3,7,1},2) --> 16    (Les 2 nombres les + grands sont 9 et 7)
sumgrand({-7,9,12,-1,-3},3) --> 20  (9 + 12 - 1 = 20)

Somme sans doublons

Ecrire une fonction qui fait la somme des éléments d’une liste, mais ignore ceux qui sont dupliqués.

Exemples : pour la liste [3, 4, 3, 6] la fonction devra renvoyer 10
et pour la liste [1, 10, 3, 10, 10] la fonction devra renvoyer 4.

Project Euler n°1

L’énoncé est ici

Exemples de solutions

Autant de “x” que de “o”

count_eq(120,ASC(T))=count_eq(111,ASC(T))

Au milieu

IFTE(N MOD 2,X,CAS(N/2*Y+X+N/2*Y))

mot pur

0=ΣLIST(ASC(T) MOD SIZE(T))

Palindrome

UPPER(T)==CHAR(revlist(ASC(UPPER(T))))

Remarquez que cette version fonctionne également (en décochant "L1") : 

revlist(sto(ASC(UPPER(T)),L1))==L1

Somme des N plus grands

ΣLIST(CAS.mid(CAS.SORT(L,"(x,y)->x>y"),1,N))

Somme sans doublons

ΣLIST(apply("x->x*(count_eq(x,L)==1)",L))

Project euler n°1

Décochez "v" :

ΣLIST(remove("x->(x MOD 3)*(x MOD 5)>0",MAKELIST(v,v,1,N)))-N

PROLOG

Tous les codes de la vidéo

Les collectionneurs

possede(cret,calc(hp71b,1984)).
possede(eric,ordi(zx81,1981)).
possede(ledudu,ordi(c64,1982)).
possede(eric,calc(hp41C,1979)).
possede(nicole,calc(pb100,1983)).
possede(audrey,ordi(apple2,1977)).
possede(cret,ordi(c64,1982)).
possede(ledudu,calc(ti57,1977)).
collection(P):-
 possede(P,ordi(_,_)),
 possede(P,calc(_,_)).

annee(B,H,X):-
    B > H,
    annee(H,B,X).

annee(B,H,X):-
    possede(_,calc(X,A)),
    A =< H,
    A >= B.

annee(B,H,X):-
    possede(_,ordi(X,A)),
    A =< H,
    A >= B. 

?- possede(P,ordi(_,1977)).
?- possede(P,T).
?- possede(_,calc(N,_)).
?- possede(X,calc(_,T)),T<1980.

Nombre d’or

Version 1 :

nombreOr(N, R) :- fc(N, 1, R).
fc(1, A, R) :- R = A.
fc(N, A, R) :-
    S is 1 + 1 / A,
    N1 is N - 1,
    fc(N1, S, R). 

Version 2 :

nombreOr(N, Num, Den) :- fc(N, 1, 1, Num, Den).
 
fc(1, Num, Den, Num, Den).
 
fc(N, NumP, DenP, Num, Den) :-
    N1 is N - 1,
    Num2 is NumP + DenP,
    fc(N1, Num2, NumP, Num, Den).

Taille d’une liste

length([],0).

length([H|T],N):-
    length(T,N1), N is N1+1. 

Triplets pythagoriciens

Version 1 :

valide(A, B, C) :-
    C2 is A*A + B*B,
    C*C =:= C2.

cherche(A, B, C, Max) :-
    between(1, Max, A),
    between(A, Max, B),
    between(B, Max, C),
    valide(A, B, C).

pythagore(Max, R) :-
    findall((A,B,C), cherche(A, B, C, Max), R). 

Version 2 :

entre(B, H, B).
entre(B, H, X) :-
    B < H,
    B1 is B + 1,
    entre(B1, H, X).

valide(A, B, C) :-
    C2 is A*A + B*B,
    C*C =:= C2.

cherche(A, B, C, Max) :-
    entre(1, Max, A),
    entre(A, Max, B),
    entre(B, Max, C),
    valide(A, B, C).

pythagore(Max) :-
    cherche(A, B, C, Max),
    L = [A,B,C],
    write(L), nl,
    fail.

Syracuse

Version 1 :

syr(N):- syr(N,0,N).
syr(1,T,M):- write(T),nl,write(M).
syr(N,T,M):-
 T2 is T + 1,
 suiv(N,T2,M).

suiv(N,T,M):-
 0 is N mod 2,
 N2 is N/2,
 syr(N2,T,M).

suiv(N,T,M):-
 N2 is 3*N+1,
 maxisuiv(N2,T,M).

maxisuiv(N,T,M):- N > M, syr(N,T,N).
maxisuiv(N,T,M):- syr(N,T,M). 

Version 2 :

syr(N):- syr(N,0,N).
syr(1,T,M):- write(T),nl,write(M).
syr(N,T,M):-
 T2 is T + 1,
 suiv(N,N2),
 M2 is max(N,M),
 syr(N2,T2,M2).

suiv(N, N2):-
 0 is N mod 2,
 N2 is N/2,!.

suiv(N, 3 * N + 1).

RPN : Recherche de permutations

On trouve, dans certains ouvrages sur les calculatrices HP, des tables donnant les combinaisons de touches pour réordonner les éléments de la pile XYZT :

Par exemple ci-dessus, ligne n°9, si ABCD sont sur la pile XYZT, la combinaison x><y, CLX et ENTER permettra d’avoir 00AC sur la pile.

L’idée m’est venue de faire la même chose pour les calculatrices HP-48/49/50g qui ont d’autres commandes pour manipuler la pile.

A partir d’éléments (e1 au niveau 1, e2 au niveau 2 etc.), on voudrait les réordonner en utilisant les commandes de base du langage RPN : DUP, SWAP, DROP, OVER, ROT, DUP2, DROP2 et PICK3.

Exemples :

La pile contient uniquement 2 éléments e1 (niveau 1) et e2 (niveau 2)et on veut les inverser. Réponse : SWAP

La pile contient 3 éléments e1 : e2 : e3 et on veut inverser les éléments des niveaux 2 et 3 pour obtenir e1 : e3 : e2. Réponse : ROT SWAP

Je vous propose un petit programme Python qui fait une recherche brute de toutes les possibilités (en partant des combinaisons à 1 commandes, puis à 2 commandes etc. Avec répétitions, par exemple OVER OVER). Il recherche jusqu’à une combinaison de 7 commandes ce qui peut demander du temps !

Tester le programme en ligne

from itertools import *

keys = "DUP", "SWAP", "DROP", "OVER", "ROT", "DUP2", "DROP2", "PICK3"

def DUP(s): return [s[0]] + s if len(s)>0 else -1
def SWAP(s): return [s[1], s[0]] + s[2:] if len(s) > 1 else -1
def DROP(s): return s[1:] if len(s) > 0 else -1
def OVER(s): return [s[1]] + s if len(s) > 1 else -1
def ROT(s): return [s[2]] + s[:2] +s[3:] if len(s) > 2 else -1
def DUP2(s): return [s[0], s[1]] + s if len(s) > 1 else -1
def DROP2(s): return s[2:] if len(s) > 2 else -1
def PICK3(s): return [s[2]] + s if len(s) > 2 else -1

def find(target, n):
 rr = [1 + v for v in range(n)]
 tr = [int(v) for v in target]
 for n in range(1, 8):
  for op in combinations_with_replacement(keys, n):
   for p in permutations(op):
    r = rr
    for c in p:
     t = eval(c)(r)
     if t == -1: break
     r = t
    if r == tr:
     rs = ':'.join([str(v) for v in rr])
     return '{} -> {} = {}'.format(rs, ':'.join(list(target)), " / ".join(p))

 return "No solution"

# 5 levels on the stack (e1:e2:e3:e4:e5). Goal : e2:e2
>>> find("22", 5)
'1:2:3:4:5 -> 2:2 = ROT / DROP2 / SWAP / ROT / DROP2 / DUP'

# 3 levels on the stack (e1:e2:e3). Goal : e1:e2:e3:e1:e2:e3
>>> find("123123", 3)
'1:2:3 -> 1:2:3:1:2:3 = PICK3 / PICK3 / PICK3'

# 3 levels on the stack. Goal : e1:e1:e2:e2:e3:e3
>>> find("112233", 3)
'1:2:3 -> 1:1:2:2:3:3 = PICK3 / ROT / ROT / DUP2 / ROT'

# 4 levels on the stack. Goal : e2:e2:e4:e4
>>> find("2244", 4)
'1:2:3:4 -> 2:2:4:4 = ROT / DROP2 / DUP2 / ROT'

Version pour la HP-50g

Cette calculatrice a plus de commandes que les HP-48G pour manipuler la pile, j’ajoute DUPDUP, UNROT et 4 ROLL (comme ROT mais avec le niveau 4)

from itertools import *

# Version HP-50g

keys = "DUP", "DUPDUP","SWAP", "DROP", "OVER", "ROT", "UNROT", "DUP2", "DROP2", "PICK3", "_4_ROLL"

def DUP(s): return [s[0]] + s if len(s)>0 else -1
def DUPDUP(s): return [s[0], s[0]] + s if len(s)>0 else -1
def SWAP(s): return [s[1], s[0]] + s[2:] if len(s) > 1 else -1
def DROP(s): return s[1:] if len(s) > 0 else -1
def OVER(s): return [s[1]] + s if len(s) > 1 else -1
def ROT(s): return [s[2]] + s[:2] +s[3:] if len(s) > 2 else -1
def UNROT(s): return s[1:3] + [s[0]] +s[3:] if len(s) > 2 else -1
def DUP2(s): return [s[0], s[1]] + s if len(s) > 1 else -1
def DROP2(s): return s[2:] if len(s) > 2 else -1
def PICK3(s): return [s[2]] + s if len(s) > 2 else -1
def _4_ROLL(s): return [s[3]] + s[:3] +s[4:] if len(s) > 3 else -1

def find(target, n):
 rr = [1 + v for v in range(n)]
 tr = [int(v) for v in target]
 for n in range(1, 8):
  for op in combinations_with_replacement(keys, n):
   for p in permutations(op):
    r = rr
    for c in p:
     t = eval(c)(r)
     if t == -1: break
     r = t
    if r == tr:
     rs = ':'.join([str(v) for v in rr])
     return '{} -> {} = {}'.format(rs, ':'.join(list(target)), " / ".join(p))

 return "No solution"

>>> find("231", 3)
'1:2:3 -> 2:3:1 = UNROT'

>>> find("4123", 4)
'1:2:3:4 -> 4:1:2:3 = _4_ROLL'

>>> find("135", 5)
'1:2:3:4:5 -> 1:3:5 = SWAP / _4_ROLL / DROP2'

>>> find("4321", 4)
'1:2:3:4 -> 4:3:2:1 = SWAP / ROT / _4_ROLL'

Version en PROLOG (en cours de réalisation)

% Règles pour les opérations sur la pile

% SWAP : Inverse les niveaux 1 et 2
swap([A,B|T], [B,A|T]).

% DUP : Duplique le niveau 1
dup([A|T], [A,A|T]).

% DROP : Supprime le niveau 1
drop([_|T], T).

% OVER : Récupère le niveau 2 et le met au niveau 1
over([A,B|T], [B,A,B|T]).

% ROT
rot([A,B,C|T], [C,A,B|T]).


% Règle pour arrêter la récursivité lorsque la pile est déjà dans l'état voulu
trouve([], Stack, Stack).

% Règles pour trouver la séquence d'opérations 

trouve([swap|Actions], Stack, FinalState) :-
    swap(Stack, NewStack),
    trouve(Actions, NewStack, FinalState).

trouve([dup|Actions], Stack, FinalState) :-
    dup(Stack, NewStack),
    trouve(Actions, NewStack, FinalState).

trouve([drop|Actions], Stack, FinalState) :-
    drop(Stack, NewStack),
    trouve(Actions, NewStack, FinalState).

trouve([over|Actions], Stack, FinalState) :-
    over(Stack, NewStack),
    trouve(Actions, NewStack, FinalState).

trouve([rot|Actions], Stack, FinalState) :-
    rot(Stack, NewStack),
    trouve(Actions, NewStack, FinalState).

Quelques tests

rot([1,4,3],X).
Réponse : X = [3, 1, 4]

swap([1,4,3],X).
Réponse : X = [4, 1, 3]

drop([1,4,3],X).
Réponse : X = [4, 3]

over([1,4,3],X).
Réponse : X = [4, 1, 4, 3]

trouve(X, [1,4,3], [4,1,3]).
X = [swap]

Version 2 proposée par chatGPT

trouve(Actions, Stack, FinalState) :-
    trouver_sequence(Actions, [], Stack, FinalState).

% Cas où l'état actuel est l'état voulu
trouver_sequence([], _, Stack, Stack).

% Cas où l'état actuel est différent de l'état voulu, continue à chercher
trouver_sequence([Action|Actions], History, CurrentState, FinalState) :-
    \+ member(CurrentState, History), % Vérifie si l'état actuel a déjà été visité pour éviter les boucles infinies
    apply_action(Action, CurrentState, NextState),
    trouver_sequence(Actions, [CurrentState|History], NextState, FinalState).

% Applique une action sur l'état actuel de la pile
apply_action(swap, [A,B|T], [B,A|T]).
apply_action(dup, [A|T], [A,A|T]).
apply_action(drop, [_|T], T).
apply_action(over, [A,B|T], [B,A,B|T]).
apply_action(rot, [A,B,C|T], [C,A,B|T]).