OFFSET
1,1
LINKS
Ray Chandler, Table of n, a(n) for n = 1..2060 (first 501 terms from Max Alekseyev)
FORMULA
a(n)=m(p), where p=A066364(n) and m(p)=lcm(p, ord_p(10), m(q)) with q going over all prime divisors of ord_p(10).
EXAMPLE
a(6)=m(5477)=22494039 since it is the smallest m such that 10^m == 1 (mod m) and 5477|m.
PROG
(PARI) { m(p) = my(f, l, q); f=factorint(p)[, 1]; l=p; for(i=1, length(f), q=znorder(Mod(10, f[i])); l=lcm(l, q); l=lcm(l, m(q)) ); l }
CROSSREFS
KEYWORD
nonn
AUTHOR
Max Alekseyev, Nov 17 2005
STATUS
approved