%I #17 Aug 16 2019 12:31:48
%S 4,18,21,33,39,72,93,99,100,159,171,177,189,207,213,231,245,249,261,
%T 275,291,297,303,333,338,357,369,381,399,400,453,471,475,477,484,495,
%U 537,539,543,561,609,633,648,657,669,681,711,717,783,801,833,861,909,927
%N Numbers n such that sigma_2(n) - n^2 is prime.
%C Numbers n such that sum of the squares of the proper (or aliquot) divisors of n is a prime number.
%H Amiram Eldar, <a href="/A201880/b201880.txt">Table of n, a(n) for n = 1..10000</a>
%F {n: A067558(n) in A000040} - _R. J. Mathar_, Dec 07 2011
%e a(3)=21 because the aliquot divisors of 21 are 1, 3, 7, the sum of whose squares is 1^2 + 3^2 + 7^2 = 59, prime.
%p A067558 := proc(n)
%p numtheory[sigma][2](n)-n^2 ;
%p end proc:
%p isA201880 := proc(n)
%p isprime(A067558(n)) ;
%p end proc:
%p for n from 1 to 1000 do
%p if isA201880(n) then
%p printf("%d,",n);
%p end if;
%p end do: # _R. J. Mathar_, Dec 07 2011
%t Select[Range[400], PrimeQ[DivisorSigma[2, #]-#^2]&]
%o (PARI) is(n)=isprime(sigma(n,2)-n^2) \\ _Charles R Greathouse IV_, Dec 06 2011
%Y Cf. A001157, A037020, A067558.
%K nonn
%O 1,1
%A _Michel Lagneau_, Dec 06 2011
|