OFFSET
0,1
MATHEMATICA
f[n_] = 10 - Mod[Prime[n + 3], 10];
g[n_] = Mod[Mod[n, 5], 4];
h[n_] = g[f[n]];
digits = 15;
aa[1] = 2; aa[n_] := aa[n] = Mod[Mod[aa[n - 1]*(1 + h[n]), 5] + 1, 5];
Table[aa[n], {n, 1, digits^2/2}]
CROSSREFS
KEYWORD
nonn,base,uned,changed
AUTHOR
Roger L. Bagula, May 14 2005
STATUS
approved