OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..1000
MAPLE
with(numtheory): A066691:=n->`if`(tau(2*n-1) = tau(2*n+1), tau(2*n-1), NULL): seq(A066691(n), n=1..1000); # Wesley Ivan Hurt, Jan 03 2016
PROG
(PARI) for (n=2, 10000, if (numdiv(2*n-1) == numdiv(2*n+1), write1("tau=tau.txt", numdiv(2*n-1), ", ")))
(PARI) { n=0; for (m=1, 10^10, if ((d=numdiv(2*m-1)) == numdiv(2*m+1), write("b066691.txt", n++, " ", d); if (n==1000, return)) ) } \\ Harry J. Smith, Mar 17 2010
(Magma) [NumberOfDivisors(2*n-1): n in [1..400] | NumberOfDivisors(2*n-1) eq NumberOfDivisors(2*n+1)]; // Vincenzo Librandi, Jan 04 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Jon Perry, Jan 11 2002
EXTENSIONS
Offset changed from 0 to 1 by Harry J. Smith, Mar 17 2010
STATUS
approved