OFFSET
1,1
COMMENTS
Partial sums of A194232.
The Mathematica code computes 50 terms, but only these 45 terms are correct.
LINKS
MATHEMATICA
Inv[a_, b_] := PowerMod[a, -1, b]; t = {1}; Do[h = 10^k-1; d = Divisors[h]; d2 = Select[d, GCD[#, h/#] == 1 &]; AppendTo[t, h]; Do[AppendTo[t, d2[[i]]*Inv[d2[[i]], h/d2[[i]]]], {i, 2, Length[d2]-1}], {k, 50}]; t = Union[t]; Table[Position[t, 10^n-1, 1, 1][[1, 1]], {n, Log[10, t[[-1]]]}]
CROSSREFS
KEYWORD
nonn,base
AUTHOR
T. D. Noe, Aug 17 2011
STATUS
approved