La (presque) question du mercredi

Ici, on fait dans le petit, le LCD qui déchire sa race, on y cause même calculatrices quand on est en manque !

Modérateur : Politburo

Répondre
Gilles59
Fonctionne à 2400 bauds
Fonctionne à 2400 bauds
Messages : 1371
Inscription : 27 oct. 2010 20:46

La (presque) question du mercredi

Message par Gilles59 » 08 mai 2019 14:56

Je dis 'presque' parce que je n'ai moi-même pas la réponse…

Existe-il un moyen optimisé (en vitesse) applicable sur une calcultrice pour trouver la suite des nombres naturels 't' tels que le résultat suivant soit un nombre entier :

Image

ce qui vient naturellement ) l'esprit c'est de tester les valeur de 't' avec une boucle mais ca devient vite très inefficace… Par exemple HP50

Code : Tout sélectionner

«
  -> t
  «
    DO 
     IF 'NOT(FP((1+√(2*t^2-2*t+1))/2)' THEN t 1. DISP END
     1 't' STO+ t 2. DISP
    UNTIL 0 END
  »
»
Casio FX-502P /602P / 603P / FX180P+ / FX4000P / TI57 / TI66 / TI74 Basicalc / TI95 Procalc / HP12C / HP15C LE / DM41L / HP 30B / HP39GII / HP 48SX USA / 49G / 49G+ / 50G / 50G NewRPL / HP Prime / Oric 1 / Amstrad CPC 6128+

Avatar de l’utilisateur
gege
Fonctionne à 14400 bauds
Fonctionne à 14400 bauds
Messages : 7016
Inscription : 31 janv. 2008 15:24
Localisation : Banlieue Paârisienne
Contact :

Re: La (presque) question du mercredi

Message par gege » 08 mai 2019 17:08

Bonjour,
En regardant super vite, j'ai l'impression qu'il n'y a pas de solutions ???
Je vérifie.
G.E.

EDIT : pff planté

Gilles59
Fonctionne à 2400 bauds
Fonctionne à 2400 bauds
Messages : 1371
Inscription : 27 oct. 2010 20:46

Re: La (presque) question du mercredi

Message par Gilles59 » 08 mai 2019 17:29

En fait l'idée est que le résultat soit un carré parfait. Les premiers 't' possibles sont :

{ 1 4 21 120 697 4060 }
Casio FX-502P /602P / 603P / FX180P+ / FX4000P / TI57 / TI66 / TI74 Basicalc / TI95 Procalc / HP12C / HP15C LE / DM41L / HP 30B / HP39GII / HP 48SX USA / 49G / 49G+ / 50G / 50G NewRPL / HP Prime / Oric 1 / Amstrad CPC 6128+

Avatar de l’utilisateur
Marge
Fonctionne à 14400 bauds
Fonctionne à 14400 bauds
Messages : 5278
Inscription : 01 oct. 2008 14:39
Localisation : En bas, tout au fond à gauche.

Re: La (presque) question du mercredi

Message par Marge » 08 mai 2019 17:51

f(4)=3 pas carré...
3 hommes, 3 demis, un 3a... Magnéto, Serge !

« Boris », c'est juste Maurice enrhumé.

Avatar de l’utilisateur
zpalm
Fonctionne à 9600 bauds
Fonctionne à 9600 bauds
Messages : 2498
Inscription : 03 mai 2008 15:33
Localisation : Grenoble

Re: La (presque) question du mercredi

Message par zpalm » 08 mai 2019 18:32

Marge a écrit :
08 mai 2019 17:51
f(4)=3 pas carré...
Oui, mais 3 est un entier, c'est ce qu'on veut.

2*t^2 -2*t+1 doit être un carré entier
donc si: 2*t^2 -2*t+1=n^2, alors: t= (1/2)*(sqrt(2*n^2-1)+1)
de plus (1+n)/2 est entier donc n est impair.

Le programme ci-dessous pour HP Prime donne les 9 premières valeurs de t:

Code : Tout sélectionner

EXPORT QDM()
BEGIN
 LOCAL n,t;
 FOR n FROM 1 TO 1000000 STEP 2 DO
  t:= (1/2)*(sqrt(2*n^2-1)+1);
  IF FP(t)==0 THEN
   PRINT("t="+t+" , f="+(1+n)/2);
  END;
 END;
 PRINT("Terminé");
END;

Avatar de l’utilisateur
bernouilli92
Fonctionne à 9600 bauds
Fonctionne à 9600 bauds
Messages : 4307
Inscription : 21 nov. 2012 14:03
Localisation : Ile de France

Re: La (presque) question du mercredi

Message par bernouilli92 » 08 mai 2019 19:26

Voici ma version pour hp 48 :

Code : Tout sélectionner

« { }
  1 35000 FOR N 
  N DUP 1 DISP 
  SQ 1 - 2 / √ CEIL
  IF DUP DUP 1 - * 2 * 1 + N SQ ==
  THEN
    DUP 2 DISP +
  ELSE 
    DROP
  END
  2 STEP 
»
Il faut trouver t entier tel qu'on ait : t(t-1) = (N^2-1)/2 avec N entier.
Comme N^2 doit être impair, N est aussi impair. On peut donc tester que les N impairs.

Sur émulateur, cela tourne beaucoup plus vite que sur une vraie hp48.
Avec le programme précédent, je trouve (en 5 secondes) { 0 4 21 120 697 4060 23661 }

0 est aussi un résultat possible.
HP, Casio, Sharp, Psion, quelques TI et divers autres

Gilles59
Fonctionne à 2400 bauds
Fonctionne à 2400 bauds
Messages : 1371
Inscription : 27 oct. 2010 20:46

Re: La (presque) question du mercredi

Message par Gilles59 » 08 mai 2019 19:52

Gilles59 a écrit :
08 mai 2019 17:29
En fait l'idée est que le résultat soit un entier. Les premiers 't' possibles sont :

{ 1 4 21 120 697 4060 }
Oups j’utilise mon smartphone et cafouille un peu. Je voulais éditer mon message. Et Bien ça calcule déjà deux fois plus vite. J’ai l’intuition qu’il doit y avoir un moyen plus rapide mais ne le trouve pas.

Édit : zpalm j’en n’ai pas testé ta solution mais ça m’a L’air t’es efficace en gain de temps.
Casio FX-502P /602P / 603P / FX180P+ / FX4000P / TI57 / TI66 / TI74 Basicalc / TI95 Procalc / HP12C / HP15C LE / DM41L / HP 30B / HP39GII / HP 48SX USA / 49G / 49G+ / 50G / 50G NewRPL / HP Prime / Oric 1 / Amstrad CPC 6128+

Avatar de l’utilisateur
dprtl
Fonctionne à 1200 bauds
Fonctionne à 1200 bauds
Messages : 388
Inscription : 27 janv. 2013 01:26
Localisation : Strasbourg
Contact :

Re: La (presque) question du mercredi

Message par dprtl » 08 mai 2019 20:05

Vous êtes trop rapides les gars ! Je dois me dépêcher de publier ma version MPOïsée pour Casio PB-1000 qui affiche très rapidement la suite entière de Gilles59, jusqu'au maximum de chiffres significatifs pour 't' sur cette calculette, soit 5406093004 :

Code : Tout sélectionner

10 A0=1:PRINT A0:A1=4:PRINT A1
20 A2=6*A1-A0-2:A0=A1:A1=A2
30 IF A2<1E10 THEN PRINT A2:GOTO 20

Avatar de l’utilisateur
C.Ret
Fonctionne à 2400 bauds
Fonctionne à 2400 bauds
Messages : 2320
Inscription : 31 mai 2008 23:43
Localisation : N 49°22 E 6°10

Re: La (presque) question du mercredi

Message par C.Ret » 08 mai 2019 20:10

Ah! Zut, je rentre du boulot et voilà j'arrive trop tard pour proposer ma solution.

Bon, c'est clair, boucler sur les carrés impairs est un progrès par rapport à la recherche séquentielle sur la variable t.

Mais c'est toujours de la force brute. Or je constate que les solutions sont de plus en plus éloignées les unes des autres, n'y aurait-il pas un moyen de ne pas 'parcourir' tous les n, mais une fraction seulement, quitte à tâtonner un peu autour des solutions possibles ??

C'est (presque) une question, je n'ai pas de solution (pour l'instant), uniquement un fin rayon jaune qui annonce l'aube d'une conjecture…


EDIT: zut je suis grillé, pas assez rapide
SHARP PC-1211 + CE-121 + CE-122. | VIC 20 Commodore 128D + Printer P-803 + SD2iec. | TI-57 LCD | TI-74 BasiCalc | TI-92 II | HP-28S + HP82240A | HP-41C + (2 memory + stat + IR) modules. | HP Prime Wireless Graphing Calculator | HP-15C | CASIO fx-602p + FA-1. .Sommaire des M.P.O.. . Sommaire du P.C.T.M. .

Gilles59
Fonctionne à 2400 bauds
Fonctionne à 2400 bauds
Messages : 1371
Inscription : 27 oct. 2010 20:46

Re: La (presque) question du mercredi

Message par Gilles59 » 08 mai 2019 20:21

dprtl a écrit :
08 mai 2019 20:05
Vous êtes trop rapides les gars ! Je dois me dépêcher de publier ma version MPOïsée pour Casio PB-1000 qui affiche très rapidement la suite entière de Gilles59, jusqu'au maximum de chiffres significatifs pour 't' sur cette calculette, soit 5406093004 :

Code : Tout sélectionner

10 A0=1:PRINT A0:A1=4:PRINT A1
20 A2=6*A1-A0-2:A0=A1:A1=A2
30 IF A2<1E10 THEN PRINT A2:GOTO 20
je ne comprends rien à ta solution mais elle semble parfaite... Et très RPN friendly à priori !

ca marche nickel ! Mais pourquoi ?

EDIT : L'émulateur NewRPL me sort "{ 21 120 697 4060 23661 137904 803761 4684660 27304197 159140520 927538921 5406093004 31509019101
183648021600 1070379110497 6238626641380 36361380737781 211929657785304 }" instanément. L'émulateur HP50 rom HP est aussi rapide en mode exact.

Le programme n'est pas optimisé en plus :
«
{} 1 4 0 ->  A0 A1 A2
« DO 6 A1 * A0 - 2 - 'A2' STO A1 'A0' STO A2 'A1' STO A2 + UNTIL 'A2>1E15' END »
»
Casio FX-502P /602P / 603P / FX180P+ / FX4000P / TI57 / TI66 / TI74 Basicalc / TI95 Procalc / HP12C / HP15C LE / DM41L / HP 30B / HP39GII / HP 48SX USA / 49G / 49G+ / 50G / 50G NewRPL / HP Prime / Oric 1 / Amstrad CPC 6128+

Avatar de l’utilisateur
C.Ret
Fonctionne à 2400 bauds
Fonctionne à 2400 bauds
Messages : 2320
Inscription : 31 mai 2008 23:43
Localisation : N 49°22 E 6°10

Re: La (presque) question du mercredi

Message par C.Ret » 08 mai 2019 21:19

Pourquoi ça marche ?

Vu la forme de la solution, cela doit pouvoir se démontrer par récurrence. En partant des hypothèses que n0 et n1 sont deux impairs dont les carrés sont exacts pour t0 et t1, démontrer que n2 sera un carré exact pour t2 = 6.t1+t0+2


Sachant que l'on a:
n0 = (1+SQRT(2.t0.(t0-1)+1))/2 soit (2.n0-1)² = 2.t0.(t0-1)+1
n1 = (1+SQRT(2.t1.(t1+1)+1))/2 soit (2.n1-1)² = 2.t1.(t1-1)+1

Montrer que
(1+SQRT(2.t2.(t2+1)+1))/2 avec t2=6.t1-t0-2 est aussi une racine entière c'est-à-dire que 2.t2.(t2-1)+1 est aussi le carré d'un entier impair.

[…]

La vraie question, est comment avoir aussi vite trouvé la conjecture 6.t1-t0-2 ???
Est-ce un cas d'école des suites ?

Bon, je suis trop fatigué ce soir: j'ai juste un petit code RPN pour HP-28S qui crée la liste des solutions t :

Code : Tout sélectionner

« { 1 4 } DUP LIST→ DROP	// Initialise la suite à partir du début la liste
  DO				// Boucle principale:
    DUP 6 * ROT 2 + -		//   sauvegarde t1 et calcule t2 ← 6.t1-t0-2 en détruisant t0 
    ROT OVER + ROT ROT		//   Ajoute t2 en fin de liste
 UNTIL DUP  1E11 > END		// Boucle jusqu'à précision maximale
 DROP2 »			// Retire t1 et t0 de la pile, laisse la liste des solutions

Je trouve
{ 1 4 21 120 697 4060 23661 137904 803761 4684660 27304197 159140520 927538921 5406093004 31509019101 183648021600 }
sans pouvoir vérifier la validité des dernières solutions sans utiliser un autre calculateur.
Dernière édition par C.Ret le 08 mai 2019 21:46, édité 1 fois.
SHARP PC-1211 + CE-121 + CE-122. | VIC 20 Commodore 128D + Printer P-803 + SD2iec. | TI-57 LCD | TI-74 BasiCalc | TI-92 II | HP-28S + HP82240A | HP-41C + (2 memory + stat + IR) modules. | HP Prime Wireless Graphing Calculator | HP-15C | CASIO fx-602p + FA-1. .Sommaire des M.P.O.. . Sommaire du P.C.T.M. .

Gilles59
Fonctionne à 2400 bauds
Fonctionne à 2400 bauds
Messages : 1371
Inscription : 27 oct. 2010 20:46

Re: La (presque) question du mercredi

Message par Gilles59 » 08 mai 2019 21:40

DPRL a réussi à MPOiser ce problème… :O
Je sors ma 602P LOL

Ca donne pour les 501 à 603P (testé sur 603P). 18 pas

Code : Tout sélectionner

1 Min00 4 Min01
LBL1
 6 x MR01 - MR00 - 2 = Min02
 X<->M01 Min00 HLT
GOTO1
Pas d'astuce. L'utilisation de X<->M qui échange un registre d'affichage et une mémoire est intéressante ici
Dernière édition par Gilles59 le 08 mai 2019 22:25, édité 3 fois.
Casio FX-502P /602P / 603P / FX180P+ / FX4000P / TI57 / TI66 / TI74 Basicalc / TI95 Procalc / HP12C / HP15C LE / DM41L / HP 30B / HP39GII / HP 48SX USA / 49G / 49G+ / 50G / 50G NewRPL / HP Prime / Oric 1 / Amstrad CPC 6128+

Gilles59
Fonctionne à 2400 bauds
Fonctionne à 2400 bauds
Messages : 1371
Inscription : 27 oct. 2010 20:46

Re: La (presque) question du mercredi

Message par Gilles59 » 08 mai 2019 21:42

C.Ret a écrit :
08 mai 2019 21:19
Je trouve
{ 1 4 21 120 697 4060 23661 137904 803761 4684660 27304197 159140520 927538921 5406093004 31509019101 183648021600 }
sans pouvoir vérifier la validité des dernières solutions sans utiliser un autre calculateur.
Ma 50g me confirme en mode exact que tes solutions sont valides ;D
Casio FX-502P /602P / 603P / FX180P+ / FX4000P / TI57 / TI66 / TI74 Basicalc / TI95 Procalc / HP12C / HP15C LE / DM41L / HP 30B / HP39GII / HP 48SX USA / 49G / 49G+ / 50G / 50G NewRPL / HP Prime / Oric 1 / Amstrad CPC 6128+

Gilles59
Fonctionne à 2400 bauds
Fonctionne à 2400 bauds
Messages : 1371
Inscription : 27 oct. 2010 20:46

Re: La (presque) question du mercredi

Message par Gilles59 » 08 mai 2019 23:04

Allez c'est l'occasion de sortir la TI57 ;D

La programmation est ici absolument identique à la Casio 602P, pas par pas.
Et bien je dois dire que l'affichage en diodes rouges à son charme ;D
Ce qui est amusant, c'est que même sur un programme aussi simple, on la voit calculer ...

Code : Tout sélectionner

1
STO0
4
STO1
Lbl1
6
x
RCL1
-
RCL0
-
2
=
STO2
Exc1
STO0
RS
GTO1
Casio FX-502P /602P / 603P / FX180P+ / FX4000P / TI57 / TI66 / TI74 Basicalc / TI95 Procalc / HP12C / HP15C LE / DM41L / HP 30B / HP39GII / HP 48SX USA / 49G / 49G+ / 50G / 50G NewRPL / HP Prime / Oric 1 / Amstrad CPC 6128+

Avatar de l’utilisateur
zpalm
Fonctionne à 9600 bauds
Fonctionne à 9600 bauds
Messages : 2498
Inscription : 03 mai 2008 15:33
Localisation : Grenoble

Re: La (presque) question du mercredi

Message par zpalm » 09 mai 2019 09:31

C.Ret a écrit :
08 mai 2019 21:19
La vraie question, est comment avoir aussi vite trouvé la conjecture 6.t1-t0-2 ???
Est-ce un cas d'école des suites ?
Je suis aussi curieux de savoir comment dprtl a trouvé cette superbe formule.

Sinon pour toutes les suites la solution de facilité pour trouver les formules est de mettre les premiers termes dans le moteur de recherche de l'OEIS (On-Line Encyclopedia of Integer Sequences), mais pour les challenges sur le forum c'est un peu de la gruge.

Répondre

Revenir vers « Tous les Pockets »