OFFSET
0,3
COMMENTS
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..10000
EXAMPLE
A={1,1,2,3,5,9,10,20,22,40,51,...} since A is the sorted union of:
A^2={1,2,5,10,20,40,67,126,203,354,571,908,1486,2250,3586,...} and
A^3={1,3,9,22,51,114,230,468,885,1674,3045,5418,9560,16341,...}.
PROG
(PARI) {a(n)=local(A=[1, 1]); for(i=1, #binary(3*n+1), A=vecsort(concat(Vec(Ser(A)^2), Vec(Ser(A)^3)))); A[n+1]}
for(n=0, 60, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Dec 09 2003
STATUS
approved