login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

a(1)=2, a(2)=3, for n >= 3, a(n) = 2*(gpd(a(n-1)) + gpd(a(n-2))) + 1, where gpd(n) is the greatest prime divisor of n.
2

%I #26 Aug 20 2017 23:18:53

%S 2,3,11,29,81,65,33,49,37,89,253,225,57,49,53,121,129,109,305,341,185,

%T 137,349,973,977,2233,2013,181,485,557,1309,1149,801,945,193,401,1189,

%U 885,201,253,181,409,1181,3181,8725,7061,1313,817,289

%N a(1)=2, a(2)=3, for n >= 3, a(n) = 2*(gpd(a(n-1)) + gpd(a(n-2))) + 1, where gpd(n) is the greatest prime divisor of n.

%C The question about the boundedness of the sequence is equivalent to the question about its eventually periodicity. For example, the sequence defined by the same recurrence relation and initial values b(1)=2 and b(2)=5 is 2, 5, 15, 21, 25, 25, 21, 25, 25, ... so is periodic for n>=4 with the period {21,25,25}.

%C Problem. (a) Do there exist initial terms a(1) and a(2) depending on a given N for which the sequence has the least period of length >= N? (b) Do there exist initial terms a(1) and a(2) for which the sequence has no period?

%C Conjecture. Problem (a) is answered in affirmative, while problem (b) is answered in the negative.

%C This sequence is eventually periodic and therefore bounded: a(61)=a(85)=85 [sic] and a(62)=a(86)=73. [_D. S. McNeil_, Dec 14 2011]

%t a[1] := 2; a[2] := 3; a[n_] := a[n] = 2(FactorInteger[a[n - 1]][[-1, 1]] + FactorInteger[a[n - 2]][[-1, 1]]) + 1; Table[a[n], {n, 50}] (* _Alonso del Arte_, Dec 14 2011 *)

%t nxt[{a_,b_}]:={b,2(FactorInteger[a][[-1,1]]+FactorInteger[b] [[-1,1]])+ 1}; Transpose[NestList[nxt,{2,3},120]][[1]] (* _Harvey P. Dale_, Dec 15 2011 *)

%Y Cf. A006530

%K nonn

%O 1,1

%A _Vladimir Shevelev_, Dec 14 2011