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”).
%I #27 May 04 2020 07:05:00
%S 1,1,1,1,1,1,1,1,1,4,2,7,4,3,4,7,5,4,6,7,8,9,7,55,22,11,55,22,55,22,
%T 55,55,33,55,55,55,55,44,55,55,55,55,505,88,66,202,262,77,121,88,99,
%U 181,151,121,101,88,505,424,181,121,151,181,131,343,202,181,141
%N a(n) is the palindrome equal to A334416(n) divided by its sum of digits.
%H <a href="https://oeis.org/wiki/Index_to_OEIS:_Section_Pac#palindromes">Index entries for sequences related to palindromes</a>
%F a(n) = A334416(n) / A007953(A334416(n)).
%e A334416(10) = 12 whose sum of digits is 3; 12/3 = 4, so a(10) = 4.
%t Select[#/Plus @@ IntegerDigits[#] & /@ Range[3000], PalindromeQ] (* _Amiram Eldar_, Apr 29 2020 *)
%o (PARI) isok(m) = iferr(my(d=digits(m/sumdigits(m))); d==Vecrev(d), E, 0);
%o apply(x->x/sumdigits(x), select(x->isok(x), [1..3000])) \\ _Michel Marcus_, Apr 29 2020
%Y Cf. A007953, A113315, A334416.
%Y Cf. A334375 (similar for primes).
%K nonn,base
%O 1,10
%A _Bernard Schott_, Apr 29 2020