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 31-40 of 146 results. Next

A317710 Uniform tree numbers. Matula-Goebel numbers of uniform rooted trees.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, 38, 39, 41, 42, 43, 46, 47, 49, 51, 53, 55, 57, 58, 59, 62, 64, 65, 66, 67, 69, 70, 73, 77, 78, 79, 81, 82, 83, 85, 86, 87, 91, 93, 94, 95, 97
Offset: 1

Views

Author

Gus Wiseman, Aug 05 2018

Keywords

Comments

A positive integer n is a uniform tree number iff either n = 1 or n is a power of a squarefree number whose prime indices are also uniform tree numbers. A prime index of n is a number m such that prime(m) divides n.

Crossrefs

Programs

  • Mathematica
    rupQ[n_]:=Or[n==1,And[SameQ@@FactorInteger[n][[All,2]],And@@rupQ/@PrimePi/@FactorInteger[n][[All,1]]]];
    Select[Range[100],rupQ]

A324525 Numbers divisible by prime(k)^k for each prime index k.

Original entry on oeis.org

1, 2, 4, 8, 9, 16, 18, 27, 32, 36, 54, 64, 72, 81, 108, 125, 128, 144, 162, 216, 243, 250, 256, 288, 324, 432, 486, 500, 512, 576, 625, 648, 729, 864, 972, 1000, 1024, 1125, 1152, 1250, 1296, 1458, 1728, 1944, 2000, 2048, 2187, 2250, 2304, 2401, 2500, 2592
Offset: 1

Views

Author

Gus Wiseman, Mar 08 2019

Keywords

Comments

These are a kind of self-describing numbers (cf. A001462, A304679).
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 prime signature of a number is the multiset of multiplicities (or exponents) in its prime factorization.
Also Heinz numbers of integer partitions where the multiplicity of k is at least k (A117144). 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 as follows. For example, 36 = prime(1) * prime(1) * prime(2) * prime(2) is a term because the prime multiplicities are {2,2}, which are greater than or equal to the prime indices {1,2}.
    1: {}
    2: {1}
    4: {1,1}
    8: {1,1,1}
    9: {2,2}
   16: {1,1,1,1}
   18: {1,2,2}
   27: {2,2,2}
   32: {1,1,1,1,1}
   36: {1,1,2,2}
   54: {1,2,2,2}
   64: {1,1,1,1,1,1}
   72: {1,1,1,2,2}
   81: {2,2,2,2}
  108: {1,1,2,2,2}
  125: {3,3,3}
  128: {1,1,1,1,1,1,1}
		

Crossrefs

Sequences related to self-description: A000002, A001462, A079000, A079254, A276625, A304360.

Programs

  • Maple
    q:= n-> andmap(i-> i[2]>=numtheory[pi](i[1]), ifactors(n)[2]):
    select(q, [$1..3000])[];  # Alois P. Heinz, Mar 08 2019
  • Mathematica
    Select[Range[1000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>k>=PrimePi[p]]&]
    seq[max_] := Module[{ps = {2}, p, s = {1}, s1, s2, emax}, While[ps[[-1]]^Length[ps] < max, AppendTo[ps, NextPrime[ps[[-1]]]]]; Do[p = ps[[k]]; emax = Floor[Log[p, max]]; s1 = Join[{1}, p^Range[k, emax]]; s2 = Select[Union[Flatten[Outer[Times, s, s1]]], # <= max &]; s = Union[s, s2], {k, 1, Length[ps]}]; s]; seq[3000] (* Amiram Eldar, Nov 23 2020 *)

Formula

Closed under multiplication.
Sum_{n>=1} 1/a(n) = Product_{k>=1} 1 + 1/(prime(k)^(k-1) * (prime(k)-1)) = 2.35782843100111139159... - Amiram Eldar, Nov 23 2020

A324697 Lexicographically earliest sequence of positive integers > 1 that are prime or whose prime indices already belong to the sequence.

Original entry on oeis.org

2, 3, 5, 7, 9, 11, 13, 15, 17, 19, 23, 25, 27, 29, 31, 33, 37, 41, 43, 45, 47, 51, 53, 55, 59, 61, 67, 69, 71, 73, 75, 79, 81, 83, 85, 89, 93, 97, 99, 101, 103, 107, 109, 113, 115, 121, 123, 125, 127, 131, 135, 137, 139, 141, 149, 151, 153, 155, 157, 163, 165
Offset: 1

Views

Author

Gus Wiseman, Mar 10 2019

Keywords

Comments

A self-describing sequence, similar to A304360.
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 sequence of terms together with their prime indices begins:
   2: {1}
   3: {2}
   5: {3}
   7: {4}
   9: {2,2}
  11: {5}
  13: {6}
  15: {2,3}
  17: {7}
  19: {8}
  23: {9}
  25: {3,3}
  27: {2,2,2}
  29: {10}
  31: {11}
  33: {2,5}
  37: {12}
  41: {13}
  43: {14}
  45: {2,2,3}
		

Crossrefs

Programs

  • Mathematica
    aQ[n_]:=Switch[n,1,False,?PrimeQ,True,,And@@Cases[FactorInteger[n],{p_,k_}:>aQ[PrimePi[p]]]];
    Select[Range[100],aQ]

A306202 Matula-Goebel numbers of rooted semi-identity trees.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 16, 17, 19, 20, 21, 22, 24, 26, 28, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 47, 48, 51, 52, 53, 55, 56, 57, 58, 59, 60, 62, 64, 65, 66, 67, 68, 70, 71, 73, 74, 76, 77, 78, 79, 80, 82, 84, 85
Offset: 1

Views

Author

Gus Wiseman, Jan 29 2019

Keywords

Comments

Definition: A positive integer belongs to the sequence iff its prime indices greater than 1 are distinct and already belong to the sequence. 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 sequence of all unlabeled rooted semi-identity trees together with their Matula-Goebel numbers begins:
   1: o
   2: (o)
   3: ((o))
   4: (oo)
   5: (((o)))
   6: (o(o))
   7: ((oo))
   8: (ooo)
  10: (o((o)))
  11: ((((o))))
  12: (oo(o))
  13: ((o(o)))
  14: (o(oo))
  15: ((o)((o)))
  16: (oooo)
  17: (((oo)))
  19: ((ooo))
  20: (oo((o)))
  21: ((o)(oo))
  22: (o(((o))))
  24: (ooo(o))
  26: (o(o(o)))
  28: (oo(oo))
  29: ((o((o))))
  30: (o(o)((o)))
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    psidQ[n_]:=And[UnsameQ@@DeleteCases[primeMS[n],1],And@@psidQ/@primeMS[n]];
    Select[Range[100],psidQ]

A324743 Number of maximal subsets of {1...n} containing no prime indices of the elements.

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 5, 8, 8, 8, 8, 12, 12, 18, 18, 19, 19, 30, 30, 54, 54, 54, 54, 96, 96, 96, 96, 96, 96, 156, 156, 244, 244, 248, 248, 248, 248, 440, 440, 440, 440, 688, 688, 1120, 1120, 1120, 1120, 1864, 1864, 1864, 1864, 1864, 1864, 3664, 3664, 3664, 3664, 3664
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(0) = 1 through a(8) = 8 maximal subsets:
  {}  {1}  {1}  {2}    {1,3}  {1,3}    {1,3}    {1,3,7}  {1,3,7}
           {2}  {1,3}  {2,4}  {1,5}    {1,5}    {1,5,7}  {1,5,7}
                       {3,4}  {3,4}    {2,4,5}  {2,4,5}  {2,4,5,8}
                              {2,4,5}  {3,4,6}  {2,5,7}  {2,5,7,8}
                                       {4,5,6}  {3,4,6}  {3,4,6,8}
                                                {3,6,7}  {3,6,7,8}
                                                {4,5,6}  {4,5,6,8}
                                                {5,6,7}  {5,6,7,8}
An example for n = 15 is {1,5,7,9,13,15}, with prime indices:
  1: {}
  5: {3}
  7: {4}
  9: {2,2}
  13: {6}
  15: {2,3}
None of these prime indices {2,3,4,6} belong to the subset, as required.
		

Crossrefs

The non-maximal case is A324741. The case for subsets of {2...n} is A324763.

Programs

  • Mathematica
    maxim[s_]:=Complement[s,Last/@Select[Tuples[s,2],UnsameQ@@#&&SubsetQ@@#&]];
    Table[Length[maxim[Select[Subsets[Range[n]],Intersection[#,PrimePi/@First/@Join@@FactorInteger/@#]=={}&]]],{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, pset(k)), d=0); for(i=1, #p, d=bitor(d, p[i]));
    my(ismax(b)=my(e=0); forstep(k=#p, 1, -1, if(bittest(b,k), e=bitor(e,p[k]), if(!bittest(e,k) && !bitand(p[k], b), return(0)) )); 1);
    ((k, b)->if(k>#p, ismax(b), my(f=!bitand(p[k], b)); if(!f || bittest(d, k), self()(k+1, b)) + if(f, self()(k+1, b+(1<Andrew Howroyd, Aug 26 2019

Extensions

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

A324753 Number of integer partitions of n containing all prime indices of their parts.

Original entry on oeis.org

1, 1, 1, 2, 2, 4, 5, 7, 8, 14, 16, 23, 29, 40, 49, 66, 81, 109, 133, 172, 211, 274, 332, 419, 511, 640, 775, 965, 1165, 1434, 1730, 2109, 2530, 3083, 3683, 4447, 5308, 6375, 7573, 9062, 10730, 12786, 15104, 17909, 21095, 24937, 29284, 34488, 40421, 47450
Offset: 0

Views

Author

Gus Wiseman, Mar 16 2019

Keywords

Comments

These could be described as transitive integer partitions.
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) = 8 integer partitions:
  (1)  (11)  (21)   (211)   (41)     (321)     (421)      (3221)
             (111)  (1111)  (221)    (411)     (2221)     (4211)
                            (2111)   (2211)    (3211)     (22211)
                            (11111)  (21111)   (4111)     (32111)
                                     (111111)  (22111)    (41111)
                                               (211111)   (221111)
                                               (1111111)  (2111111)
                                                          (11111111)
		

Crossrefs

The subset version is A324736. The strict case is A324748. The Heinz number version is A290822. An infinite version is A324698.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],SubsetQ[#,PrimePi/@First/@Join@@FactorInteger/@DeleteCases[#,1]]&]],{n,0,30}]

A358377 Numbers k such that the k-th standard ordered rooted tree is a generalized Bethe tree (counted by A003238).

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 9, 11, 16, 17, 32, 37, 43, 64, 128, 129, 137, 171, 256, 257, 293, 512, 529, 683, 1024, 1025, 2048, 2185, 2341, 2731, 4096, 8192, 10923, 16384, 16913, 18725, 32768, 32769, 32897, 34953, 43691, 65536, 65537, 131072, 131329, 149797, 174763
Offset: 1

Views

Author

Gus Wiseman, Nov 14 2022

Keywords

Comments

We define the n-th standard ordered rooted tree to be obtained by taking the (n-1)-th composition in standard order (graded reverse-lexicographic, A066099) as root and replacing each part with its own standard ordered rooted tree. This ranking is an ordered variation of Matula-Goebel numbers, giving a bijective correspondence between positive integers and unlabeled ordered rooted trees.
A generalized Bethe tree is an unlabeled rooted tree where all branches directly under the same root are equal.

Examples

			The terms together with their corresponding ordered rooted trees begin:
    1: o
    2: (o)
    3: ((o))
    4: (oo)
    5: (((o)))
    8: (ooo)
    9: ((oo))
   11: ((o)(o))
   16: (oooo)
   17: ((((o))))
   32: (ooooo)
   37: (((o))((o)))
   43: ((o)(o)(o))
   64: (oooooo)
  128: (ooooooo)
  129: ((ooo))
  137: ((oo)(oo))
  171: ((o)(o)(o)(o))
		

Crossrefs

These trees are counted by A003238.
The unordered version is A214577, also counted by A003238.
A000081 counts unlabeled rooted trees, ranked by A358378.
A358371 and A358372 count leaves and nodes in standard ordered rooted trees.
A358374 ranks ordered identity trees, counted by A032027.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join @@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    srt[n_]:=If[n==1,{},srt/@stc[n-1]];
    Select[Range[1000],FreeQ[srt[#],[_]?(!SameQ@@#&)]&]

A032027 Number of planted planar trees (n+1 nodes) where any 2 subtrees extending from the same node are different.

Original entry on oeis.org

1, 1, 1, 3, 5, 13, 35, 95, 255, 715, 2081, 6003, 17645, 52127, 155863, 468129, 1415521, 4301055, 13134789, 40275109, 123970669, 382919917, 1186475687, 3686899725, 11487023793, 35876838669, 112304155021, 352276801491
Offset: 1

Views

Author

Keywords

Examples

			From _Gus Wiseman_, Nov 15 2022: (Start)
The a(1) = 1 through a(6) = 13 ordered rooted identity trees (ranked by A358374):
  o  (o)  ((o))  ((o)o)   (((o))o)   (((o)o)o)
                 (o(o))   (((o)o))   ((o(o))o)
                 (((o)))  ((o(o)))   (o((o)o))
                          (o((o)))   (o(o(o)))
                          ((((o))))  ((((o)))o)
                                     ((((o))o))
                                     ((((o)o)))
                                     (((o))(o))
                                     (((o(o))))
                                     ((o)((o)))
                                     ((o((o))))
                                     (o(((o))))
                                     (((((o)))))
(End)
		

Crossrefs

The unordered version is A004111, ranked by A276625.
These trees (ordered rooted identity) are ranked by A358374.

Programs

  • Mathematica
    aot[n_]:=If[n==1,{{}},Join@@Table[Tuples[aot/@c],{c,Join@@Permutations/@IntegerPartitions[n-1]}]];
    Table[Length[Select[aot[n],FreeQ[#,[_]?(!UnsameQ@@#&)]&]],{n,1,10}] (* Gus Wiseman, Nov 15 2022 *)
  • PARI
    AGK(v)={apply(p->subst(serlaplace(y^0*p),y,1), Vec(prod(k=1, #v, (1 + x^k*y + O(x*x^#v))^v[k])-1, -#v))}
    seq(n)={my(v=[1]); for(i=2, n, v=concat([1], AGK(v))); v} \\ Andrew Howroyd, Sep 20 2018

Formula

Shifts left under "AGK" (ordered, elements, unlabeled) transform.

A324524 Numbers where every prime index divides its multiplicity in the prime factorization. Numbers divisible by a power of prime(k)^k for each prime index k.

Original entry on oeis.org

1, 2, 4, 8, 9, 16, 18, 32, 36, 64, 72, 81, 125, 128, 144, 162, 250, 256, 288, 324, 500, 512, 576, 648, 729, 1000, 1024, 1125, 1152, 1296, 1458, 2000, 2048, 2250, 2304, 2401, 2592, 2916, 4000, 4096, 4500, 4608, 4802, 5184, 5832, 6561, 8000, 8192, 9000, 9216
Offset: 1

Views

Author

Gus Wiseman, Mar 07 2019

Keywords

Comments

These are a kind of self-describing numbers (cf. A001462, A304679).
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 prime signature of a number is the multiset of multiplicities (or exponents) in its prime factorization.
Also Heinz numbers of integer partitions in which every part divides its multiplicity (counted by A001156). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
Also products of elements of A062457.

Examples

			The sequence of terms together with their prime indices begins as follows. For example, we have 18: {1,2,2} because 18 = prime(1) * prime(2) * prime(2).
    1: {}
    2: {1}
    4: {1,1}
    8: {1,1,1}
    9: {2,2}
   16: {1,1,1,1}
   18: {1,2,2}
   32: {1,1,1,1,1}
   36: {1,1,2,2}
   64: {1,1,1,1,1,1}
   72: {1,1,1,2,2}
   81: {2,2,2,2}
  125: {3,3,3}
  128: {1,1,1,1,1,1,1}
  144: {1,1,1,1,2,2}
  162: {1,2,2,2,2}
  250: {1,3,3,3}
  256: {1,1,1,1,1,1,1,1}
		

Crossrefs

Range of values of A090884.
Sequences related to self-description: A000002, A001462, A079000, A079254, A276625, A304360.

Programs

  • Maple
    q:= n-> andmap(i-> irem(i[2], numtheory[pi](i[1]))=0, ifactors(n)[2]):
    select(q, [$1..10000])[];  # Alois P. Heinz, Mar 08 2019
  • Mathematica
    Select[Range[1000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>Divisible[k,PrimePi[p]]]&]
    v = Join[{1}, Prime[(r = Range[10])]^r]; n = Length[v]; vmax = 10^4; s = {1}; Do[v1 = v[[k]]; rmax = Floor[Log[v1, vmax]]; s1 = v1^Range[0, rmax]; s2 = Select[Union[Flatten[Outer[Times, s, s1]]], # <= vmax &]; s = Union[s, s2], {k, 2, n}]; Length[s] (* Amiram Eldar, Sep 30 2020 *)

Formula

Closed under multiplication.
Sum_{n>=1} 1/a(n) = Product_{k>=1} 1/(1-prime(k)^(-k)) = 2.26910478689594012492... - Amiram Eldar, Sep 30 2020

A324571 Numbers whose ordered prime signature is equal to the set of distinct prime indices in decreasing order.

Original entry on oeis.org

1, 2, 9, 12, 40, 112, 125, 352, 360, 675, 832, 1008, 2176, 2401, 3168, 3969, 4864, 7488, 11776, 14000, 19584, 29403, 29696, 43776, 44000, 63488, 75600, 104000, 105984, 123201, 151552, 161051, 214375, 237600, 267264, 272000, 335872, 496125, 561600, 571392, 608000
Offset: 1

Views

Author

Gus Wiseman, Mar 08 2019

Keywords

Comments

These are a kind of self-describing numbers (cf. A001462, A304679). The increasing case is A109298.
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 ordered prime signature (A124010) is the sequence of multiplicities (or exponents) in a number's prime factorization, taken in order of the prime base.
Also Heinz numbers of the integer partitions counted by A324572. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
Each finite set of positive integers determines a unique term with those prime indices. For example, corresponding to {1,2,4,5} is 1397088 = prime(1)^5 * prime(2)^4 * prime(4)^2 * prime(5)^1.

Examples

			The sequence of terms together with their prime indices begins as follows. For example, we have 40: {1,1,1,3} because 40 = prime(1) * prime(1) * prime(1) * prime(3).
      1: {}
      2: {1}
      9: {2,2}
     12: {1,1,2}
     40: {1,1,1,3}
    112: {1,1,1,1,4}
    125: {3,3,3}
    352: {1,1,1,1,1,5}
    360: {1,1,1,2,2,3}
    675: {2,2,2,3,3}
    832: {1,1,1,1,1,1,6}
   1008: {1,1,1,1,2,2,4}
   2176: {1,1,1,1,1,1,1,7}
   2401: {4,4,4,4}
   3168: {1,1,1,1,1,2,2,5}
   3969: {2,2,2,2,4,4}
   4864: {1,1,1,1,1,1,1,1,8}
   7488: {1,1,1,1,1,1,2,2,6}
  11776: {1,1,1,1,1,1,1,1,1,9}
  14000: {1,1,1,1,3,3,3,4}
  19584: {1,1,1,1,1,1,1,2,2,7}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[1000],Reverse[PrimePi/@First/@If[#==1,{},FactorInteger[#]]]==Last/@If[#==1,{},FactorInteger[#]]&]
Previous Showing 31-40 of 146 results. Next