|
FORMULA
|
b(n) = b(n-1)+A005185[[n]] a(m)=if Product[b[n], {n, 1, m}]+1 is prime then Product[b[n], {n, 1, m}]+1
|
|
MATHEMATICA
|
Hofstadter[0] = Hofstadter[1] = 1 Hofstadter[n_Integer?Positive] := Hofstadter[n] = Hofstadter[n - Hofstadter[n - 1]] + Hofstadter[n - Hofstadter[n - 2]] a[1] = 2; a[n_] := a[n] = a[n - 1] + 2*Hofstadter[n] b = Flatten[Table[If[PrimeQ[Product[a[n], {n, 1, m}] + 1] == True, Product[a[n], {n, 1, m}] + 1, {}], {m, 1, 200}]]
|