login
Primes expressible as sigma(sigma(n)) - sigma(n), in order of their occurrence.
4

%I #20 Dec 31 2013 05:33:48

%S 3,7,17,31,31,41,23,73,127,73,89,127,463,523,241,523,157,241,523,463,

%T 211,257,131,983,379,1153,311,1153,83,983,521,4339,4339,113,8893,4339,

%U 4339,1093,4339,769,2851,8893,4339,1429,1097,4339,1093,4339,8893,4339,8893

%N Primes expressible as sigma(sigma(n)) - sigma(n), in order of their occurrence.

%H K. D. Bajpai, <a href="/A228567/b228567.txt">Table of n, a(n) for n = 1..2200</a>

%e a(9)= 127: sigma(sigma(93))-sigma(93)= 255-128= 127, which is prime.

%e a(11)= 89: sigma(sigma(98))-sigma(98)= 260-171= 89, which is prime.

%p with(numtheory):KD := proc() local a; a:= sigma(sigma(n))-sigma(n);if isprime(a) then RETURN (a); fi; end: seq(KD(),n=1..5000);

%Y Cf. A000203 (sigma(n): sum of divisors of n).

%Y Cf. A019279 (superperfect numbers: sigma(sigma(n))=2n).

%Y Cf. A033632 (numbers n: sigma(n)is prime).

%Y Cf. A051027 (sigma(sigma(n))).

%K nonn

%O 1,1

%A _K. D. Bajpai_, Nov 10 2013