login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 


Cube root of A030678.
2

%I #12 May 18 2019 00:19:36

%S 5,3,7,16,8,4,9,20,21,48,61,15,164,82,188,92,96,215,466,481,23,236,

%T 242,533,253,258,122,576,587,277,606,615,624,632,64,301,656,143,312,

%U 679,1479,322,70,707,714,72,1565,34,1593,346,1619,1632,763,357,1668

%N Cube root of A030678.

%H Robert Israel, <a href="/A030679/b030679.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A030669(A002113(n+1)) = A030678(n)^(1/3). - _Robert Israel_, Dec 25 2018

%p N:= 3: # to get extensions of all palindromes of <= N digits

%p f:= proc(n) local d, x;

%p for d from 1 do

%p x:= ceil((n*10^d)^(1/3));

%p if x^3 < (n+1)*10^d then return x fi

%p od

%p end proc:

%p digrev:= proc(n) local i, L;

%p L:= convert(n, base, 10);

%p add(L[-i]*10^(i-1), i=1..nops(L))

%p end proc:

%p Res:= seq(f(i), i=1..9):

%p for d from 2 to N do

%p if d::even then

%p m:= d/2;

%p Res:= Res, seq(f(n*10^m + digrev(n)), n=10^(m-1)..10^m-1);

%p else

%p m:= (d-1)/2;

%p Res:= Res, seq(seq(f(n*10^(m+1)+y*10^m+digrev(n)), y=0..9), n=10^(m-1)..10^m-1)

%p fi

%p od: Res; # _Robert Israel_, Dec 25 2018

%Y Cf. A002113, A030669, A030678.

%K nonn,base

%O 1,1

%A _Patrick De Geest_

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 23 00:23 EDT 2024. Contains 376140 sequences. (Running on oeis4.)