cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 51-60 of 466 results. Next

A035053 Number of connected graphs on n unlabeled nodes where every block is a complete graph.

Original entry on oeis.org

1, 1, 1, 2, 4, 9, 22, 59, 165, 496, 1540, 4960, 16390, 55408, 190572, 665699, 2354932, 8424025, 30424768, 110823984, 406734060, 1502876903, 5586976572, 20884546416, 78460794158, 296124542120, 1122346648913, 4270387848473
Offset: 0

Views

Author

Christian G. Bower, Oct 15 1998

Keywords

Comments

Equivalently, this is the number of "hypertrees" on n unlabeled nodes, i.e., connected hypergraphs that have no cycles, assuming that each edge contains at least two vertices. - Don Knuth, Jan 26 2008. See A134955 for hyperforests.
Graphs where every block is a complete graph are also called block graphs or clique tree. They can be characterized as induced-diamond-free chordal graphs. - Falk Hüffner, Jul 25 2019

Examples

			From _Gus Wiseman_, May 20 2018: (Start)
Non-isomorphic representatives of the a(5) = 9 hypertrees are the following:
  {{1,2,3,4,5}}
  {{1,5},{2,3,4,5}}
  {{1,2,5},{3,4,5}}
  {{1,2},{2,5},{3,4,5}}
  {{1,4},{2,5},{3,4,5}}
  {{1,5},{2,5},{3,4,5}}
  {{1,3},{2,4},{3,5},{4,5}}
  {{1,4},{2,5},{3,5},{4,5}}
  {{1,5},{2,5},{3,5},{4,5}}
(End)
		

References

  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 71, (3.4.14).

Crossrefs

Programs

  • Maple
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; `if`(n=0,1, add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n) end end: b:= etr(B): c:= etr(b): B:= n-> if n=0 then 0 else c(n-1) fi: C:= etr(B): a:= n-> B(n)+C(n) -add(B(k)*C(n-k), k=0..n): seq(a(n), n=0..30); # Alois P. Heinz, Sep 09 2008
  • Mathematica
    ClearAll[etr, b, a]; etr[p_] := etr[p] = Module[{b}, b[n_] := b[n] = If[n == 0, 1, Sum[ Sum[ d*p[d], {d, Divisors[j]}]*b[n-j], {j, 1, n}]/n]; b]; b[0]=0; b[n_] := b[n] = etr[etr[b]][n-1]; a[n_] := b[n] + etr[b][n] - Sum[b[k]*etr[b][n-k], {k, 0, n}]; Table[ a[n], {n, 0, 27}] (* Jean-François Alcover, Oct 09 2012, after Alois P. Heinz *)
  • PARI
    \\ here b(n) is A007563 as vector
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    b(n)={my(v=[1]);for(i=2, n, v=concat([1], EulerT(EulerT(v)))); v}
    seq(n)={my(u=b(n)); Vec(1 + x*Ser(EulerT(u))*(1-x*Ser(u)))} \\ Andrew Howroyd, May 22 2018

Formula

G.f.: A(x)=1+(C(x)-1)*(1-B(x)). B: G.f. for A007563. C: G.f. for A035052.
a(n) ~ c * d^n / n^(5/2), where d = 4.189610958393826965527036454524... (see A245566), c = 0.245899549044224207821149415964395... . - Vaclav Kotesovec, Jul 26 2014
a(n) = A304937(n) - A304937(n-1) for n>1, a(n) = 1 for n<2. - Gus Wiseman, May 22 2018

A302494 Products of distinct primes of squarefree index.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 13, 15, 17, 22, 26, 29, 30, 31, 33, 34, 39, 41, 43, 47, 51, 55, 58, 59, 62, 65, 66, 67, 73, 78, 79, 82, 83, 85, 86, 87, 93, 94, 101, 102, 109, 110, 113, 118, 123, 127, 129, 130, 134, 137, 139, 141, 143, 145, 146, 149, 155, 157, 158, 163
Offset: 1

Views

Author

Gus Wiseman, Apr 08 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.

Examples

			Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of set systems.
01: {}
02: {{}}
03: {{1}}
05: {{2}}
06: {{},{1}}
10: {{},{2}}
11: {{3}}
13: {{1,2}}
15: {{1},{2}}
17: {{4}}
22: {{},{3}}
26: {{},{1,2}}
29: {{1,3}}
30: {{},{1},{2}}
31: {{5}}
33: {{1},{3}}
34: {{},{4}}
39: {{1},{1,2}}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,SquareFreeQ[#]&&And@@SquareFreeQ/@PrimePi/@FactorInteger[#][[All,1]]]&]
  • PARI
    is(n) = if(bigomega(n)!=omega(n), return(0), my(f=factor(n)[, 1]~); for(k=1, #f, if(!issquarefree(primepi(f[k])) && primepi(f[k])!=1, return(0)))); 1 \\ Felix Fröhlich, Apr 10 2018

A306005 Number of non-isomorphic set-systems of weight n with no singletons.

Original entry on oeis.org

1, 0, 1, 1, 3, 4, 12, 19, 51, 106, 274, 647, 1773, 4664, 13418, 38861, 118690, 370588, 1202924, 4006557, 13764760, 48517672, 175603676, 651026060, 2471150365, 9590103580, 38023295735, 153871104726, 635078474978, 2671365285303, 11444367926725, 49903627379427
Offset: 0

Views

Author

Gus Wiseman, Jun 16 2018

Keywords

Comments

A set-system is a finite set of finite nonempty sets (edges). The weight is the sum of cardinalities of the edges. Weight is generally not the same as number of vertices.

Examples

			Non-isomorphic representatives of the a(6) = 12 set-systems:
  {{1,2,3,4,5,6}}
  {{1,2},{3,4,5,6}}
  {{1,5},{2,3,4,5}}
  {{3,4},{1,2,3,4}}
  {{1,2,3},{4,5,6}}
  {{1,2,5},{3,4,5}}
  {{1,3,4},{2,3,4}}
  {{1,2},{1,3},{2,3}}
  {{1,2},{3,4},{5,6}}
  {{1,2},{3,5},{4,5}}
  {{1,3},{2,4},{3,4}}
  {{1,4},{2,4},{3,4}}
		

Crossrefs

Programs

  • 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)) - Vec(sum(j=1, #q, if(t%q[j]==0, q[j])) + 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

Formula

a(n) = A283877(n) - A330053(n). - Gus Wiseman, Dec 09 2019

Extensions

Terms a(11) and beyond from Andrew Howroyd, Sep 01 2019

A045782 Number of factorizations of n for some n (image of A001055).

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 9, 11, 12, 15, 16, 19, 21, 22, 26, 29, 30, 31, 36, 38, 42, 45, 47, 52, 56, 57, 64, 66, 67, 74, 77, 92, 97, 98, 101, 105, 109, 118, 135, 137, 139, 141, 162, 165, 171, 176, 181, 189, 195, 198, 203, 212, 231, 249, 250, 254, 257, 267, 269, 272, 289
Offset: 1

Views

Author

Keywords

Comments

Also the image of A318284. - Gus Wiseman, Jan 11 2020

Crossrefs

Factorizations are A001055 with image this sequence and complement A330976.
Strict factorizations are A045778 with image A045779 and complement A330975.
The least number with exactly a(n) factorizations is A045783(n).
The least number with exactly n factorizations is A330973(n).

Programs

  • Mathematica
    terms = 61; m0 = 10^5; dm = 10^4;
    f[1, ] = 1; f[n, k_] := f[n, k] = Sum[f[n/d, d], {d, Select[Divisors[n], 1 < # <= k &]}];
    Clear[seq]; seq[m_] := seq[m] = Sort[Tally[Table[f[n, n], {n, 1, m}]][[All, 1]]][[1 ;; terms]]; seq[m = m0]; seq[m += dm]; While[Print[m]; seq[m] != seq[m - dm], m += dm];
    seq[m] (* Jean-François Alcover, Oct 04 2018 *)

Formula

The Luca et al. paper shows that the number of terms with a(n) <= x is x^{ O( log log log x / log log x )}. - N. J. A. Sloane, Jun 12 2009

Extensions

Name edited by Gus Wiseman, Jan 11 2020

A319721 Number of non-isomorphic antichains of multisets of weight n.

Original entry on oeis.org

1, 1, 4, 8, 24, 50, 148, 349, 1014, 2717, 8114
Offset: 0

Views

Author

Gus Wiseman, Sep 26 2018

Keywords

Comments

In an antichain, no part is a proper submultiset of any other. The weight of an antichain is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(3) = 8 antichains:
1: {{1}}
2: {{1,1}}
   {{1,2}}
   {{1},{1}}
   {{1},{2}}
3: {{1,1,1}}
   {{1,2,2}}
   {{1,2,3}}
   {{1},{2,2}}
   {{1},{2,3}}
   {{1},{1},{1}}
   {{1},{2},{2}}
   {{1},{2},{3}}
		

Crossrefs

A368094 Number of non-isomorphic set-systems of weight n contradicting a strict version of the axiom of choice.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 5, 12, 36, 97, 291
Offset: 0

Views

Author

Gus Wiseman, Dec 23 2023

Keywords

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. Weight is generally not the same as number of vertices.
The axiom of choice says that, given any set of nonempty sets Y, it is possible to choose a set containing an element from each. The strict version requires this set to have the same cardinality as Y, meaning no element is chosen more than once.

Examples

			Non-isomorphic representatives of the a(5) = 1 through a(7) = 12 set-systems:
  {{1},{2},{3},{2,3}}  {{1},{2},{1,3},{2,3}}    {{1},{2},{1,2},{3,4,5}}
                       {{1},{2},{3},{1,2,3}}    {{1},{3},{2,3},{1,2,3}}
                       {{2},{3},{1,3},{2,3}}    {{1},{4},{1,4},{2,3,4}}
                       {{3},{4},{1,2},{3,4}}    {{2},{3},{2,3},{1,2,3}}
                       {{1},{2},{3},{4},{3,4}}  {{3},{1,2},{1,3},{2,3}}
                                                {{1},{2},{3},{1,3},{2,3}}
                                                {{1},{2},{3},{2,4},{3,4}}
                                                {{1},{2},{3},{4},{2,3,4}}
                                                {{1},{3},{4},{2,4},{3,4}}
                                                {{1},{4},{5},{2,3},{4,5}}
                                                {{2},{3},{4},{1,2},{3,4}}
                                                {{1},{2},{3},{4},{5},{4,5}}
		

Crossrefs

The case of unlabeled graphs is A140637, complement A134964.
The case of labeled graphs is A367867, complement A133686.
The labeled version is A367903, ranks A367907.
The complement is counted by A368095, connected A368410.
Repeats allowed: A368097, ranks A355529, complement A368098, ranks A368100.
Minimal multiset partitions of this type are ranked by A368187.
The connected case is A368409.
Factorizations of this type are counted by A368413, complement A368414.
Allowing repeated edges gives A368421, complement A368422.
A000110 counts set partitions, non-isomorphic A000041.
A003465 counts covering set-systems, unlabeled A055621.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A058891 counts set-systems, unlabeled A000612, connected A323818.
A283877 counts non-isomorphic set-systems, connected A300913.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]] /@ Cases[Subsets[set],{i,_}];
    mpm[n_]:=Join@@Table[Union[Sort[Sort/@(#/.x_Integer:>s[[x]])]& /@ sps[Range[n]]], {s,Flatten[MapIndexed[Table[#2,{#1}]&,#]]& /@ IntegerPartitions[n]}];
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{i,p[[i]]},{i,Length[p]}])], {p,Permutations[Union@@m]}]]];
    Table[Length[Union[brute/@Select[mpm[n], UnsameQ@@#&&And@@UnsameQ@@@# && Select[Tuples[#], UnsameQ@@#&]=={}&]]],{n,0,8}]

A305148 Number of integer partitions of n whose distinct parts are pairwise indivisible.

Original entry on oeis.org

1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 9, 12, 12, 17, 20, 22, 28, 35, 39, 48, 55, 65, 79, 90, 105, 121, 143, 166, 190, 219, 254, 290, 332, 382, 436, 493, 567, 637, 729, 824, 931, 1052, 1186, 1334, 1504, 1691, 1894, 2123, 2380, 2664, 2968, 3319, 3704, 4119, 4586, 5110
Offset: 0

Views

Author

Gus Wiseman, May 26 2018

Keywords

Examples

			The a(9) = 7 integer partitions are (9), (72), (54), (522), (333), (3222), (111111111).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Select[Tuples[Union[#],2],UnsameQ@@#&&Divisible@@#&]=={}&]],{n,20}]

Extensions

More terms from Alois P. Heinz, May 26 2018

A330098 Number of distinct multisets of multisets that can be obtained by permuting the vertices of the multiset of multisets with MM-number n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 2, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Dec 09 2019

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798. The multiset of multisets with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.
a(n) is a divisor of A303975(n)!.

Examples

			The vertex-permutations of {{1,2},{2,3,3}} are:
  {{1,2},{1,3,3}}
  {{1,2},{2,3,3}}
  {{1,3},{1,2,2}}
  {{1,3},{2,2,3}}
  {{2,3},{1,1,2}}
  {{2,3},{1,1,3}}
so a(4927) = 6.
		

Crossrefs

Positions of 1's are A330232.
Positions of first appearances are A330230 and A330233.
The BII-number version is A330231.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    graprms[m_]:=Union[Table[Sort[Sort/@(m/.Rule@@@Table[{p[[i]],i},{i,Length[p]}])],{p,Permutations[Union@@m]}]];
    Table[Length[graprms[primeMS/@primeMS[n]]],{n,100}]

A367901 Number of sets of subsets of {1..n} contradicting a strict version of the axiom of choice.

Original entry on oeis.org

1, 2, 9, 195, 63765, 4294780073, 18446744073639513336, 340282366920938463463374607341656713953, 115792089237316195423570985008687907853269984665640564039457583610129753447747
Offset: 0

Views

Author

Gus Wiseman, Dec 05 2023

Keywords

Comments

The axiom of choice says that, given any set of nonempty sets Y, it is possible to choose a set containing an element from each. The strict version requires this set to have the same cardinality as Y, meaning no element is chosen more than once.

Examples

			The a(2) = 9 sets of sets:
  {{}}
  {{},{1}}
  {{},{2}}
  {{},{1,2}}
  {{},{1},{2}}
  {{},{1},{1,2}}
  {{},{2},{1,2}}
  {{1},{2},{1,2}}
  {{},{1},{2},{1,2}}
		

Crossrefs

The version for simple graphs is A367867, covering A367868.
The complement is counted by A367902, no singletons A367770, ranks A367906.
The version without empty edges is A367903, ranks A367907.
For a unique choice (instead of none) we have A367904, ranks A367908.
A000372 counts antichains, covering A006126, nonempty A014466.
A003465 counts covering set-systems, unlabeled A055621.
A058891 counts set-systems, unlabeled A000612.
A059201 counts covering T_0 set-systems.
A323818 counts covering connected set-systems, unlabeled A323819.
A326031 gives weight of the set-system with BII-number n.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Subsets[Range[n]]], Select[Tuples[#],UnsameQ@@#&]=={}&]],{n,0,3}]

Formula

a(n) = 2^2^n - A367902(n). - Christian Sievers, Aug 01 2024

Extensions

a(5)-a(8) from Christian Sievers, Aug 01 2024

A035310 Let f(n) = number of ways to factor n = A001055(n); a(n) = sum of f(k) over all terms k in A025487 that have n factors.

Original entry on oeis.org

1, 4, 12, 47, 170, 750, 3255, 16010, 81199, 448156, 2579626, 15913058, 102488024, 698976419, 4976098729, 37195337408, 289517846210, 2352125666883, 19841666995265, 173888579505200, 1577888354510786, 14820132616197925, 143746389756336173, 1438846957477988926
Offset: 1

Views

Author

Keywords

Comments

Ways of partitioning an n-multiset with multiplicities some partition of n.
Number of multiset partitions of strongly normal multisets of size n, where a finite multiset is strongly normal if it covers an initial interval of positive integers with weakly decreasing multiplicities. The (weakly) normal version is A255906. - Gus Wiseman, Dec 31 2019

Examples

			a(3) = 12 because there are 3 terms in A025487 with 3 factors, namely 8, 12, 30; and f(8)=3, f(12)=4, f(30)=5 and 3+4+5 = 12.
From _Gus Wiseman_, Dec 31 2019: (Start)
The a(1) = 1 through a(3) = 12 multiset partitions of strongly normal multisets:
  {{1}}  {{1,1}}    {{1,1,1}}
         {{1,2}}    {{1,1,2}}
         {{1},{1}}  {{1,2,3}}
         {{1},{2}}  {{1},{1,1}}
                    {{1},{1,2}}
                    {{1},{2,3}}
                    {{2},{1,1}}
                    {{2},{1,3}}
                    {{3},{1,2}}
                    {{1},{1},{1}}
                    {{1},{1},{2}}
                    {{1},{2},{3}}
(End)
		

Crossrefs

Sequence A035341 counts the ordered cases. Tables A093936 and A095705 distribute the values; e.g. 81199 = 30 + 536 + 3036 + 6181 + 10726 + 11913 + 14548 + 13082 + 21147.
Row sums of A317449.
The uniform case is A317584.
The case with empty intersection is A317755.
The strict case is A317775.
The constant case is A047968.
The set-system case is A318402.
The case of strict parts is A330783.
Multiset partitions of integer partitions are A001970.
Unlabeled multiset partitions are A007716.

Programs

  • Maple
    with(numtheory):
    g:= proc(n, k) option remember;
          `if`(n>k, 0, 1) +`if`(isprime(n), 0,
          add(`if`(d>k, 0, g(n/d, d)), d=divisors(n) minus {1, n}))
        end:
    b:= proc(n, i, l)
          `if`(n=0, g(mul(ithprime(t)^l[t], t=1..nops(l))$2),
          `if`(i<1, 0, add(b(n-i*j, i-1, [l[], i$j]), j=0..n/i)))
        end:
    a:= n-> b(n$2, []):
    seq(a(n), n=1..10);  # Alois P. Heinz, May 26 2013
  • Mathematica
    g[n_, k_] := g[n, k] = If[n > k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d > k, 0, g[n/d, d]], {d, Divisors[n] ~Complement~ {1, n}}]]; b[n_, i_, l_] := If[n == 0, g[p = Product[Prime[t]^l[[t]], {t, 1, Length[l]}], p], If[i < 1, 0, Sum[b[n - i*j, i-1, Join[l, Array[i&, j]]], {j, 0, n/i}]]]; a[n_] := b[n, n, {}]; Table[Print[an = a[n]]; an, {n, 1, 13}] (* Jean-François Alcover, Dec 12 2013, after Alois P. Heinz *)
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    D(p, n)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); my(u=EulerT(v)); Vec(1/prod(k=1, n, 1 - u[k]*x^k + O(x*x^n))-1, -n)/prod(i=1, #v, i^v[i]*v[i]!)}
    seq(n)={my(s=0); forpart(p=n, s+=D(p,n)); s} \\ Andrew Howroyd, Dec 30 2020
  • Python
    from sympy.core.cache import cacheit
    from sympy import divisors, isprime, prime
    from operator import mul
    @cacheit
    def g(n, k):
        return (0 if n > k else 1) + (0 if isprime(n) else sum(g(n//d, d) for d in divisors(n)[1:-1] if d <= k))
    @cacheit
    def b(n, i, l):
        if n==0:
            p = reduce(mul, (prime(t + 1)**l[t] for t in range(len(l))))
            return g(p, p)
        else:
            return 0 if i<1 else sum([b(n - i*j, i - 1, l + [i]*j) for j in range(n//i + 1)])
    def a(n):
        return b(n, n, [])
    for n in range(1, 11): print(a(n)) # Indranil Ghosh, Aug 19 2017, after Maple code
    

Extensions

More terms from Erich Friedman.
81199 from Alford Arnold, Mar 04 2008
a(10) from Alford Arnold, Mar 31 2008
a(10) corrected by Alford Arnold, Aug 07 2008
a(11)-a(13) from Alois P. Heinz, May 26 2013
a(14) from Alois P. Heinz, Sep 27 2014
a(15) from Alois P. Heinz, Jan 10 2015
Terms a(16) and beyond from Andrew Howroyd, Dec 30 2020
Previous Showing 51-60 of 466 results. Next