login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A274962 Numbers n such that sigma(n) and sigma(n) + 2 are both primes. 3
2, 2401, 19356878641, 46904541018721, 119601542190001, 360371335935601, 16472757578830081, 26835157974988801, 59879777952495601, 147669280778756881, 170589096345900241, 219660193449998401, 1103765757989399761, 1515946818108402241, 2044393722679974961, 2608728003079029841, 2805689752523610241, 3071293995460971361, 4537323492222149281, 9583348094642219041, 9982134924573725761 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Intersection of A249763 and A023194.
The next term, if it exists, must be greater than 5*10^12.
Each term > 2 is a square.
From Chai Wah Wu, Jul 13 2016: (Start)
Every term > 2 is of the form p^(2m) with p prime and m > 1. Proof: from the discussion in A023194, a term is of the form p^(2m). An odd term cannot be of the form n = p^2. If p = 6k+1, then sigma(n) = 36k^2 + 18k + 3 is composite. If p = 6k-1, then sigma(n) + 2 = 36k^2 - 6k + 3 is composite. Finally, 4 is not a term.
This could be the reason this sequence is so much sparser than A274963.
(End)
Terms cannot be of the form 2^(2m) since sigma(2^(2m)) + 2 = 2^(2m+1) + 1 is divisible by 3. - Altug Alkan, Jul 14 2016
Terms cannot be of the form 3^(2m) since sigma(3^(3m)) + 2 = 3(3^(2m) + 1)/2 is divisible by 3, i.e., all terms are of the form (6*k+1)^(2m) or (6*k-1)^(2m) - Chai Wah Wu, Aug 06 2016
Terms cannot be of the form p^6 since if p = 6*k+1, then sigma((6*k+1)^6) + 2 = 9*(5184*k^6 + 6048*k^5 + 3024*k^4 + 840*k^3 + 140*k^2 + 14*k + 1) and if p = 6*k-1 then sigma((6*k-1)^6) + 2 = 3*(15552*k^6 - 12960*k^5 + 4752*k^4 - 936*k^3 + 108*k^2 - 6*k + 1). Also note that terms cannot be of the form p^8 since if p = 6*k-1 then sigma((6*k-1)^8) = (1 - 6*k + 36*k^2)*(1 - 18*k + 432*k^2 - 4104*k^3 + 19440*k^4 - 46656*k^5 + 46656*k^6) and if p = 6*k+1 then sigma((6*k+1)^8) = 9*(186624*k^8 + 279936*k^7 + 186624*k^6 + 72576*k^5 + 18144*k^4 + 3024*k^3 + 336*k^2 + 24*k + 1). The least term that is of the form p^10 is 2089^10. So this partially explains why numbers of the form p^4 appear in this sequence most of the time in limited range. - Altug Alkan, Jul 15 2016
From Chai Wah Wu, Jul 20 2016: (Start)
If p^m > 2 is a term, then m == 4 mod 6 and p == 1 mod 6. Proof: Let q(k) be sigma(p^m) expressed as a polynomial in k. If p = 6k-1, then q(k) = 1 + (6k-1) + (6k-1)^2 + ... + (6k-1)^m.
The constant term of q(k) is 1-1+1-1+...-1+1 = 1 whereas the other coefficients are multiples of 6, i.e., q(k) = 1 + 6k*(...), thus sigma(p^m) + 2 is a multiple of 3.
Suppose p = 6k+1, then q(k) = 1 + (6k+1) + (6k+1)^2 + ... + (6k+1)^m. The constant term is m+1 and the other coefficients are multiples of 6, i.e., q(k) = (m+1) + 6k*(...).
This means that if m = 6r+2, then sigma(p^m) is a multiple of 3 and if m = 6r, then sigma(p^m) + 2 is a multiple of 3. End of Proof.
The following table lists the minimal k for r <= 4.
r | smallest k such that (6k+1)^(6r+4) is a term (A275237)
------------------------------------------------------------
0 | 1
1 | 348
2 | 436
3 | 6018
4 | 5880
For every prime p = 6k+1, does there exist r >= 0 such that(6k+1)^(6r+4) is a term?
(End)
Altug Alkan found that sigma((6k+1)^34) (i.e., the r = 5 case) is always composite (see comment in A275237). - Chai Wah Wu, Jul 21 2016
LINKS
EXAMPLE
2401 is in the sequence because sigma(2401) = 2801 and sigma(2401) + 2 = 2803 are both primes.
PROG
(Magma) [n: n in[1..10^7] | IsPrime(SumOfDivisors(n)) and IsPrime(SumOfDivisors(n)+2)]
(Python)
from sympy import isprime, divisor_sigma
A274962_list = [2]+[n for n, s in ((d**2, divisor_sigma(d**2)) for d in range(1, 10**3)) if isprime(s) and isprime(s+2)] # Chai Wah Wu, Jul 13 2016
(PARI) isok(n) = isprime(s=sigma(n)) && isprime(s+2); \\ Michel Marcus, Jul 14 2016
CROSSREFS
Sequence in context: A280312 A276650 A135234 * A199948 A261382 A281692
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Jul 12 2016
EXTENSIONS
a(4)-a(21) from Chai Wah Wu, Jul 13 2016
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 09:59 EDT 2024. Contains 371268 sequences. (Running on oeis4.)