OFFSET
1,2
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..222
Martin Guerra and Doron Zeilberger, Maple program
MAPLE
# load programs from linked file:
seq(Pis([[0$5, 1], [1$6]], n, 400000), n=1..16);
PROG
(Python)
from math import lcm
from functools import lru_cache
from sympy import factorint
@lru_cache(maxsize=None)
def A381508(n):
if n == 1:
return 1
f = factorint(n).items()
if len(f) > 1:
return lcm(*(A381508(a**b) for a, b in f))
else:
k, x = 1, (0, 0, 0, 0, 1, 1)
while x != (0, 0, 0, 0, 0, 1):
k += 1
x = x[1:]+(sum(x) % n, )
return k # Chai Wah Wu, Apr 25 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Martin Guerra and Doron Zeilberger, Apr 24 2025
EXTENSIONS
a(17)-a(41) from Alois P. Heinz, Apr 25 2025
STATUS
approved
