Schachcomputer.info Community

Schachcomputer.info Community (https://www.schachcomputer.info/forum/index.php)
-   Schach und künstliche Intelligenz, Knobeleien, Denkspiele / Chess and artificial intelligence (https://www.schachcomputer.info/forum/forumdisplay.php?f=54)
-   -   Mathe: Das 3n+1 / 5n+1 Problem (https://www.schachcomputer.info/forum/showthread.php?t=6905)

dreihirn 14.08.2023 14:44

AW: Das 3n+1 / 5n+1 Problem
 
Hallo Thomas,
danke für das Rechnen.

Deine Beobachtungen passen in das allgemeine Schema.
Und vielleicht findet sich für 5n+-1 wegen seiner Glattheit
sogar ein Beweis der Konvergenz gegen 1 oder 7.

Zur Erinnerung für alle:
* Beim klassischen 3n+1-Problem gab es schon starke Ausreißer,
wo also die Folge lange umhersauste und zwischendurch auch
sehr grosse Werte annahm.

* Beim 3n+1/5n+1 waren die Ausreißer und grossen Werte
noch extremer.

3n+1 hat nach jedem *3-Schritt im Durchschnitt zwei Halbierungs-
schritte, pro Runde also im Schnitt eine Reduzierung um Faktor 3/4.
3n+1/5n+1 hat in jeder Runde *3*5 und dazu im Durchschnitt
zwei + zwei Halbierungsschritte; pro Runde also im Schnitt eine
Reduzierung um Faktor 15/16. Das ist größer als 3/4, aber noch
kleiner als 1.

5n+-1 hat nach jedem *5-Schritt mindestens zwei Halbierungsschritte,
im Durchschnitt aber 3 Halbierungsschritte, die sich ergeben aus
2 Schritte mit 1/2
3 Schritte mit 1/4
4 Schritte mit 1/8
5 Schritte mit 1/16
Die 3 ergibt sich als als gwichtete Summe aus den Einzelfällen.
5/8 < 3/4 < 15/16 < 1, also sollte 5n+-1 das glatteste
Reduzieren haben.
Andererseits würde der Worst Case: nach jedem *5 nur zwei
Halbierungsschritte ein beliebig grosses Ansteigen erlauben.

**************************************
Für mögliche Beweise scheint mir:
Ein Beweis für 3n+1 würde ziemlich sicher einen Beweis für 5n+-1
nach sich ziehen. Im Gegensatz könnte es einen Beweis für 5n+-1
geben, der nicht "fast-automatisch" einen 3n+1-Beweis nach sch zieht.

Viele Grüße, Ingo.

dreihirn 18.08.2023 09:11

AW: Das 3n+1 / 5n+1 Problem
 
Mir ist jetzt eine allgemeine Vermutung klar geworden.
Betrachtet wird ein System, was auf ungeraden Zahlen
k(1), k(2), ... k(s) basiert. Dabei ist s eine natürliche Zahl
ungleich 0.

Es wird bei einer ungeraden natürlichen Zahl n losgerechnet.
Jede Runde hat s Etappen.
1. Aus n wird k(1)*n +1. Dann wird runterhalbiert.
2. Sei m das Ergebnis der vorherigen Etappe. Bilde k(2)*m +1
und halbiere dann runter.
...
Etappe s: Sei r das Ergebnis der vorherigen Etappe. Bilde
k(s)*r + 1 und halbiere runter.

Das Ergebnis ist der Ausgangswert der nächsten Runde.



Sei B= k(1) * k(2) * ... * k(s).
Der entscheidende Wert ist jetzt C = B / (4^s).

Vermutung:
Ist C < 1, läuft jeder Startwert in einen von endlich vielen Zykeln.
Ist C > 1, läuft fast jeder Startwert (im Sinne von asymptotischer Dichte 1)
nach unendlich.
C=1 kann nicht vorkommen, da B ungerade sein muss.

Begründung: Runterhalbieren besteht im Durchschnitt aus zwei Schritten.
Gleichzeitig wird in einer Etappe mit k(i) multippliziert. Also liefert die
Etappe im Schnitt Faktor k(i) / 4.

Beispiele:
s=1 und k(1)=3 ist der klassische Collatz-Fall.
s=2 und k(1)=3, k(2)=5 steht oben im Thread, von Gilgamesch gerechnet.

Allgemein dürften Instanzen mit kleinem C-Wert schnelles Abstürzen
in die Zykel geben, und C-Werte nahe 1 ganz langsames. Ein vermutetes
Beispiel mit sehr hohen Zwischenwerten
s=3, k= (3,3,7). Dafür ist der C-Wert
3*3*7 / 64 = 63/64.

Viele Grüße, Ingo.

Gilgamesch 18.08.2023 23:04

AW: Das 3n+1 / 5n+1 Problem
 
Hallo Ingo,

die Vermutung ist wohl genau so richtig.

Ich hatte gewisse Schwierigkeiten bei der Berechnung von (3,3,7) und mußte auf einen unsigned Integer Wert mit 256 Bit wechseln.
Auf jeden Fall scheinen alle Werte in einen Zyklus mit einem kleinem Fixpunkt zu laufen.
Die höchsten Werte können aber gelegentlich sehr groß werden, bei dem Startwert 2799 ist selbst mit 256 Bit Schluss.

Einige Beispiele, Spalten:
- Startwert
- Runden zum Erreichen des Fixpunktes (stimmt nur ungefähr)
- Schritte zum Erreichen des Fixpunktes (stimmt nur ungefähr)
- Fixpunkt (kleinster vorkommen Wert der sich wiederholt)
- Höchster vorkommender Wert

Code:

    1      2    10        1                                                                        8
    3      2    14        1                                                                      16
    5      2    12        1                                                                      16
    7      8    62      15                                                                    7680
    9      4    25      39                                                                      624
    11      4    28      35                                                                      372
    13      2    16        1                                                                      40
    15      7    56      23                                                                    7680
    17      5    35      39                                                                      624
    19      4    26      39                                                                      624
    21      2    14        1                                                                      64
    23      3    20      35                                                                      372
    25      3    23      29                                                                      232
    27      3    23      31                                                                      328
    29      8    64      15                                                                    7680
    31    11    91      15                                                                    7680
...
  199    219  1968        1                                                          30143685317248
  201    41    363      15                                                              1257744720
  203      6    53      29                                                                    1604
  205      3    26      29                                                                      616
  207    43    381      15                                                              1257744720
  209    44    390      15                                                              1257744720
  211    223  2004        1                                                          30143685317248
...
  473    19    168      35                                                                  397524
  475      6    51      43                                                                    5620
  477    10    93        1                                                                  110132
  479  1507  13525      15                                5633474152358653058455030362573301198800
  481    11    95      15                                                                    7680
  483      8    68      15                                                                    7680
  485      8    68      15                                                                    7680
...
  1649      3    29      29                                                                    4948
  1651      7    65      29                                                                    7432
  1653      4    38      29                                                                    4960
  1655  5688  51062      35    532883555724939922698357411067488379361088412281800423855405840357376
  1657      4    38      29                                                                    13056
  1659    221  1989        1                                                      562413204037782208
  1661    276  2477      35                                                        34810817575894536
  1663    96    864      29                                                                322328464
...
  1961    21    188      35                                                                  397524
  1963      9    79      15                                                                    15464
  1965      9    79      15                                                                    7680
  1967  5159  46312      15    141331469837768875083859002482956122113450198439139246988547504397056
  1969      9    79      15                                                                    7680
  1971      9    79      15                                                                    8872
  1973      9    79      15                                                                    7680
  1975    110    987      35                                                            543379976224
  1977  1509  13545      15                                5633474152358653058455030362573301198800
  1979    12    113        1                                                                  346240
  1981    22    197      35                                                                  397524
  1983    22    197      35                                                                  397524
  1985    13    115      15                                                                    8800
...
  2787      4    33      43                                                                    12544
  2789    34    311        1                                                                29812096
  2791    34    311        1                                                                29812096
  2793    34    311        1                                                                29812096
  2795      4    33      43                                                                    22016
  2797    45    404      35                                                                  2569008
2799  Overflow

Viele Grüße,
Thomas

dreihirn 19.08.2023 10:40

AW: Das 3n+1 / 5n+1 Problem
 
Hallo Thomas,
danke für das Rechnen.

Zitieren:

Zitat von Gilgamesch (Beitrag 118968)
...
Die höchsten Werte können aber gelegentlich sehr groß werden, bei dem Startwert 2799 ist selbst mit 256 Bit Schluss.

...
1655 5688 51062 35 53288355572493992269835741106748837936108841228180 0423855405840357376
1657 4 38 29 13056
1659 221 1989 1
...

Das Zwischen-Maximum für 1655 ist ja echt beeindruckend gross.
Und daneben dann so bescheidene Verläufe für 1657 und 1659..

Dank und Gruß, Ingo.

Gilgamesch 19.08.2023 17:58

AW: Das 3n+1 / 5n+1 Problem
 
1 Anhang/Anhänge
Anbei ein Zip File mit der langen 1655 - (3,3,7) Folge.
Etwas Auffälliges ausser der Länge kann ich nicht erkennen.

dreihirn 20.08.2023 09:29

AW: Das 3n+1 / 5n+1 Problem
 
3 Anhang/Anhänge
Hallo Thomas,
danke für die Daten.

Zitieren:

Zitat von Gilgamesch (Beitrag 118998)
Etwas Auffälliges ausser der Länge
kann ich nicht erkennen.

Aber das reicht doch schon. Ich habe mal bei Schritgröße 2
drei Screenshots gemacht: vom Anfang der Folge, irgendwo
aus der Mitte, und am Schluss. Natürlich erkennt man da die
Ziffern nicht mehr, aber die Längen der Zahlen.

Man sieht die Ähnlichkeit zu einem Casino mit "positivem
Erwartungswert", wo in jeder Runde log(3) bzw log(7) zu
zahlen ist und dann zufällig ein Gewinn von log(2) oder
2*log(2) oder 3*log(2) oder ... ausgezahlt wird, so ähnlich wie
bei den Gewinnleitern der Gauselmann-Geldautomaten.
(Wobei natürlich die Gauselmann-Automaten aus Sicht der
Spieler einen negativen Erwartungswert haben.)

Die große Länge der Folge zeigt, dass man auch bei gutem
Erwartungswert (hier 63/64 < 1) manchmal lange Durststrecken
überwinden muss, eher man wirklich ins Plus kommt.

Dank und Gruß, Ingo.

Gilgamesch 20.08.2023 18:38

AW: Das 3n+1 / 5n+1 Problem
 
Jetzt werden die Zahlen wirklich groß.

Ich hatte schon länger nach einer C++ Library gesucht, die beliebig große natürliche Zahlen zulässt und endlich eine nicht zu komplexe Implementierung gefunden.
Damit habe ich die Overflows der 256 Bit Variante (3,3,7) nochmal überprüft.

Ergebnis:

Code:

Startwert:  2799
Runden:    14966
Schritte: 134261
Fixpunkt:      1
Höchster Wert:
4576380734615710479690106325689609211241863537762981964042455051075724332822215383625091530700741156
546681974982555045360453632 (127 Stellen)


Code:

Startwert: 14975
Runden:    51639
Schritte: 463293
Fixpunkt:    35
Höchster Wert:
9989323727783341288713747130739468253471810232357536641687347870748657993683879904466884139557940582
5553160183246255023839530657063922381301089297022146122314293811553434648714267582962923939302464618
9655985406114827467681702953665986109731506425351328900741453810959854609376  (276 Stellen)

Grüße,
Thomas

dreihirn 20.08.2023 21:52

AW: Das 3n+1 / 5n+1 Problem
 
Hallo Thomas,

abgefahren!
Wenn das der alte Collatz noch erlebt hätte.

Zitieren:

Zitat von Gilgamesch (Beitrag 119019)
Code:

Startwert:  2799
Runden:    14966
Schritte: 134261


Code:

Startwert: 14975
Runden:    51639
Schritte: 463293


Interessant ist der Quotient Schritte / Runden.
Bei 2799 ist er 8,971;
bei 14975 ist er 8,972.

8,97 circa= 9, und das passt gut.
In jeder Runde drei *-Schritte und durchschnittlich
sechs Halbierungen.

Mit Dank und Gruss, Ingo.

Wandersleben 20.08.2023 22:06

AW: Das 3n+1 / 5n+1 Problem
 
Zitieren:

Zitat von Gilgamesch (Beitrag 119019)
Code:

Startwert: 14975
Runden:    51639
Schritte: 463293
Fixpunkt:    35
Höchster Wert:
9989323727783341288713747130739468253471810232357536641687347870748657993683879904466884139557940582
5553160183246255023839530657063922381301089297022146122314293811553434648714267582962923939302464618
9655985406114827467681702953665986109731506425351328900741453810959854609376  (276 Stellen)


Vielen dank, Thomas!
Mathematik von ihrer schönsten seite!

Viele grüße
Horst

dreihirn 21.08.2023 04:43

AW: Das 3n+1 / 5n+1 Problem
 
Hallo Horst, hallo Thomas,

Zitieren:

Zitat von Wandersleben (Beitrag 119031)
Vielen dank, Thomas!
Mathematik von ihrer schönsten seite!

die Zahl ist schon monströs gross. Nur bei Primzahlen der
Form 2^m - 1 habe ich bisher größere Zahlen erlebt.


Vielleicht liefert folgende Variante des
Collatz-Problems noch grössere Zahlen für
kleine Startwerte:

n+1 / 3n+1 / 5n+1 / 17n+1
dazwischen jeweils runterhalbieren.
Also hat jede Runde vier Etappen.

Der Clou ist das Produkt 1*3*5*17 = 255.
255 / (4^4) = 255 / 256, also nur ganz knapp unter 1.

Spekulation: Es sollte 3- oder 4-stellige Startwerte
im Dezimalsystem geben, bei denen zwar Konvergenz eintritt,
aber zwischendurch mit 400-stelligen Werten.

Viele Grüße, Ingo.

Wandersleben 21.08.2023 13:41

AW: Das 3n+1 / 5n+1 Problem
 
Moin, Ingo

da werden wir Thomas eine eigene windstromtrasse von der Nordsee aus legen lassen müssen, denn nach deinem prinzip kann man die zyklen immer weiter ausbauen.
Es macht aber enorm spaß in den zahlenreihen zu stöbern, zumal ich mit Q-Basic nicht weit gekommen bin. :o

Viele grüße auch an Thomas
sendet Horst

dreihirn 21.08.2023 16:41

AW: Das 3n+1 / 5n+1 Problem
 
Hallo Horst,

Zitieren:

Zitat von Wandersleben (Beitrag 119046)
da werden wir Thomas eine eigene windstromtrasse
von der Nordsee aus legen lassen müssen, denn nach
deinem prinzip kann man die zyklen immer weiter ausbauen.

Ich habe einen anderen Kumpel in der Szene, wo Raumfahrt-
Trajektorien berechnet werden. Es ist Dr. Dietmar Wolz.

https://althofer.de/space-trajectories.html

Als er vor einigen Jahren mal wirklich sehr arg gerechnet hat,
habe ich Dietmar 707,- Euro Stromgeld überwiesen. Das war ein
Monatsverbraucht von ihm.


Zitieren:

zumal ich mit Q-Basic nicht weit gekommen bin. :o
Tipp: Mach die Arithmetik wie in der Grundschule,
für jede Ziffer eine Speicherstelle. Damit solltest Du
zumindest 3n+1 5n+1 für 95 hinkriegen.

Der Experte Jürgen Dankert hat es übrigens in seinem
Collatz-Applet genau so gemacht.
https://www.juergendankert.de/spezma.../collatzm.html
Da kann man Zahlen mit beliebig vielen Stellen eingeben.

Viele Grüße, Ingo.

Gilgamesch 21.08.2023 17:36

AW: Das 3n+1 / 5n+1 Problem
 
Hallo zusammen,

die (1,3,5,17) Sequenz verhält sich wirklich ungewöhnlich.
Über weite Strecken (z.B. 4000 Werte) geht die Sequenz recht schnell zu einem (niedrigen) Fixpunkt bzw. einer Wiederholung über, ohne dabei extrem hohe Wert zu erreichen.
Ab und zu "explodiert" die Sequenz aber regelrecht und erreicht extrem hohe Schrittzahlen und Werte.

Das Verhalten könnte vielleicht auch mit der sogenannten "Glide" Funktion bzw. mit den "Modular restrictions" zu tun haben.
Für viele Startwerte fällt eine Sequenz (egal welche, d.h. auch die normale Collatz Sequenz) in wenigen Schritten unter den Startwert.
z.B. für (Startwert mod 256) fallen viele Restklassen (immer) fast sofort unter den Startwert und sind deswegen nicht interessant.

Code:

Startwert:  4165
Runden:    180902
Schritte: 2169817
Fixpunkt:    748
Höchster Wert:

2161143116297822246586194619977508617335044059394982400266868709463203826758789684321154213886571484
2228598324912387854660523318189767877241929630107753167804425983346453349482982502636150662249264775
2246597978552709459354960708982804986752801679761998930824934057771597972364194568504305802666521191
9434019245847927522149242497296232485701911808641372732069755093846187580228232787975884 (388 Stellen)

Fixpunkt-Schleife:

180903  2169818  1  748
180903  2169819  2  374
180903  2169820  2  187
180903  2169821  3  562
180903  2169822  2  281
180903  2169823  5  1406
180903  2169824  2  703
180903  2169825  17 11952
180903  2169826  2  5976
180903  2169827  2  2988
180903  2169828  2  1494
180903  2169829  2  747


Code:

Startwert:  6021
Runden:    172580
Schritte: 2069999
Fixpunkt:    2260
Höchster Wert:

5002486244923240312001743017936892183504525045175114686760177504439694985925177642231418841638452238
4420897266822081684487027504394784367736643993207339613861899663069829502356060560319158033166297379
9556318974804775045635970211811150887820530931434906754222727758428335449847793295309568320562692662
0399648614783855598078787791050982586618350972418048252851818368843319939450277431441121279081936047
3347327942895811507588245166492526382484418480794524970004664633879837338331940388136662349394172332
33497312639867311674001436257115440512990507742 (547 Stellen)

Fixpunkt-Schleife:

172581  2070000  1  2260
172581  2070001  2  1130
172581  2070002  2  565
172581  2070003  3  1696
172581  2070004  2  848
172581  2070005  2  424
172581  2070006  2  212
172581  2070007  2  106
172581  2070008  2  53
172581  2070009  5  266
172581  2070010  2  133
172581  2070011  17 2262
172581  2070012  2  1131

172582  2070013  1  1132
172582  2070014  2  566
172582  2070015  2  283
172582  2070016  3  850
172582  2070017  2  425
172582  2070018  5  2126
172582  2070019  2  1063
172582  2070020  17 18072
172582  2070021  2  9036
172582  2070022  2  4518
172582  2070023  2  2259

Viele Grüße, Thomas

dreihirn 21.08.2023 18:54

AW: Das 3n+1 / 5n+1 Problem
 
Hallo Thomas,

Hammer! Einfach nur Hammer!

Aus Neugier frage ich mal: wieviel Rechenzeit
auf welcher Hardware hast Du für das 1-3-5-17-
Problem investiert?

Restlos begeistert, Ingo.

Gilgamesch 21.08.2023 20:18

AW: Das 3n+1 / 5n+1 Problem
 
Hallo Ingo,

so groß war der reine Rechenaufwand gar nicht mal, die Library scheint recht gut optimiert zu sein.
Ein Durchlauf einer solchen langen Sequenz hat etwas über 8 Minuten auf einem CPU Core gedauert.

Ich habe einen Haswell Core i7-4790K @ 4.00 GHz benutzt.
8 Minuten / 2.000.000 Steps => 960.000 Takzyklen / Schritt
Da die Zahlen mühelos in den first Level Cache der CPU (32 KB) passen, kommt man mit 960.000 Takzyklen pro Schritt offenbar hin.

Insgesamt habe ich dann ca. 12 Durchläufe gebraucht bis alles gepasst hat, also 96 Minuten Rechenzeit auf einem Core.
Da würde auf den schnelleren Rechnern und mit Multi Tasking noch mehr gehen.
Das Anpassen des Sourcecodes und Kontrollieren der Ausgabe ist da aufwendiger, geht aber alles noch in vertretbarer Zeit.


Viele Grüße, Thomas

dreihirn 24.08.2023 18:26

AW: Das 3n+1 / 5n+1 Problem
 
Hallo,

an dem 3n+1 Problem sitze ich ja seit 2011
immer wieder mal. Jetzt habe ich mich dazu
durchgerungen, ein paar Preise für theoretische
Beweise auszusetzen:

https://www.althofer.de/collatz-prizes.html

Die sind natürlich nicht nur für Leute aus diesem
Forum gedacht, aber eben auch.

Laufzeit der Preise bis Ende 2037.

Frohes Probieren,

Ingo.

dreihirn 10.09.2023 17:55

AW: Das 3n+1 / 5n+1 Problem
 
Hallo,

gerade habe ich im Internet gefunden, dass es ein jabanisches KI-
Startub-Unternehmen gibt, was im Juli 2021 einen 120-Millionen-
Yen-Breis für die Lösung des klassischen Collatz-Broblems bis
Juli 2031 ausgesetzt hat. Die 120 Mio Yen entsbrechen etwa
800.000 US-Dollar, wogegen mein 1.000-Euro Breis bis
Ende 2037 natürlich Beanuts ist.

Hier die Webseiten des Breises und der Firma:
https://mathprize.net/posts/collatz-conjecture/
https://bakuage.com/en/
https://bakuage.com/en/about/index.html

Als Kabital der Firma sind 3 Mio Yen angegeben, also circa
nur 20.000 US-Dollar. Ich kenne mich mit den jabanischen
Marktregularien nicht aus, aber wenn es eine GmbH mit nur
circa 20.000 US-Dollar Einlage wäre, dann sollte sich ein
Löser des Broblems nicht zu viel Hoffnungen machen, dadurch
knabber Dollar-Millionär zu werden.

Viele Grüße, Ingo.

Zum Vergleich die Regeln meines bescheidenen Breises:
https://www.althofer.de/collatz-prizes.html

dreihirn 12.01.2024 10:13

AW: Das 3n+1 / 5n+1 Problem
 
Hallo,

Ende November 2023 war die Tagung "Advances in Computer Games 2023",
bei der ich Ergebnisse von Michael Hartisch, Thomas Zipproth und mir
vorstellen durfte.

Die Tagung war onlline. Jetzt hängen alle Vorträge und Diskussionsteile
bei Youtube online, unter
https://www.youtube.com/playlist?lis...hOINjLHeuf9q7-

Die Präsentation zum 3n+-1-Spiel und zum 3n+1 / 5n+1 Problem
findet sich in den ersten 22 Minuten des Videos von Tag 3:
https://www.youtube.com/watch?v=-BWh...uf9q7-&index=9

Chairman der Sitzung war Bruno Bouzy, was die Gelegenheit gibt, direkt
nebeneinander französisches und westfälisches Englisch zu hören.
Keine Angst: Brunos und mein Englisch ist einfach, und der Grossteil
des Textes findet sich auch auf den gezeigten Folien.

Viele Grüsse, Ingo.


Alle Zeitangaben in WEZ +1. Es ist jetzt 08:52 Uhr.

Powered by vBulletin (Deutsch)
Copyright ©2000 - 2025, Jelsoft Enterprises Ltd.
©Schachcomputer.info