Après avoir longtemps cherché une réponse "mathématique" à cette énigme, je me suis tourné vers l'informatique. "Mais comment résoudre un problème avec des variables réelles (ou presque) sur un ordinateur" me direz vous : j'en donne une explication après. Après une journée entière de grattage de tête, mon programme me renvoyait une solution pour N=17 et aucune pour N=18. Alors, mais sans entière celconviction, je réponds
le N maximum demandé est 17. En revanche, voilà un exemple de suite dont je suis sûr :
0.0294 0.5420 0.8516 0.2697 0.7101 0.4226 0.9258 0.1716 0.6202
0.3431 0.7735 0.0871 0.4853 0.9706 0.2071 0.6569 0.3550
Plusieurs remarques avant d'expliquer mon algorithme :
*
Pour
fini, peu importe que les nombres cherchés soient décimaux ou juste réels. En effet, en subdivisant

en

segments

de tailles égales, on se rend compte que si

est le

-ème terme solution, et si

alors tout réel du segment

associé peut remplacer

dans la solution.

étant un ouvert non vide, et l'ensemble des décimaux étant dense dans les réels, on pourra toujours trouver déduire d'une solution "réelle" une solution "décimale".
* Ainsi, le premier terme de ma solution est 0.0294 et fait partie de l'intervalle de Farey

: donc il peut être remplacé par n'importe quel terme de cet intervalle. A l'ordre 17 qui est je l'espère l'ordre maximum,
le nombre de suites solutions est infini (et dénombrable je pense).
Mon algorithme :
* Avec ma première remarque, j'ai rendu le nombre de variables fini (les

). Cependant, travailler avec N! est un peu lourd : je me suis donc servi des intervalles de Farey, qui sont eux au nombre de

.
Je travaille donc ensuite non pas sur l'ensemble des nombres décimaux ou réels de

, mais seulement sur l'ensemble des intervalles de Farey, ou plutôt
sur des nombres décimaux représentatifs de chaque intervalle de Farey (car si

est un nombre de l'intervalle de Farey n°

alors, tous les nombres de cet intervalle peuvent remplacer

).
* Etant donné les

premiers termes de la suite solution, un nombre

ne peut être le

-ème terme que si :
Pour tout

(et

bien sûr) et pour tout

:

n'est pas dans le même

-ème intervalle que

.
C'est une condition nécessaire mais non suffisante bien sûr. Cela se traduit informatiquement par la condition :
\neq int(j*u_i))
, avec
)
= partie_entiere de X.
* On peut également remarquer qu'il suffit de vérifier la condition précédente non pas sur tous les

, mais seuleument pour tous les

puis pour tous les

premiers

.
* Avec ces quelques remarques, on cherche alors par "exploration des possibilités". Etant donné qu'il y a de l'ordre de

nombres candidats pour la solution d'ordre

, le calcul s'effectue assez rapidement.
Voilà désolé pour la longueur du message, mais je tenais à expliquer mon raisonnement, car
celui qui a passé du temps sur cette énigme ne peut passer sous silence sa logique. Un bravo pour ceux qui auront démontré (rigoureusement, pas comme moi) le résultat. Et un merci à Minkus pour ce genre d'énigme comme j'aimerais en avoir plus souvent !