OFFSET
1,1
COMMENTS
The sequence simulated: cc = Table[N[(Pi/2)*(b[n] + 7)], {n, 1, 3*digits}] behaves a lot like the rational approximation of the zeta zeros.
FORMULA
f[n] from either {2, 1, 0} or {0, 1, 2} a(n) = a(n-1) + 1+f[n]
Empirical g.f.: x*(x^9+2*x^8+x^7+x^6+2*x^5+3*x^4+3*x^3+2*x^2+x+2) / ((x-1)^2*(x^2+x+1)*(x^6+x^3+1)). [Colin Barker, Dec 02 2012]
MATHEMATICA
digits = 67 f = Flatten[Table[If[Mod[n, 2] == 0, {2, 1, 0}, {0, 1, 2}], {n, 1, digits}]] b[1] = 2; b[n_] := b[n] = b[n - 1] + 1 + f[[n]] aa = Table[b[n], {n, 1, 3*digits}]
CROSSREFS
KEYWORD
nonn,uned
AUTHOR
Roger L. Bagula, May 11 2005
STATUS
approved