OFFSET
0,4
LINKS
PROG
(PARI)
permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
C(d, r)={sum(i=1, #r, my(t=r[i]); if(d%t==0, t))}
E(v) = {prod(i=2, #v, prod(j=1, i-1, my(g=gcd(v[i], v[j])); C(v[i]*v[j]/g, v)^g)) * prod(i=1, #v, my(t=v[i]); C(t, v)^((t-1)\2)*(C(t, v)-1)*if(t%2, 1, C(t/2, v)))}
a(n) = {my(s=0); forpart(p=n, s += permcount(p)*E(p)); s/n!} \\ Andrew Howroyd, Dec 10 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Christian G. Bower, Feb 15 1998 and May 15 1998
STATUS
approved
