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.

Showing 1-10 of 42 results. Next

A303386 Number of aperiodic factorizations of n > 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 2, 3, 1, 4, 1, 4, 2, 2, 1, 7, 1, 2, 2, 4, 1, 5, 1, 6, 2, 2, 2, 7, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 12, 1, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 7, 2, 5, 1, 4, 2, 5, 1, 16, 1, 2, 4, 4, 2, 5, 1, 12, 3, 2, 1, 11, 2, 2, 2, 7, 1, 11, 2, 4, 2, 2, 2, 19, 1, 4, 4, 7, 1, 5, 1, 7, 5
Offset: 2

Views

Author

Gus Wiseman, Apr 23 2018

Keywords

Comments

An aperiodic factorization of n is a finite multiset of positive integers greater than 1 whose product is n and whose multiplicities are relatively prime.

Examples

			The a(36) = 7 aperiodic factorizations are (2*2*9), (2*3*6), (2*18), (3*3*4), (3*12), (4*9), and (36). Missing from this list are (2*2*3*3) and (6*6).
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],GCD@@Length/@Split[#]===1&]],{n,2,100}]
  • PARI
    A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s));
    A052409(n) = { my(k=ispower(n)); if(k, k, n>1); }; \\ From A052409
    A303386(n) = if(1==n,n,my(r); sumdiv(A052409(n),d, ispower(n,d,&r); moebius(d)*A001055(r))); \\ Antti Karttunen, Sep 25 2018

Formula

a(n) = Sum_{d|A052409(n)} mu(d) * A001055(n^(1/d)), where mu = A008683.

Extensions

More terms from Antti Karttunen, Sep 25 2018

A303431 Aperiodic tree numbers. Matula-Goebel numbers of aperiodic rooted trees.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 12, 13, 15, 18, 20, 22, 24, 26, 29, 30, 31, 33, 37, 39, 40, 41, 44, 45, 47, 48, 50, 52, 54, 55, 58, 60, 61, 62, 65, 66, 71, 72, 74, 75, 78, 79, 80, 82, 87, 88, 89, 90, 93, 94, 96, 99, 101, 104, 108, 109, 110, 111, 113, 116, 117, 120, 122
Offset: 1

Views

Author

Gus Wiseman, Apr 23 2018

Keywords

Comments

A positive integer is an aperiodic tree number iff either it is equal to 1 or it belongs to A007916 (numbers that are not perfect powers, or numbers whose prime multiplicities are relatively prime) and all of its prime indices are also aperiodic tree numbers, where a prime index of n is a number m such that prime(m) divides n.

Examples

			Sequence of aperiodic rooted trees begins:
01 o
02 (o)
03 ((o))
05 (((o)))
06 (o(o))
10 (o((o)))
11 ((((o))))
12 (oo(o))
13 ((o(o)))
15 ((o)((o)))
18 (o(o)(o))
20 (oo((o)))
22 (o(((o))))
24 (ooo(o))
26 (o(o(o)))
29 ((o((o))))
30 (o(o)((o)))
31 (((((o)))))
33 ((o)(((o))))
		

Crossrefs

Programs

  • Mathematica
    zapQ[1]:=True;zapQ[n_]:=And[GCD@@FactorInteger[n][[All,2]]===1,And@@zapQ/@PrimePi/@FactorInteger[n][[All,1]]];
    Select[Range[100],zapQ]

A303546 Number of non-isomorphic aperiodic multiset partitions of weight n.

Original entry on oeis.org

1, 3, 9, 29, 90, 285, 909, 2984, 9935, 34113, 119368, 428923, 1574223, 5915235, 22699730, 89000042, 356058539, 1453069854, 6044132793, 25612564200, 110503626702, 485161228675, 2166488899641, 9835209480533, 45370059225227
Offset: 1

Views

Author

Gus Wiseman, Apr 26 2018

Keywords

Comments

A multiset is aperiodic if its multiplicities are relatively prime. For this sequence neither the parts nor their multiset union are required to be aperiodic, only the multiset of parts.

Examples

			Non-isomorphic representatives of the a(3) = 9 aperiodic multiset partitions are:
  {{1,1,1}}, {{1,2,2}}, {{1,2,3}},
  {{1},{1,1}}, {{1},{2,2}}, {{1},{2,3}}, {{2},{1,2}},
  {{1},{2},{2}}, {{1},{2},{3}}.
		

Crossrefs

Formula

a(n) = Sum_{d|n} mu(d) * A007716(n/d).

A306200 Number of unlabeled rooted semi-identity trees with n nodes.

Original entry on oeis.org

0, 1, 1, 2, 4, 8, 18, 41, 98, 237, 591, 1488, 3805, 9820, 25593, 67184, 177604, 472177, 1261998, 3388434, 9136019, 24724904, 67141940, 182892368, 499608724, 1368340326, 3756651116, 10336434585, 28499309291, 78727891420, 217870037932, 603934911859, 1676720329410
Offset: 0

Views

Author

Gus Wiseman, Jan 29 2019

Keywords

Comments

A rooted tree is a semi-identity tree if the non-leaf branches of the root are all distinct and are themselves semi-identity trees.

Examples

			The a(1) = 1 through a(7) = 8 trees:
  o  (o)  (oo)   (ooo)    (oooo)     (ooooo)
          ((o))  ((oo))   ((ooo))    ((oooo))
                 (o(o))   (o(oo))    (o(ooo))
                 (((o)))  (oo(o))    (oo(oo))
                          (((oo)))   (ooo(o))
                          ((o(o)))   (((ooo)))
                          (o((o)))   ((o)(oo))
                          ((((o))))  ((o(oo)))
                                     ((oo(o)))
                                     (o((oo)))
                                     (o(o(o)))
                                     (oo((o)))
                                     ((((oo))))
                                     (((o(o))))
                                     ((o)((o)))
                                     ((o((o))))
                                     (o(((o))))
                                     (((((o)))))
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,
          add(b(n-i*j, i-1)*binomial(a(i), j), j=0..n/i))
        end:
    a:= n-> `if`(n=0, 0, b(n-1$2)):
    seq(a(n), n=0..35);  # Alois P. Heinz, Jan 29 2019
  • Mathematica
    ursit[n_]:=Join@@Table[Select[Union[Sort/@Tuples[ursit/@ptn]],UnsameQ@@DeleteCases[#,{}]&],{ptn,IntegerPartitions[n-1]}];
    Table[Length[ursit[n]],{n,10}]
    (* Second program: *)
    b[n_, i_] := b[n, i] = If[n == 0 || i == 1, 1,
         Sum[b[n - i*j, i - 1]*Binomial[a[i], j], {j, 0, n/i}]];
    a[n_] := If[n == 0, 0, b[n - 1, n - 1]];
    a /@ Range[0, 35] (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)

Extensions

More terms from Alois P. Heinz, Jan 29 2019

A317712 Number of uniform rooted trees with n nodes.

Original entry on oeis.org

1, 1, 2, 4, 8, 15, 35, 72, 169, 388, 934, 2234, 5508, 13557, 33883, 85017, 215091, 546496, 1396524, 3582383, 9228470, 23852918, 61857180, 160871716, 419516462, 1096671326, 2873403980, 7544428973, 19847520789, 52308750878, 138095728065, 365153263313, 966978876376
Offset: 1

Views

Author

Gus Wiseman, Aug 05 2018

Keywords

Comments

An unlabeled rooted tree is uniform if the multiplicities of the branches directly under any given node are all equal.

Examples

			The a(5) = 8 uniform rooted trees:
  ((((o))))
  (((oo)))
  ((o(o)))
  ((ooo))
  (o((o)))
  (o(oo))
  ((o)(o))
  (oooo)
		

Crossrefs

Programs

  • Mathematica
    purt[n_]:=Join@@Table[Select[Union[Sort/@Tuples[purt/@ptn]],SameQ@@Length/@Split[#]&],{ptn,IntegerPartitions[n-1]}];
    Table[Length[purt[n]],{n,10}]
  • PARI
    WeighT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v,n,(-1)^(n-1)/n))))-1,-#v)}
    seq(n)={my(v=[1]); for(n=2, n, my(t=WeighT(v)); v=concat(v, sumdiv(n-1, d, t[d]))); v} \\ Andrew Howroyd, Aug 28 2018

Formula

a(n) ~ c * d^n / n^(3/2), where d = 2.774067238136373782458114960391469140405537808253... and c = 0.43338208953061974806801546569720246018271214... - Vaclav Kotesovec, Sep 07 2019

Extensions

Term a(21) and beyond from Andrew Howroyd, Aug 28 2018

A324765 Number of recursively anti-transitive rooted trees with n nodes.

Original entry on oeis.org

1, 1, 2, 3, 6, 11, 26, 52, 119, 266, 618, 1432, 3402, 8093, 19505, 47228, 115244, 282529, 696388, 1723400
Offset: 1

Views

Author

Gus Wiseman, Mar 17 2019

Keywords

Comments

An unlabeled rooted tree is recursively anti-transitive if no branch of a branch of a terminal subtree is a branch of the same subtree.

Examples

			The a(1) = 1 through a(6) = 11 recursively anti-transitive rooted trees:
  o  (o)  (oo)   (ooo)    (oooo)     (ooooo)
          ((o))  ((oo))   ((ooo))    ((oooo))
                 (((o)))  (((oo)))   (((ooo)))
                          ((o)(o))   ((o)(oo))
                          (o((o)))   (o((oo)))
                          ((((o))))  (oo((o)))
                                     ((((oo))))
                                     (((o)(o)))
                                     ((o((o))))
                                     (o(((o))))
                                     (((((o)))))
		

Crossrefs

Programs

  • Mathematica
    nallt[n_]:=Select[Union[Sort/@Join@@(Tuples[nallt/@#]&/@IntegerPartitions[n-1])],Intersection[Union@@#,#]=={}&];
    Table[Length[nallt[n]],{n,10}]

A317707 Number of powerful rooted trees with n nodes.

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 11, 13, 22, 29, 46, 57, 94, 115, 180, 230, 349, 435, 671, 830, 1245, 1572, 2320, 2894, 4287, 5328, 7773, 9752, 14066, 17547, 25328, 31515, 45010, 56289, 79805, 99467, 140778, 175215, 246278, 307273, 429421, 534774, 745776, 927776, 1287038
Offset: 1

Views

Author

Gus Wiseman, Aug 05 2018

Keywords

Comments

An unlabeled rooted tree is powerful if either it is a single node or a single node with a single powerful tree as a branch, or if the branches of the root all appear with multiplicities greater than 1 and are themselves powerful trees.

Examples

			The a(7) = 11 powerful rooted trees:
  ((((((o))))))
  (((((oo)))))
  ((((ooo))))
  ((((o)(o))))
  (((oooo)))
  ((ooooo))
  (((o))((o)))
  ((oo)(oo))
  ((o)(o)(o))
  (oo(o)(o))
  (oooooo)
		

Crossrefs

Programs

  • Maple
    h:= proc(n, k, t) option remember; `if`(k=0, binomial(n+t, t),
          `if`(n=0, 0, add(h(n-1, k-j, t+1), j=2..k)))
        end:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(b(n-i*j, i-1)*h(a(i), j, 0), j=0..n/i)))
        end:
    a:= proc(n) option remember; `if`(n<2, n, b(n-1$2)+a(n-1)) end:
    seq(a(n), n=1..50);  # Alois P. Heinz, Aug 31 2018
  • Mathematica
    purt[n_]:=If[n==1,{{}},Join@@Table[Select[Union[Sort/@Tuples[purt/@ptn]],Or[Length[#]==1,Min@@Length/@Split[#]>1]&],{ptn,IntegerPartitions[n-1]}]];
    Table[Length[purt[n]],{n,10}]
    (* Second program: *)
    h[n_, k_, t_] := h[n, k, t] = If[k == 0, Binomial[n + t, t], If[n == 0, 0, Sum[h[n - 1, k - j, t + 1], {j, 2, k}]]];
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0, Sum[b[n - i*j, i - 1]* h[a[i], j, 0], {j, 0, n/i}]]];
    a[n_] := a[n] = If[n < 2, n, b[n - 1, n - 1] + a[n - 1]];
    Array[a, 50] (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)

Extensions

a(27)-a(45) from Alois P. Heinz, Aug 31 2018

A317708 Number of aperiodic relatively prime trees with n nodes.

Original entry on oeis.org

1, 1, 1, 2, 4, 10, 20, 48, 108, 255, 595, 1435, 3434, 8372, 20419, 50289, 124289, 309122, 771508, 1934462
Offset: 1

Views

Author

Gus Wiseman, Aug 05 2018

Keywords

Comments

An unlabeled rooted tree is aperiodic and relatively prime iff either it is a single node or a single node with a single aperiodic relatively prime branch, or the branches directly under any given node have empty intersection (relatively prime) and also have relatively prime multiplicities (aperiodic) and are themselves aperiodic relatively prime trees.

Examples

			The a(6) = 10 aperiodic relatively prime trees:
  (((((o)))))
  (((o(o))))
  ((o((o))))
  ((oo(o)))
  (o(((o))))
  (o(o(o)))
  ((o)((o)))
  (oo((o)))
  (o(o)(o))
  (ooo(o))
		

Crossrefs

Programs

  • Mathematica
    rurt[n_]:=If[n==1,{{}},Join@@Table[Select[Union[Sort/@Tuples[rurt/@ptn]],Or[Length[#]==1,And[Intersection@@#=={},GCD@@Length/@Split[#]==1]]&],{ptn,IntegerPartitions[n-1]}]];
    Table[Length[rurt[n]],{n,10}]

A320813 Number of non-isomorphic multiset partitions of an aperiodic multiset of weight n such that there are no singletons and all parts are themselves aperiodic multisets.

Original entry on oeis.org

1, 0, 1, 2, 5, 13, 33, 104, 293, 938, 2892
Offset: 0

Views

Author

Gus Wiseman, Nov 08 2018

Keywords

Comments

Also the number of nonnegative integer matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, in which (1) the row sums are all > 1, (2) the positive entries in each row are relatively prime, and (3) the column-sums are relatively prime.
A multiset is aperiodic if its multiplicities are relatively prime.
The weight of a multiset partition 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(2) = 1 through a(5) = 13 multiset partitions:
  {{1,2}}  {{1,2,2}}  {{1,2,2,2}}    {{1,1,2,2,2}}
           {{1,2,3}}  {{1,2,3,3}}    {{1,2,2,2,2}}
                      {{1,2,3,4}}    {{1,2,2,3,3}}
                      {{1,2},{3,4}}  {{1,2,3,3,3}}
                      {{1,3},{2,3}}  {{1,2,3,4,4}}
                                     {{1,2,3,4,5}}
                                     {{1,2},{1,2,2}}
                                     {{1,2},{2,3,3}}
                                     {{1,2},{3,4,4}}
                                     {{1,2},{3,4,5}}
                                     {{1,3},{2,3,3}}
                                     {{1,4},{2,3,4}}
                                     {{2,3},{1,2,3}}
		

Crossrefs

This is the case of A320804 where the underlying multiset is aperiodic.

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]}]]];
    aperQ[m_]:=Length[m]==0||GCD@@Length/@Split[Sort[m]]==1;
    Table[Length[Union[brute /@ Select[mpm[n],And[Min@@Length/@#>1,aperQ[Join@@#]&&And@@aperQ /@ #]&]]],{n,0,7}] (* Gus Wiseman, Jan 19 2024 *)

Extensions

Definition corrected by Gus Wiseman, Jan 19 2024

A303708 Number of aperiodic factorizations of n using elements of A007916 (numbers that are not perfect powers).

Original entry on oeis.org

0, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 3, 1, 2, 2, 0, 1, 3, 1, 3, 2, 2, 1, 4, 0, 2, 0, 3, 1, 5, 1, 0, 2, 2, 2, 3, 1, 2, 2, 4, 1, 5, 1, 3, 3, 2, 1, 5, 0, 3, 2, 3, 1, 4, 2, 4, 2, 2, 1, 9, 1, 2, 3, 0, 2, 5, 1, 3, 2, 5, 1, 8, 1, 2, 3, 3, 2, 5, 1, 5, 0, 2, 1, 9, 2, 2, 2, 4, 1, 9, 2
Offset: 1

Views

Author

Gus Wiseman, Apr 29 2018

Keywords

Comments

An aperiodic factorization of n is a finite multiset of positive integers greater than 1 whose product is n and whose multiplicities are relatively prime.
The positions of zeros in this sequence are the prime powers A000961.

Examples

			The a(144) = 8 aperiodic factorizations are (2*2*2*3*6), (2*2*2*18), (2*2*3*12), (2*3*24), (2*6*12), (2*72), (3*48) and (6*24). Missing from this list are (12*12), (2*2*6*6) and (2*2*2*2*3*3).
		

Crossrefs

Programs

  • Mathematica
    radQ[n_]:=Or[n===1,GCD@@FactorInteger[n][[All,2]]===1];
    facsr[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsr[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],radQ]}]];
    Table[Length[Select[facsr[n],GCD@@Length/@Split[#]===1&]],{n,100}]

Formula

a(n) = Sum_{d in A007916, d|A052409(n)} mu(d) * A303707(n^(1/d)).
Showing 1-10 of 42 results. Next