OFFSET
1,1
PROG
(PARI) half(n) = if (n % 2, n\2, n/2 - 1);
a285524(n) = {v = vector(half(n), d, d!*(d + 1)!*(2^(n-2*d-1)*stirling(n-d, d+1, 2))); w = vecsort(v, , 1); w[#w]; }
lista(nn) = {v = vector(nn, n, a285524(n+3)); for (k=1, #v-4, if ((v[k] == v[k+1]) && (v[k] == v[k+2]) && (v[k] == v[k+3]), print1(v[k], ", ")); ); }
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Michel Marcus, Apr 20 2017
STATUS
approved