login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Indices k where A057176(k) = 1.
2

%I #23 Mar 22 2024 19:52:14

%S 1,4,64,400,489,519,2164,3589,8703,84761,358837,1463825,1668392,

%T 20471837,31960443

%N Indices k where A057176(k) = 1.

%C By definition k is in the sequence if and only if k divides A057176(k-1). - _Farideh Firoozbakht_, Aug 04 2003

%C No more terms < 3870000. - _David Wasserman_, Mar 31 2005

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

%p 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;

%Y Cf. A057176.

%K hard,more,nonn

%O 1,2

%A _W. Edwin Clark_ and _Farideh Firoozbakht_ at the suggestion of _Leroy Quet_, Aug 03 2003

%E More terms from _David Wasserman_, Mar 31 2005

%E a(14)-a(15) from _Michael S. Branicky_, Mar 22 2024