A086838 Indices k where A057176(k) = 4.

3, 6, 9, 13, 35, 66, 115, 145, 297, 402, 491, 521, 1217, 2166, 3397, 3591, 4641, 4829, 6810, 6829, 7978, 8705, 17676, 21165, 50722, 52884, 84763, 121344, 121482, 162070, 358839, 560596, 1463827, 1593475, 1668394, 2086707, 2164032, 2517791, 4364421, 4970661, 19572673, 20471839, 31904004, 31960445, 51887007

Offset: 1

Author

Farideh Firoozbakht, Aug 08 2003

Keywords

• nonn

Comments

By definition of the sequence A087045(n)+2 and A086809(n)+1 for every n > 0 are in the sequence. There are no further terms up to 200000

No more terms < 3870000. - David Wasserman, Mar 31 2005

Crossrefs

Cf. A057176, A087045, A086809.

Programs

Mathematica

a[m_] := Sum[a[j], {j, 0, Mod[a[m-1], m]}]; vv={}; Do[ a[n]=If[n==0, 1, b]; v=a[n+1]; b=v; If[v==4, vv=Insert[vv, n+1, -1]; Print[vv]], {n, 0, 200000-1}]

Extensions

More terms from David Wasserman, Mar 31 2005

a(39)-a(45) from Michael S. Branicky, Mar 23 2024