A087045 Indices k where A057176(k) = 1.

1, 4, 64, 400, 489, 519, 2164, 3589, 8703, 84761, 358837, 1463825, 1668392, 20471837, 31960443

Offset

1,2

Comments

By definition k is in the sequence if and only if k divides A057176(k-1). - Farideh Firoozbakht, Aug 04 2003

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

Links

Table of n, a(n) for n=1..15.

Example

a(0) = 0 since A057176(0) = 1.

Maple

N := 10^5: A := array(0..N): A[0] := 1; for m from 1 to N do A[m] := add(A[j], j=0..modp(A[m-1], m)): od: a := NULL: for i from 0 to N do if A[i]=1 then a := a, i; fi; od; a;

Crossrefs

Cf. A057176.

Sequence in context: A165516 A064935 A030098 * A169801 A103751 A053959

Adjacent sequences: A087042 A087043 A087044 * A087046 A087047 A087048

Keyword

hard,more,nonn

Author

W. Edwin Clark and Farideh Firoozbakht at the suggestion of Leroy Quet, Aug 03 2003

Extensions

More terms from David Wasserman, Mar 31 2005

a(14)-a(15) from Michael S. Branicky, Mar 22 2024

Status

approved