|
%I #12 Jul 10 2015 00:36:20
%S 7,13,43,31,211,157,307,757,307,1483,1723,3541,1723,5113,3907,8191,
%T 3541,8011,1723,6007,5113,5113,14281,5113,14281,8011,3541,28057,20593,
%U 20593,5113,37831,28057,17293,14281,8011,12433,28057,20593,14281,24181,10303,46441
%N Primes in A259184.
%C These primes are neither sorted nor uniqued. They are listed in the order found in A259184.
%H Robert Price, <a href="/A259186/b259186.txt">Table of n, a(n) for n = 1..1737</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) = A259184(A259185(n)).
%p with(numtheory): A259186:=n->`if`(isprime(1 - sigma(n) + sigma(n)^2), 1 - sigma(n) + sigma(n)^2, NULL): seq(A259186(n), n=1..200); # _Wesley Ivan Hurt_, Jul 09 2015
%t Select[Table[1 - DivisorSigma[1, n] + DivisorSigma[1, n]^2, {n, 10000}], PrimeQ]
%t Select[Table[Cyclotomic[6, DivisorSigma[1, n]], {n, 10000}], PrimeQ]
%Y Cf. A000203, A259184, A259185.
%K easy,nonn
%O 1,1
%A _Robert Price_, Jun 20 2015
|
