OFFSET
1,1
LINKS
Jon E. Schoenfield, Table of n, a(n) for n = 1..50
Jon E. Schoenfield, Magma code
Robert Tanniru, PARI Code
EXAMPLE
Fifth root of 33459 = 8.033459...
MATHEMATICA
f[n_] := Block[{l = Floor[ Log[10, n] + 1], rd = RealDigits[n^(1/5), 10, 24], id = IntegerDigits[n]}, rdd = Drop[ rd[[1]], rd[[2]]]; While[ rdd[[1]] == 0, rdd = Drop[ rdd, 1]]; Take[ rdd, l] == id]; Do[ If[ StringPosition[ ToString[ N[ n^(1/5), 24]], ToString[ n]] != {}, If[ f[n], Print[ n]]], {n, 2, 170000000}] (* Robert G. Wilson v, Jul 30 2004 *)
PROG
(PARI)
/* Uses PARI functions provided in link
* Note: does not predict 639 due to simplification error and
* 877-884 due to checking only first solutions to the Grafting Equation.
* Sample run uses a = [0, 16], b=10, p=5, direct=True */
GetAllGIs(0, 16, 10, 5, 1) \\ Robert Tanniru, Nov 20 2013
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Paul Lusch, Sep 06 2002
EXTENSIONS
Edited and extended by Robert G. Wilson v, Jul 31 2004
a(22)-a(25) by Robert Tanniru, Nov 20 2013
More terms from Jon E. Schoenfield, Aug 17 2014
STATUS
approved