login
Kaprekar triples: q such that q = x + y + z and q^3 = x*10^2n + y*10^n + z, where z < 10^n and n is the number of digits in q. q is not a power of 10 (except q=1).
(Formerly M4478)
9

%I M4478 #53 Oct 29 2023 21:39:49

%S 1,8,45,297,2322,2728,4445,4544,4949,5049,5455,5554,7172,27100,44443,

%T 55556,60434,77778,143857,208494,226071,279720,313390,324675,329967,

%U 346060,368928,395604,422577,427868,461539,472823,478115,488214,494208

%N Kaprekar triples: q such that q = x + y + z and q^3 = x*10^2n + y*10^n + z, where z < 10^n and n is the number of digits in q. q is not a power of 10 (except q=1).

%C The initial term a(1) = 1 is somewhat conventional: it is the only term with x = y = 0 and q = z = 10^k, which is explicitly allowed only for k = 0 and forbidden for k > 0. In all other cases, 0 < x, y, z < q, and q^3 has the same number of digits as x*10^2n. - _M. F. Hasler_, Aug 24 2017

%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

%D David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, NY, 1986, p. 151.

%H Jack Brennen and Hans Havermann, <a href="/A006887/b006887.txt">Table of n, a(n) for n = 1..1000</a> (first 200 terms from Giovanni Resta)

%H Futility Closet's "Math Notes", <a href="https://www.futilitycloset.com/2017/08/21/math-notes-116/">Shows the cubes of a(9) to a(13)</a>

%H Hans Havermann, <a href="http://gladhoboexpress.blogspot.ca/2017/08/cube-wonders.html">Cube wonders</a>

%H Douglas E. Iannucci and Bertrum Foster, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL8/Iannucci/iannucci45.html">Kaprekar Triples</a>, Journal of Integer Sequences, Vol. 8 (2005), Article 05.4.8.

%H R. Munafo, <a href="http://www.mrob.com/pub/seq/kaprekar.html">Kaprekar Sequences</a>

%e 1 = 0 + 0 + 1 and 1^3 = (00)1 (cf. comment),

%e 8 = 5 + 1 + 2 and 8^3 = 512,

%e 45 = 9 + 11 + 25, and 45^3 = 91125. - _M. F. Hasler_, Aug 24 2017

%t ok[n_] := n==1 || Block[{k = 10^IntegerLength[n], m = n^3}, n == Mod[m, k] + Floor[ m/k^2] + Mod[Floor[m/k], k] && ! IntegerQ@ Log10@ n]; Select[ Range@ 500000, ok] (* _Giovanni Resta_, Aug 23 2017 *)

%o (PARI) m=1;for(n=1,6,for(q=m+(n>1),-1+m*=10,q==sumdigits(q^3,m)&&print1(q","))) \\ _M. F. Hasler_, Aug 24 2017

%Y Cf. A291461.

%K nonn,base

%O 1,2

%A _Robert Munafo_

%E Entry revised by Larry Reeves (larryr(AT)acm.org), Apr 25 2001 and Dec 08 2002