OFFSET
0,3
LINKS
FORMULA
For a list n(1), n(2), n(3), ..., let fixF[n] = prod{i, j >= 1}(sum{d|[ i, j ]}(d*n(d))^((i, j)*n(i)*n(j)-(i=j)n(i))).
a(n) = Sum_{1*s_1+2*s_2+...=n} (fix A[s_1, s_2, ...]/(1^s_1*s_1!*2^s_2*s2!*...)) where fixA[s_1, s_2, ...] = Product_{i>=j>=1} f(i, j, s_i, s_j) where f(i, j, s_i, s_j) = {i=j} (Sum_{d|i} (d*s_d))^(i*s_i^2-s_i) or {i != j} (Sum_{d|lcm(i, j)} (d*s_d))^(2*gcd(i, j)*s_i*s_j).
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=1, #v, prod(j=1, #v, my(g=gcd(v[i], v[j])); C(v[i]*v[j]/g, v)^g)) / prod(i=1, #v, C(v[i], 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, May 15 1998 and Dec 03 2003
STATUS
approved
