A306361 Numbers k divisible by A101337(k) (narcissistic function).
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 50, 100, 110, 111, 153, 200, 221, 370, 371, 407, 500, 702, 1000, 1010, 1011, 1020, 1100, 1101, 1110, 1121, 1122, 1634, 2000, 2322, 4104, 5000, 8208, 9474, 10000, 10010, 10011, 10100, 10101, 10110, 11000, 11001, 11010, 11022, 11100, 11122, 11220, 12012, 12110, 12210, 12320, 14550
Offset: 1
Examples
For k = 20, 20 / (2^2 + 0^2) = 5; for k = 221, 221 / (2^3 + 2^3 + 1^3) = 13.
Programs
-
PARI
isok(n) = frac(n/A101337(n)) == 0; \\ Michel Marcus, Feb 11 2019
-
PARI
select( is_A306361=t->!(t%A101337(t)), [0..9999]) \\ M. F. Hasler, Nov 18 2019
Comments