Calcul des valeurs de Grundy pour des jeux octaux

De Wiki du LAMA (UMR 5127)
Aller à la navigation Aller à la recherche

Étudiant : Mathieu BRUNOT

Tuteur : Valenti Gledel

Jeu de Impartiaux

Description

Un jeu impartial est un jeu à 2 joueurs où les joueurs jouent à tour de rôle et leurs mouvements possibles ne dépendent pas du joueur qui joue, il n'y a donc pas de hasard ni d'égalité dans une partie.

Le perdant étant celui qui ne peut plus jouer, aussi appelé version normale.

En voici des exemples:

  • Jeu de Cram
    Les joueurs mettent des dominos sur une grille. Le joueur qui met le dernier domino gagne.
    Jeu cram.png
  • Jeu de Nim
    Il se joue avec un ou plusieurs tas d'objets exemple des allumettes.
    Chacun des joueurs joue à tour de rôle en choisissant un tas et en retirant autant d’allumettes que souhaité dans le tas, pour que son adversaire n’ait plus la possibilité de jouer, c’est ainsi que le jeu prend fin.
    Jeu de nim.jpg
  • Jeu des bâtonnets Fort Boyard
    Ce jeu consiste à n’avoir que le droit de retirer de 1 à 3 bâtonnets à chaque tour. Hormis la différence est que le joueur perdant est celui qui prend le dernier bâtonnet. Cette condition de victoire est appelé version misère.
    Batonnets-fort-boyard.jpg
  • Stratégie Gagnante

    En effet, l’ensemble des jeux impartiaux sont résolubles à l’aide d’outils théorique afin de connaitre l'action à faire pour être le joueur en position gagnante.

    Prenons le jeu de Fort Boyard mais cette fois-ci en version normale. Dans celui-ci la stratégie gagnante est de laisser un tas à l'adversaire qui serait un multiple de quatre.

    • Pour le trouver il suffit de regarder les positions perdantes dans chaque cas:
      1. Si le nombre de bâtonnets est n=1, le joueur courant gagne car il peut enlever le dernier jeton et laisser l'adversaire avec un tas vide.
      2. Si n=2, le joueur courant gagne car il peut retirer les derniers bâtonnets.
      3. Si n=3, le joueur courant gagne car il peut retirer les derniers bâtonnets.
      4. Si n=4, le joueur courant perd car il va forcément laisser son adversaire dans une position gagnante.
      5. Plus généralement, toute position où est une position gagnante, et toute position où est une position perdante.

    Ainsi on pourrait se demander la stratégie gagnant d'autres jeux tel que le jeu de Nim ou le jeu de Cram

    Valeurs de Grundy

    Théorème de Sprague-Grundy

    Le théorème de Sprague-Grundy nous permet de trouver la stratégie gagnante d’un jeu impartial fini en version normale (le perdant est le joueur ne pouvant plus jouer).

    La valeur de Grundy ou d’un tas se calcule récursivement par ses règles :

    - Un tas vide a pour valeur de Grundy 0.
    - Le nombre Grundy d’une position non nulle est le plus petit entier positif ou nul qui n’est pas dans la liste des valeurs de Grundy des positions que l’on peut atteindre à partir de notre position.
    

    Une position perdante d’un jeu celle qui donne une valeur de Grundy de 0.

    Calcul des valeurs pour le jeu de Nim

    Cas à 1 tas

    La valeur de Grundy du jeu de Nim classique d’un tas de taille n (un entier) est celle du nombre de bâtonnets dans ce tas.

    On l’explique car toutes les tailles de tas strictement inférieures à n sont atteignables à partir de celui-ci. Alors le minimum exclu(mex) des valeurs atteignables est n.

    On appelle ce nombre précis dans le cas du jeu de Nim, un nimber.

    Donc par exemple dans un tas de taille 9 du jeu Nim, on aurait une valeur de Grundy de 9, ceci est simple à comprendre.

    Cas à plusieurs tas

    Dans le cas de plusieurs tas, la valeur de Grundy de tous les jeux est le XOR de toutes les valeurs de Grundy de chaque tas.

    Par exemple :

    Voici un jeu de Nim composé de trois tas de bâtonnets: Exemple jeu de nim.png

    Ce jeu a pour valeur de Grundy :

    Application

    Mettons à l’œuvre nos connaissances afin de gagner une partie dans le jeu de Nim de l'exemple précédent.

    Imaginons que nous soyons les premiers joueurs à jouer.

    Nous nous trouvons dans une position gagnante car la valeur de Grundy du Jeu est de .

    Pour pouvoir ramener celui-ci à une valeur de 0 pour que notre adversaire soit dans une position perdante, il nous suffit de changer un tas où la valeur de grundy de celui-ci possède un bit à 1 au même indice que le bit à 1 avec le plus grand exposant de 2 de la valeur de Grundy du jeu.

    Dans notre cas le premier bit à 1 de 15 est celui qui est au 3eme exposant de 2. Le nombre qu’on doit changer est 24. Pour savoir à quelle valeur on doit le baisser il suffit de faire le XOR de 24 et de la valeur de Grundy du jeu qui est 15. On obtient .

    • On a choisi de retirer 1 bâton au tas de 24 pour avoir une valeur de Grundy totale de 0.
    • L'adversaire est dans une position perdante, quoiqu’il fasse il ne peut pas passer en position gagnante.
    • Il nous reste plus qu'à continuer ainsi pour s'assurer de gagner.

    Jeux Octaux

    Les jeux octaux sont aussi des jeux impartiaux, assez proches du jeu de Nim. Ils se jouent aussi avec des tas d’objet. Leurs caractéristiques sont que les coups autorisés sont définis à l’aide d’un système numérique de 0 à 7 d’où le nom jeu octal.

    Règles

    Les règles d'un jeu octal précisent les actions possibles lorsqu'un joueur retire n objets d'un tas donné, et sont encodées par une suite de chiffres entre 0 et 7 (d'où le nom jeu octal) :

    Notation 0.k1 k2 k3 k4

    où chaque chiffre kn indique en combien de tas le joueur doit scinder le tas dans lequel il a retiré n objets. Le chiffre kn est la somme de :

    • si le joueur peut laisser 0 tas, 0 sinon;
    • si le joueur peut laisser 1 tas, 0 sinon;
    • si le joueur peut laisser 2 tas, 0 sinon

    Par exemple, si kn=0, les joueurs n'ont pas le droit de retirer n objets d'un tas. Si kn=3, on retrouve la règle classique du jeu de Nim, à savoir que le joueur peut retirer n objets d'un tas de taille supérieure ou égale n.

    Représentation en jeu octal

    Comme vous l'avez compris il y a une infinité de jeux octaux.

    On retrouve les jeux qu'on a introduit dans ces catégories. Jeu de Fort Boyard = 0.333
    Jeu de Fort Cram = 0.07

    Néanmoins, le jeu de Nim n'est pas un jeu octal. Dans Nim, les mouvements sont extrêmement flexibles : un joueur peut retirer n'importe quel nombre d'objets d'une pile choisie. Cette flexibilité n'est pas présente dans les jeux octaux, où les mouvements sont prédéfinis. Un jeu "octal" qui pourrait représenter le jeu de Nim serait 0.333... un sauf que ce n'est pas un jeu octal.

    Valeurs de Grundy pour les jeux octaux

    Méthode

    JEU_VALIDE = (0,1,2,3,4,5,6,7)
    TAS_0 = (1,3,5,7) #les valeurs laissant 0 tas
    TAS_1 = (2,3,6,7) #les valeurs laissant 1 tas
    TAS_2 = (4,5,6,7) #les valeurs laissant 2 tas
    
    def verif_jeu(nbr_octal:str)->str:
        """ renvoie nombre en str et enleve le 0. """
        assert type(nbr_octal) == str, "tu dois donner une chaine de caractères en argument"
        for ele in nbr_octal:
            ele = int(ele)
            assert (ele in JEU_VALIDE), "ton jeu octal n'est pas conforme"
        return nbr_octal
    
    def calc_mex(un_set:set)->int:
        """ renvoie la plus petite valeur entière qui n'appartient pas au subset. """
        mex = 0
        while (mex in un_set):
            mex += 1
        return mex
    
    def tab_kn_dn(nbr_octal:str)->list:
        """ renvoie un tableau de tableau avec [kn:str,dn:int] pour chaque action d'un jeu octal. """
        nb = verif_jeu(nbr_octal)
        tab_res = []
        
        for i in range(0, len(nb)):
            tab_res += [[nb[i], i+1]]
        
        return tab_res
    
    def liste_action_possibles(taille_tab:int, kn:str, dn:int)->list:
        """ renvoie un tableau de bool sur la possibilité des 3 actions. """
        assert taille_tab > 0 , 'ton tableau est vide'
        
        k_int = int(kn)
        action = [False,False,False] #tableau des 3 actions est nous dis si elles sont possibles
        if not dn > taille_tab: #lance les conditions seuleument la taille retiré est convenable
            if k_int in TAS_0:
                if dn == taille_tab:
                    action[0] = True
            if k_int in TAS_1:
                if dn < taille_tab:
                    action[1] = True
            if k_int in TAS_2:
                if dn+1 < taille_tab:
                    action[2] = True
        
        return action
        
    def liste_des_separations_tas_en_2(taille_tab:int, dn:int)->list:
        """ renvoie tous les couples de possibilités d'une separation en deux d'un tas
        dans un tableau de set. """
        res = []
        somme_des_taille_couple = taille_tab - dn
        for i in range(1, (somme_des_taille_couple//2)+1):
            val1 = i
            val2 = somme_des_taille_couple-i
            couple = (val1,val2) #crée le couple sous forme de tuple
            res += [couple]
        return res
    
    def set_des_valeurs_grundy_atteignables(jeu_actions:list, tab_grundy:list, instant_i:int):
        """ renvoie un set des valeurs atteignables pour le coup pour calculer la 
        valeur suivant du tab_grundy"""
        res = set()
        for actions in jeu_actions:
            les_actions_possibles = liste_action_possibles(instant_i, actions[0], actions[1])
            if les_actions_possibles[0] == True:
                res.add(0)
                
            if les_actions_possibles[1] == True:
                res.add(tab_grundy[instant_i - actions[1]])           
                    
            if les_actions_possibles[2] == True:
                liste_separations_possibles = liste_des_separations_tas_en_2(instant_i, actions[1])
                for tab in liste_separations_possibles:
                    res.add(tab_grundy[tab[0]] ^ tab_grundy[tab[1]])
        return res
    
    def tab_val_grundy(taille_tab:int, nbr_octal:str)->list:
        """ renvoie un tableau de valeur de Sprague Grundy pour une taille_tab donné. """
        tab_grundy = [0]
        jeu_actions = tab_kn_dn(nbr_octal)
        for i in range(1, taille_tab+1):
            set_valeurs_atteignables = set_des_valeurs_grundy_atteignables(jeu_actions, tab_grundy, i)              
            tab_grundy += [calc_mex(set_valeurs_atteignables)]
        return tab_grundy
        
    def Grundy(taille_tab:int , nbr_octal:float):
        """ renvoie la valeur de grundy du nombre donné"""
        tab_grundy = tab_val_grundy(taille_tab, nbr_octal)
        return tab_grundy[taille_tab]
    

    Périodicité

    def est_periode(p:int,s:int,tab_grundy:list,n:int)->bool:
        """ test la périodicité"""
        res = True
        m = n-(p+s)
        for i in range(m):
            if s+p+i == len(tab_grundy):
                res = False
                break
            if tab_grundy[s+i] != tab_grundy[s+p+i]:
                res = False
                break
        return res
    
    def test_periodicite(n:int,k:int,tab:list)->list:
        """ renvoie un tableau si on trouve la période sinon renvoie None"""
        for p in range(1,n+1):#periode
            s = n-p
            if est_periode(p,s,tab,n*2+k):
                return (s,p,tab[s:s+p])
        return None
    
    def periodicite(nbr_octal:str)->list:
        """ renvoie la période , taille de la période, taille de la prépériode d'un tableau avec les valeurs de grundy d'un jeu"""
        taille_test = 30000
        tab_test = tab_val_grundy(taille_test, nbr_octal)[1:taille_test] #tableau où chercher la périodicité
        k = len(nbr_octal)
        for n in range(k+2,len(tab_test)+1,2):
            res = test_periodicite(n,k,tab_test)
            if res != None:
                return res 
        return "Periode non trouvé pour .{}".format(nbr_octal)
    

    Méthode rare/fréquentes

    Masque

    def trie_dico_en_tab_croissant(dico:dict)->list:
        """ renvoie un tableau de tableau trié par ordre croissant 
        selon les valeurs des clefs"""
        res = []
        dico_copy = dico.copy()
        while dico_copy != {}:
            clef_min = min(dico_copy, key=dico_copy.get)#cherche la clef avec la val minimale
            res += [[clef_min,dico_copy[clef_min]]]
            dico_copy.pop(clef_min)
        return res
            
    def cherche_id_gap_rare_freq(tab_frequences:list, nbr_echantillon:int)->int:
        """ renvoie l'id du dernier element rare
        d'un tableau trie de frequence de valeur """
        freq = 0.1 * nbr_echantillon #frequence choisi pour determiner les valeurs rares
        for i in range(len(tab_frequences)):
            if tab_frequences[i] > freq:
                id_gap = i-1
                break
        return id_gap
    
    def cherche_rare_commun(tab_grundy:list)->list:
        """ cherche les valeures rare d'un tableau en regardant leurs nombre d'apparitions
        à partir d'un dico qui a pour clef le nombre et pour valeur son nombre d'apparition"""
        nombre_ele = len(tab_grundy)
        dico_frequence = {}
        for ele in tab_grundy:
            dico_frequence[ele] = dico_frequence.get(ele,0) + 1 #creer le dico de toutes les valeurs en clef et en valeurs leurs fréquences
        
        tab_trie = trie_dico_en_tab_croissant(dico_frequence)
        print(tab_trie)
        tab_clef = [tab[0] for tab in tab_trie]
        tab_frequences = [tab[1] for tab in tab_trie]
        id_gap = cherche_id_gap_rare_freq(tab_frequences,nombre_ele)
        
        rare = {0}
        commun = set()
        
        for i in range(len(tab_clef)):
            if i <= id_gap:
                rare.add(tab_clef[i])
            else:
                if tab_clef[i] != 0:
                    commun.add(tab_clef[i])
                
        
        res = [rare,commun]
        return res
    
    def indice_des_bits_a_un(masque:str)->list:
        """ renvoie un tableau contenant les indices des bits à 1"""
        res = []
        for i in range(len(masque)):
            if masque[i] == '1':
                res.append(i)
        return res
    
    def est_rare(ele_bin:str,ind_1_du_masque:list)->bool:
        """ verfie si le nombre est rare"""
        cpt_1 = 0
        for ind in ind_1_du_masque:
            if ele_bin[ind] == '1':#verifie que l'element à l'indice ind est bien 1 à l'indice des 1 sur le masque
                cpt_1 += 1
        return cpt_1 % 2 == 0
    
    def rares_dans_masque(masque:str,val_max:int)->set:
        """Donne toutes les valeurs rares d'un masque dans un set"""
        res = set()
        n_bit = len(masque)
        ind_des_1 = indice_des_bits_a_un(masque)
        for nbr in range(val_max+1):
            nbr_bin = mets_nombre_en_bin_sur_x_bits(nbr,n_bit)
            if est_rare(nbr_bin,ind_des_1):
                res.add(nbr)
        
        return res
    
    def cherche_masque_pour_un_tab_grundy(tab_grundy:list)->str:
        """ va chercher un masque de n bit qui convient le mieux pour le tableau de grundy donné"""
        tab_grundy = tab_grundy[1:] #enlève la première valeur qui est 0
        val_max = max(tab_grundy)
        taille_bit = len(bin(val_max)[2:])
        masques_possibles = 2**taille_bit
        tab_rare_commun= cherche_rare_commun(tab_grundy)
        set_rare_cherche = tab_rare_commun[0]
        set_frequent = tab_rare_commun[1]
        res = 0
        
        print("\nset_rare : {}".format(set_rare_cherche))
        print("set_frequent : {}\n".format(set_frequent))
        
        for i in range(1,masques_possibles):
            masque = mets_nombre_en_bin_sur_x_bits(i,taille_bit)
            set_rare_masque = rares_dans_masque(masque,val_max)
            
            contient_rare = set_rare_cherche.issubset(set_rare_masque)
            inter = set_frequent.intersection(set_rare_masque)
            ne_contient_pas_frequent = inter == set() 
            if contient_rare and ne_contient_pas_frequent: #verifie si les rares cherchés sont biens dans le masque et qu'il ne contient pas de valeurs frequentes
                res = masque
                break
        return f"\nle masque qui convient est {res}"
    

    Code

    def liste_des_id_des_rares_dans_tab(tab_grundy:list, listes_rares:set):
        """ renvoie un tableau de tous les id des valeurs rares dans le tableau """
        res = []
        for i in range(1,len(tab_grundy)):
            if tab_grundy[i] in listes_rares:
                res + [i]
        return res
    
    def liste_des_separations_tas_en_2_pour_avoir_val_commune(taille_tab:int, dn:int, liste_rares:set)->list:
        """ renvoie tous les couples de possibilités d'une separation en deux d'un tas
        contenant uniquement une valeur rare"""
        res = []
        somme_des_taille_couple = taille_tab - dn
        
        for val1_rare in liste_rares:
            val2 = somme_des_taille_couple-val1_rare
            if not val2 in liste_rares: #verifie que la valeur ne soit pas rares
                couple = (val1_rare,val2) #crée le couple sous forme de tuple
            res += [couple]
        
        return res
    
    def set_des_valeurs_grundy_atteignables_rares(jeu_actions:list, tab_grundy:list, instant_i:int, liste_rares:set):
        """ renvoie un set des valeurs atteignables avec des couples contenant des valeurs rares"""
        res = set()
        for actions in jeu_actions:
            les_actions_possibles = liste_action_possibles(instant_i, actions[0], actions[1])
            if les_actions_possibles[0] == True:
                res.add(0)
                
            if les_actions_possibles[1] == True:
                res.add(tab_grundy[instant_i - actions[1]])           
                    
            if les_actions_possibles[2] == True:
                global liste_separations_possibles_rares_commun
                liste_separations_possibles_rares_commun = liste_des_separations_tas_en_2_pour_avoir_val_commune(instant_i, actions[1], liste_rares)
                for couple in liste_separations_possibles_rares_commun:
                    res.add(tab_grundy[couple[0]] ^ tab_grundy[couple[1]])
        return res
    
    def tab_val_grundy_methode_masque(taille_tab_souhaite:int, nbr_octal:str, masque:str)->list:
        taille_tab_inital = 10000
        assert taille_tab_souhaite > taille_tab_inital, f"La taille de tableau recherché doit dépasser {taille_tab_inital}"
        tab_grundy = tab_val_grundy(taille_tab_inital, nbr_octal)
        jeu_actions = tab_kn_dn(nbr_octal)
        listes_rares = rares_dans_masque(masque,max(tab_grundy)) #un set contant toutes les valeurs rares du masque
        liste_id_rares = liste_des_id_des_rares_dans_tab(tab_grundy, listes_rares)
    
        for taille in range(taille_tab_inital+1, taille_tab_souhaite+1):#Trouver les valeurs suivantes avec la méthode du masque
            set_valeurs_atteignables = set_des_valeurs_grundy_atteignables_rares(jeu_actions, tab_grundy, taille, liste_id_rares) 
            mex = calc_mex(set_valeurs_atteignables)
            if not mex in listes_rares:#si le mex est commun alors
                tab_grundy += [mex]
            
            else: #on va rajouter des valeurs jusqu'à avoir un mex commun sinon on continue jusqu'au bout
                for actions in jeu_actions:
                    les_actions_possibles = liste_action_possibles(taille, actions[0], actions[1])
                    if les_actions_possibles[2] == True:
                        somme_des_taille_couple = taille - actions[1]
                        for i in range(1, (somme_des_taille_couple//2)+1):
                            val1 = i
                            val2 = somme_des_taille_couple-i
                            if all(val1 not in couple for (couple) in liste_separations_possibles_rares_commun):#verifie que la valeur n'a pas déjà été testé
                                set_valeurs_atteignables.add(tab_grundy[val1] ^ tab_grundy[val2])
                                mex = calc_mex(set_valeurs_atteignables)
                                if not mex in listes_rares:
                                    break
                    if not mex in listes_rares:
                        break
                tab_grundy += [mex]
                
        return tab_grundy