OFFSET
0,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..50 (terms 0..30 from Max Alekseyev)
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}
cross(u, v) = {sum(i=1, #u, sum(j=1, #v, gcd(u[i], v[j])))}
R(n, m, u)={if(n==0, 1, sum(k=if(m==1, n, 0), n\m, my(s=0); forpart(p=k, s+=self()(n-m*k, m-1, concat(u, Vec(p)))*2^cross(p, u)*permcount(p)); s/k!))}
a(n)={R(n, n, [])} \\ Andrew Howroyd, Sep 18 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Max Alekseyev following a suggestion from Franklin T. Adams-Watters, Sep 16 2020
STATUS
approved