%I #20 Sep 08 2022 08:46:13
%S 3,757,262657,64008001,46656216001,30841155073,729027001,46656216001,
%T 30841155073,225200075257,885843322057,46656216001,41407378175593,
%U 885843322057,25002115044733,1126163480473,85766130261001,191102989824001,85766130261001,41407378175593
%N Primes of the form 1 + sigma(n)^3 + sigma(n)^6.
%C These primes are neither sorted nor uniqued. They are listed in the order found in A259369.
%H Robert Price, <a href="/A259371/b259371.txt">Table of n, a(n) for n = 1..756</a>
%H OEIS Wiki, <a href="https://oeis.org/wiki/Cyclotomic Polynomials at x=n, n! and sigma(n)">Cyclotomic Polynomials at x=n, n! and sigma(n)</a>
%F a(n) = A259369(A259370(n)).
%p with(numtheory): b:=n->1+sigma(n)^3+sigma(n)^6: A259371:=n->`if`( isprime(b(n)), b(n), NULL): seq(A259371(n), n=1..200); # _Wesley Ivan Hurt_, Jun 29 2015
%t Select[Table[1 + DivisorSigma[1, n]^3 + DivisorSigma[1, n]^6, {n, 10000}], PrimeQ]
%t Select[Table[Cyclotomic[9, DivisorSigma[1, n]], {n, 10000}], PrimeQ]
%o (Magma) [a: n in [1..150] | IsPrime(a) where a is 1 + SumOfDivisors(n)^3 + SumOfDivisors(n)^6]; // _Vincenzo Librandi_, Jun 26 2015
%Y Cf. A000203, A259369, A259370.
%K easy,nonn
%O 1,1
%A _Robert Price_, Jun 25 2015
|