login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A351651
a(n) is the quotient obtained when digsum(m^2) is divided by digsum(m), with digsum = sum of digits = A007953 and m = A351650(n).
1
1, 2, 3, 1, 1, 2, 3, 4, 1, 1, 2, 3, 4, 3, 2, 3, 4, 3, 2, 3, 1, 1, 2, 1, 3, 2, 2, 2, 1, 2, 1, 1, 2, 3, 4, 2, 2, 3, 4, 5, 3, 3, 4, 5, 3, 3, 2, 4, 3, 2, 2, 1, 2, 2, 2, 1, 1, 1, 2, 1, 1, 2, 3, 4, 3, 3, 2, 3, 4, 5, 3, 2, 4, 5, 2, 2, 3, 3, 3, 3, 2, 2, 2, 3, 2, 1, 3, 4, 3, 4, 5
OFFSET
1,2
COMMENTS
All positive integers are terms of this sequence (see A280012).
a(n) = 1 iff m = A351650(n) is a term of A058369 \ {0}.
a(n) = digsum(n) if m = A351650(n) is a term of A061909 \ {0}.
FORMULA
a(n) = A004159(A351650(n)) / A007953(A351650(n)).
EXAMPLE
A351650(8) = 13, then digsum(13) = 1+3 = 4 while digsum(13^2) = digsum(169) = 1+6+9 = 16; hence, a(8) = 16/4 = 4.
MATHEMATICA
Select[Total[IntegerDigits[#^2]]/Total[IntegerDigits[#]]& /@ Range[300], IntegerQ] (* Amiram Eldar, Feb 16 2022 *)
PROG
(PARI) lista(nn) = {my(list = List(), q); for (n=1, nn, if (denominator(q=sumdigits(n^2)/sumdigits(n))==1, listput(list, q)); ); Vec(list); } \\ Michel Marcus, Feb 16 2022
KEYWORD
nonn,base
AUTHOR
Bernard Schott, Feb 16 2022
EXTENSIONS
More terms from Michel Marcus, Feb 16 2022
STATUS
approved