login
Numbers k such that sigma(k) + tau(k) is a prime.
11

%I #50 Dec 06 2022 07:56:51

%S 1,2,8,128,162,512,32768,41472,101250,125000,1414562,3748322,5120000,

%T 6837602,8000000,13530402,24234722,35701250,66724352,75031250,

%U 78125000,86093442,91125000,171532242,177058562,226759808,233971712,617831552,664301250,686128968

%N Numbers k such that sigma(k) + tau(k) is a prime.

%C The terms involve powers of small primes. - _Jud McCranie_, Nov 29 2001

%C From _Kevin P. Thompson_, Jun 20 2022: (Start)

%C Theorem: Terms that are greater than one must be twice a square.

%C Proof: Since sigma(k) is odd if and only if k is a square or twice a square, and tau(k) is odd if and only if k is a square, then an odd sum only occurs when k is twice a square, in which case sigma(k) is odd and tau(k) is even. So, these are the only candidates for sigma(k) + tau(k) being prime.

%C Theorem: No terms are congruent to 4 or 6 (mod 10).

%C Proof: Since no square ends in 2, 3, 7, or 8, and each term > 1 is twice a square, no term ends in 4 or 6. (End)

%H Amiram Eldar, <a href="/A064205/b064205.txt">Table of n, a(n) for n = 1..5000</a> (first 34 terms from Harry J. Smith, terms 35..276 from Kevin P. Thompson)

%e 128 is a term since sigma(128) + tau(128) = 255 + 8 = 263, which is prime.

%t Do[ If[ PrimeQ[ DivisorSigma[1, n] + DivisorSigma[0, n]], Print[n]], {n, 1, 10^7}]

%o (PARI) { n=0; for (m=1, 10^9, if (isprime(sigma(m) + numdiv(m)), write("b064205.txt", n++, " ", m); if (n==100, break)) ) } \\ _Harry J. Smith_, Sep 10 2009

%o (Python)

%o from itertools import count, islice

%o from sympy import isprime, divisor_sigma as s, divisor_count as t

%o def agen(): # generator of terms

%o yield 1

%o yield from (k for k in (2*i*i for i in count(1)) if isprime(s(k)+t(k)))

%o print(list(islice(agen(), 30))) # _Michael S. Branicky_, Jun 20 2022

%Y Cf. A007503 (sigma+tau), A065061, A055813.

%K nonn

%O 1,2

%A _Jason Earls_, Sep 21 2001

%E More terms from _Robert G. Wilson v_, Nov 12 2001

%E More terms from _Labos Elemer_, Nov 22 2001

%E More terms from _Jud McCranie_, Nov 29 2001

%E a(28) from _Harry J. Smith_, Sep 10 2009