|
|
A285087
|
|
Numbers n such that the number of partitions of n^2-1 is prime.
|
|
8
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Because asymptotically A000041(n^2-1) ~ exp(Pi*sqrt(2/3*(n^2-1))) / (4*sqrt(3)*(n^2-1)), the sum of the prime probabilities ~1/log(A000041(n^2-1)) is diverging and there are no obvious restrictions on primality; therefore, this sequence may be conjectured to be infinite.
a(5) > 50000.
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
13 is in the sequence because A000041(13^2-1) = 228204732751 is a prime.
|
|
PROG
|
(PARI) for(n=1, 2000, if(ispseudoprime(numbpart(n^2-1)), print1(n, ", ")))
(Python)
from itertools import count, islice
from sympy import isprime, npartitions
def A285087_gen(startvalue=1): # generator of terms >= startvalue
return filter(lambda n: isprime(npartitions(n**2-1)), count(max(startvalue, 1)))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,hard,more
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|