login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A283877 Number of non-isomorphic set-systems of weight n. 206
1, 1, 2, 4, 9, 18, 44, 98, 244, 605, 1595, 4273, 12048, 34790, 104480, 322954, 1031556, 3389413, 11464454, 39820812, 141962355, 518663683, 1940341269, 7424565391, 29033121685, 115921101414, 472219204088, 1961177127371, 8298334192288, 35751364047676, 156736154469354 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
COMMENTS
A set-system is a finite set of finite nonempty sets. The weight of a set-system is the sum of cardinalities of its elements.
LINKS
FORMULA
Euler transform of A300913.
EXAMPLE
Non-isomorphic representatives of the a(4)=9 set-systems are:
((1234)),
((1)(234)), ((3)(123)), ((12)(34)), ((13)(23)),
((1)(2)(12)), ((1)(2)(34)), ((1)(3)(23)),
((1)(2)(3)(4)).
PROG
(PARI)
WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, (-1)^(n-1)/n))))-1, -#v)}
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}
K(q, t, k)={WeighT(Vec(sum(j=1, #q, my(g=gcd(t, q[j])); g*x^(q[j]/g)) + O(x*x^k), -k))}
a(n)={if(n==0, 1, my(s=0); forpart(q=n, my(g=sum(t=1, n, subst(x*Ser(K(q, t, n\t)/t), x, x^t) )); s+=permcount(q)*polcoef(exp(g - subst(g, x, x^2)), n)); s/n!)} \\ Andrew Howroyd, Jan 16 2024
CROSSREFS
Sequence in context: A000678 A362037 A368422 * A331850 A081490 A292478
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 17 2017
EXTENSIONS
a(0) = 1 prepended and terms a(11) and beyond from Andrew Howroyd, Sep 01 2019
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 23 02:53 EDT 2024. Contains 371906 sequences. (Running on oeis4.)