Poate ca e bine sa se vada explicit si divizorii cu n cifre ai lui
10...001 scris cu (n+1) cifre,
folosind cele de mai sus in mod esential...
Folosind
sage
for n in [ 1..40 ]:
k = 10^n+1
if k % 7 == 0 :
print "%s : 7 =...\n...= %s\n" % ( k, k/7 )
1001 : 7 =...
...= 143
1000000001 : 7 =...
...= 142857143
1000000000000001 : 7 =...
...= 142857142857143
1000000000000000000001 : 7 =...
...= 142857142857142857143
1000000000000000000000000001 : 7 =...
...= 142857142857142857142857143
1000000000000000000000000000000001 : 7 =...
...= 142857142857142857142857142857143
1000000000000000000000000000000000000001 : 7 =...
...= 142857142857142857142857142857142857143
(143 este mereu la sfarsit. Numarul 142857 este cel ce inmultit cu 7...)
Daca teorema mica a lui Fermat este cunoscuta, se intelege usor perioda 6 pentru puterile n pentru care avem divizibiltate cu 7 pentru 10^n + 1 .
Daca nu, acest fapt (si calculul primelor zecimale ale lui 1/7) reprezinta o pregatire buna pentru ziua in care ii va veni vremea teoremei...