|
%I #18 Dec 20 2022 11:22:31
%S 3,13,2593,2426113
%N Primes of the form 1 + Product_{j=1..k} b(j), where b(n) = b(n-1) + 2*A005185(n) and b(1) = 2.
%C Next term contains 886 digits. - _G. C. Greubel_, Dec 19 2022
%H G. C. Greubel, <a href="/A108583/b108583.txt">Table of n, a(n) for n = 0..4</a>
%t Hofstadter[n_]:= Hofstadter[n]= If[n<2, 1, Hofstadter[n-Hofstadter[n- 1]] +Hofstadter[n-Hofstadter[n-2]]]; (* A005185 *)
%t b[n_]:= b[n]= If[n==1, 2, b[n-1] +2*Hofstadter[n]]; (* A108585 *)
%t p[n_]:= p[n]= Product[b[j], {j,n}];
%t Select[Table[p[n] +1, {n, 500}], PrimeQ]
%Y Cf. A005185, A108585.
%K nonn,less
%O 0,1
%A _Roger L. Bagula_, Jul 05 2005
%E Edited by _G. C. Greubel_, Dec 19 2022
|
