A066691 - OEIS

%I #16 Sep 08 2022 08:45:04

%S 2,2,2,2,2,4,2,4,2,2,4,4,4,2,2,4,2,4,4,2,4,2,4,4,2,2,4,4,4,4,4,4,2,4,

%T 2,6,4,4,2,2,4,4,4,2,4,4,4,4,4,2,4,4,4,4,4,2,6,2,4,4,2,4,4,6,4,4,2,4,

%U 4,4,4,2,4,4,4,2,6,2,4,6,2,2,8,4,4,4,4,4,4,4,4,4,4,4,4,2,4,4,2,2,6,4,2,4,2

%N Value of tau(2n-1) when tau(2n-1) = tau(2n+1).

%H Harry J. Smith, <a href="/A066691/b066691.txt">Table of n, a(n) for n = 1..1000</a>

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

%o (PARI) for (n=2, 10000, if (numdiv(2*n-1) == numdiv(2*n+1), write1("tau=tau.txt", numdiv(2*n-1), ", ")))

%o (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

%o (Magma) [NumberOfDivisors(2*n-1): n in [1..400] | NumberOfDivisors(2*n-1) eq NumberOfDivisors(2*n+1)]; // _Vincenzo Librandi_, Jan 04 2016

%Y Cf. A000005 (tau).

%K nonn

%O 1,1

%A _Jon Perry_, Jan 11 2002

%E Offset changed from 0 to 1 by _Harry J. Smith_, Mar 17 2010