login
n^a(n) is the smallest power of n (with a(n) > 1) which starts with n.
10

%I #28 Oct 08 2023 04:53:18

%S 2,8,18,6,24,10,20,11,22,2,26,14,10,8,92,6,166,5,19,11,60,39,48,51,94,

%T 42,467,86,14,66,58,3,5268,33,58,10,45,70,23,6,205,70,31,15,50,84,248,

%U 70,72,94,107,170,30,10139,182,136,87,318,49,10,178,237,295,99,17,206

%N n^a(n) is the smallest power of n (with a(n) > 1) which starts with n.

%C A070327(n) = n^a(n) for n > 1. - _Reinhard Zumkeller_, Sep 29 2011

%C The largest of the first 10000 terms is a(1955) = 119589011. - _Giovanni Resta_, May 27 2016

%H Giovanni Resta, <a href="/A051248/b051248.txt">Table of n, a(n) for n = 1..10000</a> (first 275 terms from Reinhard Zumkeller)

%e 2^8 = 256 begins with 2, and 2^k does not begin with 2 for any smaller integer k > 1, so a(2) = 8.

%t f[n_]:=Module[{i=2},While[Take[IntegerDigits[n^i],IntegerLength[n]] != IntegerDigits[n],i++];i]; Array[f,70] (* _Harvey P. Dale_, Oct 04 2011 *)

%o (Haskell)

%o import Data.List (isPrefixOf)

%o a051248 n = 2 + length

%o (takeWhile (not . (show n `isPrefixOf`) . show) $ iterate (* n) (n^2))

%o -- _Reinhard Zumkeller_, Sep 29 2011

%Y Cf. A045537, A104782.

%K nonn,base,nice

%O 1,1

%A _Michel ten Voorde_

%E More terms from _Jud McCranie_