OFFSET
0,6
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..50
FORMULA
Euler transform of A342556.
PROG
(PARI) \\ See also A369932 for a more efficient program.
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}
edges(v, t) = {prod(i=2, #v, prod(j=1, i-1, my(g=gcd(v[i], v[j])); t(v[i]*v[j]/g)^g )) * prod(i=1, #v, my(c=v[i]); t(c)^((c-1)\2)*if(c%2, 1, t(c/2)))}
seq(n)={my(s=0); forpart(p=2*n, s+=permcount(p)*prod(i=1, #p, 1-x^p[i])*edges(p, w->1+x^w + O(x*x^n))); Vec(s/(2*n)!)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Jan 30 2024
STATUS
approved