OFFSET
1,7
COMMENTS
LINKS
T. D. Noe, Table of n, a(n) for n=1..1000
MATHEMATICA
nn=1000; a=b=1; t=Join[{1}, Table[c=((2n-1)b+3(n-1)a)/n; a=b; b=c; c, {n, 2, nn}]]; Table[p=Prime[i]; cnt=0; Do[If[Mod[t[[j]], p]==0, cnt++ ], {j, p}]; cnt, {i, PrimePi[nn]}]
CROSSREFS
KEYWORD
nonn
AUTHOR
T. D. Noe, Oct 24 2005
STATUS
approved