OFFSET
1,2
COMMENTS
736 terms < 10^7 were found.
All terms are squares because their number of divisors is odd (see formula field in A000005: a(n) is odd iff n is square).
EXAMPLE
1 is a term because tau(1) = 1, tau(1 + 1 - 1) = 1 and 1 modulo 2 is 1.
16 is a term because tau(16) = 5, tau(256 + 16 - 1) = 2 and 5 modulo 4 is 1.
50 is not a term because tau(50) = 6, tau(2500 + 50 -1) = 2 and 6 modulo 4 is 2.
PROG
(PARI) isok(k)=my(d_1=numdiv(k), d_2=numdiv(k^2+k-1)); d_1%(2*d_2)==1;
for(k=1, 1000, if(isok(k), print1(k", ")))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Claude H. R. Dequatre, Dec 07 2024
STATUS
approved