OFFSET
1,3
COMMENTS
If n is a term, then 2*n is also a term, and vice versa. See A391584 for odd terms.
LINKS
PROG
(PARI) is_A391585(k) = (!k || (0==lift(Pol(binary(3*k)*Mod(1, 2)) % Pol(binary(k))*Mod(1, 2))));
(PARI)
up_to = 10001;
A391585list(up_to_n) = { my(v=vector(up_to_n), k=1, i=0); v[1] = 0; while(k<#v, i++; my(P3i=Pol(binary(3*i))*Mod(1, 2), P_i=Pol(binary(i))*Mod(1, 2)); if(0==lift(P3i % P_i), k++; v[k] = i)); (v); };
v391585 = A391585list(up_to);
A391585(n) = v391585[n];
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Antti Karttunen, Dec 18 2025
STATUS
approved
