|
%I #40 Apr 28 2016 21:08:11
%S 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,
%T 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,
%U 4,4,4,4,4,4,4
%N a(1) = 2, a(2) = 3, a(n) = A020639(a(n-2)) + A006530(a(n-1)).
%C 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, ...
%F a(n) = lpf(a(n-2)) + gpf(a(n-1)), where lpf(n) is the least prime dividing n and gpf(n) is the greatest prime dividing n.
%e a(1) = 13; a(2) = 46.
%e lpf(13) = 13; gpf(46) = 23.
%e a(3) = 13 + 23 = 36.
%t a[1] = 2; a[2] = 3; a[n_] := a[n] = FactorInteger[a[n - 2]][[1, 1]] +
%t FactorInteger[a[n - 1]][[-1, 1]]; Array[a, {120}] (* _Michael De Vlieger_, Apr 12 2016 *)
%o (PARI) spf(n) = if (n==1, 1, vecmin(factor(n)[,1]));
%o gpf(n) = if (n==1, 1, vecmax(factor(n)[,1]));
%o 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
%Y Cf. A020639, A006530, A074320.
%K nonn
%O 1,1
%A _Cody M. Haderlie_, Apr 10 2016
|
