|
%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
|
