

A271621


a(1) = 2, a(2) = 3, a(n) = A020639(a(n2)) + A006530(a(n1)).


1



2, 3, 5, 8, 7, 9, 10, 8, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,1


COMMENTS

Any sequence generated by this formula and any values for a(1) and a(2) will have a finite number of terms not equal to 4; i.e., all such sequences will eventually end up at 4 (and all following terms will be 4; 4 is the only term that can appear more than twice in a row in a sequence because it is the only number equal to the sum of its least and greatest prime factors). Example: a(1) = 77713; a(2) = 16; the sequence is: 77713, 16, 77715, 159, 56, 10, 7, 9, 10, 8, 4, 4, 4, ...


LINKS



FORMULA

a(n) = lpf(a(n2)) + gpf(a(n1)), where lpf(n) is the least prime dividing n and gpf(n) is the greatest prime dividing n.


EXAMPLE

a(1) = 13; a(2) = 46.
lpf(13) = 13; gpf(46) = 23.
a(3) = 13 + 23 = 36.


MATHEMATICA

a[1] = 2; a[2] = 3; a[n_] := a[n] = FactorInteger[a[n  2]][[1, 1]] +


PROG

(PARI) spf(n) = if (n==1, 1, vecmin(factor(n)[, 1]));
gpf(n) = if (n==1, 1, vecmax(factor(n)[, 1]));
lista(nn) = {print1(x=2, ", "); print1(y=3, ", "); for (n=1, nn, ny = spf(x) + gpf(y); print1(ny, ", "); x = y; y = ny; ); } \\ Michel Marcus, Apr 15 2016


CROSSREFS



KEYWORD

nonn


AUTHOR



STATUS

approved



