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-6 of 6 results.

A324844 Number of unlabeled rooted trees with n nodes where the branches of no non-leaf branch of any terminal subtree form a submultiset of the branches of the same subtree.

Original entry on oeis.org

1, 1, 2, 3, 7, 13, 32, 71, 170, 406, 1002, 2469, 6204, 15644, 39871, 102116, 263325, 682079, 1775600, 4640220
Offset: 1

Views

Author

Gus Wiseman, Mar 18 2019

Keywords

Examples

			The a(1) = 1 through a(6) = 13 rooted trees:
  o  (o)  (oo)   (ooo)    (oooo)     (ooooo)
          ((o))  ((oo))   ((ooo))    ((oooo))
                 (((o)))  (o(oo))    (o(ooo))
                          (((oo)))   (((ooo)))
                          ((o)(o))   ((o)(oo))
                          (o((o)))   ((o(oo)))
                          ((((o))))  (o((oo)))
                                     (oo((o)))
                                     ((((oo))))
                                     (((o)(o)))
                                     ((o((o))))
                                     (o(((o))))
                                     (((((o)))))
		

Crossrefs

The Matula-Goebel numbers of these trees are given by A324845.

Programs

  • Mathematica
    submultQ[cap_,fat_]:=And@@Function[i,Count[fat,i]>=Count[cap,i]]/@Union[List@@cap];
    rallt[n_]:=Select[Union[Sort/@Join@@(Tuples[rallt/@#]&/@IntegerPartitions[n-1])],And@@Table[!submultQ[b,#],{b,DeleteCases[#,{}]}]&];
    Table[Length[rallt[n]],{n,10}]

A324744 Number of maximal subsets of {1...n} containing no element whose prime indices all belong to the subset.

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 4, 5, 6, 8, 8, 11, 11, 22, 22, 22, 22, 28, 28, 44, 44, 52, 52, 76, 76, 88, 88, 96, 96, 184, 184, 240, 240, 264, 264, 296, 296, 592, 592, 592, 592, 728, 728, 1456, 1456, 1456, 1456, 2912, 2912, 3168, 3168, 3168, 3168, 5568, 5568, 5568, 5568
Offset: 0

Views

Author

Gus Wiseman, Mar 15 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.

Examples

			The a(1) = 1 through a(8) = 6 maximal subsets:
  {1}  {1}  {2}    {1,3}  {1,3}    {1,3,6}    {3,4,6}    {1,3,6,7}
       {2}  {1,3}  {2,4}  {1,5}    {1,5,6}    {1,3,6,7}  {1,5,6,7}
                   {3,4}  {3,4}    {3,4,6}    {1,5,6,7}  {3,4,6,8}
                          {2,4,5}  {2,4,5,6}  {2,4,5,6}  {3,6,7,8}
                                              {2,5,6,7}  {2,4,5,6,8}
                                                         {2,5,6,7,8}
		

Crossrefs

The non-maximal case is A324738. The case for subsets of {2...n} is A324762.

Programs

  • Mathematica
    maxim[s_]:=Complement[s,Last/@Select[Tuples[s,2],UnsameQ@@#&&SubsetQ@@#&]];
    Table[Length[maxim[Select[Subsets[Range[n]],!MemberQ[#,k_/;SubsetQ[#,PrimePi/@First/@FactorInteger[k]]]&]]],{n,0,10}]
  • PARI
    pset(n)={my(b=0, f=factor(n)[, 1]); sum(i=1, #f, 1<<(primepi(f[i])))}
    a(n)={my(p=vector(n, k, if(k==1, 1, pset(k))), d=0); for(i=1, #p, d=bitor(d, p[i]));
    my(ismax(b)=for(k=1, #p, if(!bittest(b,k) && bitnegimply(p[k], b), my(e=bitor(b, 1<#p, ismax(b), my(f=bitnegimply(p[k], b)); if(!f || bittest(d, k), self()(k+1, b)) + if(f, self()(k+1, b+(1<Andrew Howroyd, Aug 27 2019

Extensions

Terms a(16) and beyond from Andrew Howroyd, Aug 27 2019

A324738 Number of subsets of {1...n} containing no element > 1 whose prime indices all belong to the subset.

Original entry on oeis.org

1, 2, 3, 5, 8, 13, 26, 42, 72, 120, 232, 376, 752, 1128, 2256, 4512, 8256, 13632, 27264, 42048, 82944, 158976, 313344, 497664, 995328, 1700352, 3350016, 5815296, 11630592, 17491968, 34983936, 56954880, 108933120, 210788352, 418258944, 804667392, 1609334784
Offset: 0

Views

Author

Gus Wiseman, Mar 13 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.

Examples

			The a(0) = 1 through a(6) = 26 subsets:
  {}  {}   {}   {}     {}     {}       {}
      {1}  {1}  {1}    {1}    {1}      {1}
           {2}  {2}    {2}    {2}      {2}
                {3}    {3}    {3}      {3}
                {1,3}  {4}    {4}      {4}
                       {1,3}  {5}      {5}
                       {2,4}  {1,3}    {6}
                       {3,4}  {1,5}    {1,3}
                              {2,4}    {1,5}
                              {2,5}    {1,6}
                              {3,4}    {2,4}
                              {4,5}    {2,5}
                              {2,4,5}  {2,6}
                                       {3,4}
                                       {3,6}
                                       {4,5}
                                       {4,6}
                                       {5,6}
                                       {1,3,6}
                                       {1,5,6}
                                       {2,4,5}
                                       {2,4,6}
                                       {2,5,6}
                                       {3,4,6}
                                       {4,5,6}
                                       {2,4,5,6}
		

Crossrefs

The maximal case is A324744. The case of subsets of {2...n} is A324739. The strict integer partition version is A324749. The integer partition version is A324754. The Heinz number version is A324759. An infinite version is A324694.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],!MemberQ[#,k_/;SubsetQ[#,PrimePi/@First/@FactorInteger[k]]]&]],{n,0,10}]
  • PARI
    pset(n)={my(b=0,f=factor(n)[,1]); sum(i=1, #f, 1<<(primepi(f[i])))}
    a(n)={my(p=vector(n,k,if(k==1, 1, pset(k))), d=0); for(i=1, #p, d=bitor(d, p[i]));
    ((k,b)->if(k>#p, 1, my(t=self()(k+1,b)); if(bitnegimply(p[k], b), t+=if(bittest(d,k), self()(k+1, b+(1<Andrew Howroyd, Aug 16 2019

Extensions

Terms a(21) and beyond from Andrew Howroyd, Aug 16 2019

A324759 Heinz numbers of integer partitions containing no part > 1 whose prime indices all belong to the partition.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 13, 16, 17, 19, 20, 21, 22, 23, 25, 26, 27, 29, 31, 32, 33, 34, 35, 37, 39, 40, 41, 43, 44, 46, 47, 49, 50, 51, 52, 53, 57, 58, 59, 61, 62, 63, 64, 65, 67, 68, 71, 73, 74, 77, 79, 80, 81, 82, 83, 85, 86, 87, 88, 89, 91, 92, 93
Offset: 1

Views

Author

Gus Wiseman, Mar 17 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 Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
   1: {}
   2: {1}
   3: {2}
   4: {1,1}
   5: {3}
   7: {4}
   8: {1,1,1}
   9: {2,2}
  10: {1,3}
  11: {5}
  13: {6}
  16: {1,1,1,1}
  17: {7}
  19: {8}
  20: {1,1,3}
  21: {2,4}
  22: {1,5}
  23: {9}
  25: {3,3}
  26: {1,6}
		

Crossrefs

The subset version is A324738, with maximal case A324744. The strict integer partition version is A324749. The integer partition version is A324754. An infinite version is A324694.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],!MemberQ[DeleteCases[primeMS[#],1],k_/;SubsetQ[primeMS[#],primeMS[k]]]&]

A324754 Number of integer partitions of n containing no part > 1 whose prime indices all belong to the partition.

Original entry on oeis.org

1, 1, 2, 2, 4, 3, 7, 8, 11, 12, 19, 19, 30, 34, 46, 50, 71, 76, 104, 119, 151, 171, 225, 247, 315, 360, 446, 504, 629, 703, 867, 986, 1192, 1346, 1636, 1837, 2204, 2500, 2965, 3348, 3980, 4475, 5276, 5963, 6973, 7852, 9194, 10335, 12009, 13536, 15650, 17589
Offset: 0

Views

Author

Gus Wiseman, Mar 16 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.
For example, (6,2) is such a partition because the prime indices of 6 are {1,2}, which do not all belong to the partition. On the other hand, (5,3) is not such a partition because the prime indices of 5 are {3}, and 3 belongs to the partition.

Examples

			The a(1) = 1 through a(8) = 11  integer partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (111)  (22)    (311)    (33)      (43)       (44)
                    (31)    (11111)  (42)      (52)       (62)
                    (1111)           (51)      (61)       (71)
                                     (222)     (331)      (422)
                                     (3111)    (511)      (611)
                                     (111111)  (31111)    (2222)
                                               (1111111)  (3311)
                                                          (5111)
                                                          (311111)
                                                          (11111111)
		

Crossrefs

The subset version is A324738, with maximal case A324744. The strict case is A324749. The Heinz number version is A324759. An infinite version is A324694.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],!MemberQ[#,k_/;SubsetQ[#,PrimePi/@First/@FactorInteger[k]]]&]],{n,0,30}]

A324845 Matula-Goebel numbers of rooted trees where the branches of no non-leaf branch of any terminal subtree form a submultiset of the branches of the same subtree.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 14, 16, 17, 19, 20, 21, 22, 23, 25, 27, 29, 31, 32, 33, 34, 35, 38, 40, 43, 44, 46, 49, 50, 51, 53, 57, 58, 59, 62, 63, 64, 67, 68, 69, 70, 71, 73, 76, 77, 79, 80, 81, 83, 85, 86, 87, 88, 92, 93, 95, 97, 98, 99, 100, 103, 106
Offset: 1

Views

Author

Gus Wiseman, Mar 18 2019

Keywords

Examples

			The sequence of terms together with their Matula-Goebel numbers begins:
   1: o
   2: (o)
   3: ((o))
   4: (oo)
   5: (((o)))
   7: ((oo))
   8: (ooo)
   9: ((o)(o))
  10: (o((o)))
  11: ((((o))))
  14: (o(oo))
  16: (oooo)
  17: (((oo)))
  19: ((ooo))
  20: (oo((o)))
  21: ((o)(oo))
  22: (o(((o))))
  23: (((o)(o)))
  25: (((o))((o)))
  27: ((o)(o)(o))
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    qaQ[n_]:=And[And@@Table[!Divisible[n,x],{x,DeleteCases[primeMS[n],1]}],And@@qaQ/@primeMS[n]];
    Select[Range[100],qaQ]
Showing 1-6 of 6 results.