A161597 Numbers such that TITO(n) = n, where TITO(n) = A161594(n).

1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 13, 17, 19, 22, 23, 26, 29, 31, 33, 37, 39, 41, 43, 44, 46, 47, 53, 55, 59, 61, 62, 66, 67, 69, 71, 73, 77, 79, 82, 83, 86, 88, 89, 93, 97, 99, 101, 103, 107, 109, 113, 121, 127, 131, 137, 139, 143, 149, 151, 157, 163, 167, 169, 173, 179

Offset: 1

Author

J. H. Conway & Tanya Khovanova, Jun 14 2009

Keywords

• nonn
• base

Comments

TITO(p) = p, for any prime p.

Links

• Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
• T. Khovanova, Turning Numbers Inside Out [From _Tanya Khovanova_, Jul 07 2009]

Crossrefs

Complement of A161598; nonprimes: A161600.

Programs

Haskell

a161597 n = a161597_list !! (n-1)
a161597_list = filter (\x -> a161594 x == x) [1..]
-- Reinhard Zumkeller, Oct 14 2011

Mathematica

reversepower[{n_, k_}] := FromDigits[Reverse[IntegerDigits[n]]]^k f[n_] := FromDigits[ Reverse[IntegerDigits[Times @@ Map[reversepower, FactorInteger[n]]]]] Select[Range[200], f[ # ] == # &]

PARI

is(n)={n==A161594(n)} \\ M. F. Hasler, May 11 2015

Extensions

Edited by N. J. A. Sloane, Jun 23 2009

Offset corrected by Reinhard Zumkeller, Oct 14 2011