|
EXAMPLE
|
Pairs (n, k) such that Eulerian(n, k) is prime are (4, 2), (10, 2), (8, 4), (14, 2), (16, 2), (12, 3), (15, 3), (26, 2), (50, 2), (56, 2), (70, 2), (51, 3), (27, 9), (72, 3), (116, 2), (87, 3), (183, 3).
|
|
PROG
|
(PARI) Eulerian(n, k)=sum(j=0, k, (-1)^j*(k-j)^n*binomial(n+1, j));
lista(nn) = {my(list=List()); for (n=1, nn, for (k=1, n, if (ispseudoprime(p=Eulerian(n, k)), listput(list, p)); ); ); Vec(Set(list)); } \\ Michel Marcus, May 25 2022
|