Satz von Clement über Primzahlzwillinge
p(>2) und p+2 sind genau dann Primzahlen, wenn mod p(p+2) ist.
Das Problem Mod(n!,m):
Um die schnell wachsende Fakultät mit Mod in den Griff zu bekommen,
wird die Gesetzmäßigkeit Mod(a*b,m) = Mod(a,m)*Mod(b,m) genutzt,
wobei die Mod-Routine mit floor() formuliert wird, weil der Mod-Befehl in GGB
nicht mit Variablen klappt:
modfak(k,x) =
Iteration(modfak,{1,1},n-1) liefert dann Mod(n!,m).
Damit kann man die obergrenze auch auf 5000 hinauf schieben.
siehe https://de.wikipedia.org/wiki/Satz_von_Clement