A161598 - OEIS

%I #15 Nov 04 2013 13:13:09

%S 10,12,14,15,16,18,20,21,24,25,27,28,30,32,34,35,36,38,40,42,45,48,49,

%T 50,51,52,54,56,57,58,60,63,64,65,68,70,72,74,75,76,78,80,81,84,85,87,

%U 90,91,92,94,95,96,98,100,102,104,105,106,108,110,111,112,114,115,116

%N Numbers such that TITO(n) is not equal to n, where TITO(n) = A161594(n).

%C There are no prime numbers in the sequence: A010051(a(n)) = 0.

%H Reinhard Zumkeller, <a href="/A161598/b161598.txt">Table of n, a(n) for n = 1..10000</a>

%H T. Khovanova, <a href="http://blog.tanyakhovanova.com/?p=144">Turning Numbers Inside Out</a> [From _Tanya Khovanova_, Jul 07 2009]

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

%o (Haskell)

%o a161598 n = a161598_list !! (n-1)

%o a161598_list = filter (\x -> a161594 x /= x) [1..]

%o -- _Reinhard Zumkeller_, Oct 14 2011

%Y Complement of A161597.

%K nonn,base

%O 1,1

%A _J. H. Conway_ & _Tanya Khovanova_, Jun 14 2009

%E Edited by _N. J. A. Sloane_, Jun 23 2009

%E Offset corrected by _Reinhard Zumkeller_, Oct 14 2011