|
|
A247838
|
|
Numbers n such that sigma(sigma(n)) is prime.
|
|
4
|
|
|
3, 2667, 3937, 57337, 172011, 253921, 677207307, 1073602561, 732959441001382539, 750688035198863979, 1000923107604038521, 1108158528150703969
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Numbers n such that A051027(n) is a prime p.
Prime 3 is the only prime p such that sigma(sigma(p)) is a prime q.
Conjecture: Subsequence of A046528 (numbers that are a product of distinct Mersenne primes).
Corresponding values of primes p: 7, 8191, 8191, 131071, 524287, 524287, ... (A247822). Conjecture: values of primes p is equal to Mersenne primes (A000668).
732959441001382539, 750688035198863979, 1000923107604038521, 1108158528150703969 and 196751176038481899983340171 are terms. - Jaroslav Krizek, Mar 25 2015
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
2667 is a term because sigma(sigma(2667)) = sigma(4096) = 8191 (i.e., prime).
|
|
MAPLE
|
|
|
MATHEMATICA
|
Select[Range[260000], PrimeQ[DivisorSigma[1, DivisorSigma[1, #]]]&] (* The program generates the first six terms of the sequence. *) (* Harvey P. Dale, Jan 18 2024 *)
|
|
PROG
|
(Magma) [n: n in [1..10000000] | IsPrime(SumOfDivisors(SumOfDivisors(n)))]
(PARI) isok(n) = isprime(sigma(sigma(n))); \\ Michel Marcus, Oct 01 2014
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|