OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..47
FORMULA
EXAMPLE
First few terms of A001917 (has offset 2) are 1, 1, 2, 1, 1, 2, 1, 2, 1, 6, so a(1) to a(3) are 2, 4, 11.
PROG
(Magma) W:=[]; r:=0; for n in [2..100000] do p:=NthPrime(n); a:=(p-1)/Modorder(2, p); if r lt a then r:=a; Append(~W, n); end if; end for; print W;
(Python)
from itertools import islice
from sympy import nextprime, n_order
def agen():
record, v, p = -1, 1, 3
while True:
if v > record: record = v; yield record
v, p = (p-1)//n_order(2, p), nextprime(p)
print(list(islice(agen(), 20))) # Michael S. Branicky, Oct 09 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Klaus Brockhaus, Dec 09 2008
EXTENSIONS
a(27)-a(37) from Amiram Eldar, Mar 08 2019
STATUS
approved