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 41-50 of 284 results. Next

A316470 Matula-Goebel numbers of unlabeled rooted RPMG-trees, meaning the Matula-Goebel numbers of the branches of any non-leaf node are relatively prime.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 14, 16, 18, 24, 26, 28, 32, 36, 38, 42, 48, 52, 54, 56, 64, 72, 74, 76, 78, 84, 86, 96, 98, 104, 106, 108, 112, 114, 122, 126, 128, 144, 148, 152, 156, 162, 168, 172, 178, 182, 192, 196, 202, 208, 212, 214, 216, 222, 224, 228, 234, 244, 252
Offset: 1

Views

Author

Gus Wiseman, Jul 04 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. A number is in the sequence iff it is 1 or its prime indices are relatively prime and already belong to the sequence.

Examples

			The sequence of all RPMG-trees preceded by their Matula-Goebel numbers begins:
   1: o
   2: (o)
   4: (oo)
   6: (o(o))
   8: (ooo)
  12: (oo(o))
  14: (o(oo))
  16: (oooo)
  18: (o(o)(o))
  24: (ooo(o))
  26: (o(o(o)))
  28: (oo(oo))
  32: (ooooo)
  36: (oo(o)(o))
  38: (o(ooo))
  42: (o(o)(oo))
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[1000],Or[#==1,And[GCD@@primeMS[#]==1,And@@#0/@primeMS[#]]]&]

A324736 Number of subsets of {1...n} containing all prime indices of the elements.

Original entry on oeis.org

1, 2, 3, 4, 7, 9, 15, 22, 43, 79, 127, 175, 343, 511, 851, 1571, 3141, 4397, 8765, 13147, 25243, 46843, 76795, 115171, 230299, 454939, 758203, 1516363, 2916079, 4356079, 8676079, 12132079, 24264157, 45000157, 73800253, 145685053, 291369853, 437054653, 728424421
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.
Also the number of subsets of {1...n} containing no prime indices of the non-elements up to n.

Examples

			The a(0) = 1 through a(6) = 15 subsets:
  {}  {}   {}     {}       {}         {}           {}
      {1}  {1}    {1}      {1}        {1}          {1}
           {1,2}  {1,2}    {1,2}      {1,2}        {1,2}
                  {1,2,3}  {1,4}      {1,4}        {1,4}
                           {1,2,3}    {1,2,3}      {1,2,3}
                           {1,2,4}    {1,2,4}      {1,2,4}
                           {1,2,3,4}  {1,2,3,4}    {1,2,6}
                                      {1,2,3,5}    {1,2,3,4}
                                      {1,2,3,4,5}  {1,2,3,5}
                                                   {1,2,3,6}
                                                   {1,2,4,6}
                                                   {1,2,3,4,5}
                                                   {1,2,3,4,6}
                                                   {1,2,3,5,6}
                                                   {1,2,3,4,5,6}
An example for n = 18 is {1,2,4,7,8,9,12,16,17,18}, whose elements have the following prime indices:
   1: {}
   2: {1}
   4: {1,1}
   7: {4}
   8: {1,1,1}
   9: {2,2}
  12: {1,1,2}
  16: {1,1,1,1}
  17: {7}
  18: {1,2,2}
All of these prime indices {1,2,4,7} belong to the subset, as required.
		

Crossrefs

The strict integer partition version is A324748. The integer partition version is A324753. The Heinz number version is A290822. An infinite version is A324698.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],SubsetQ[#,PrimePi/@First/@Join@@FactorInteger/@DeleteCases[#,1]]&]],{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]));
    ((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 15 2019

Extensions

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

A006508 a(n+1) = a(n)-th composite number, with a(0) = 1.

Original entry on oeis.org

1, 4, 9, 16, 26, 39, 56, 78, 106, 141, 184, 236, 299, 374, 465, 570, 696, 843, 1014, 1212, 1441, 1708, 2014, 2365, 2769, 3226, 3749, 4343, 5016, 5774, 6630, 7596, 8676, 9897, 11259, 12784, 14482, 16383, 18502, 20847, 23458, 26354, 29562, 33112, 37041, 41370
Offset: 0

Views

Author

Keywords

Comments

Generated by a sieve: start with natural numbers, remove those terms which occupy positions which are prime, leaving 1,4,6,8,9,10,12,14,15,16,18,...; remove those terms whose positions are primes plus one; leaving 1,4,9,12,15,16,18,...; remove those whose positions are primes plus two; continue. - Robert G. Wilson v, Jan 06 2008
Number of terms <= 10^k, k=0..: 1, 3, 8, 18, 34, 54, 80, 110, 147, 188, 235, 287, 345, 407, 475, ..., . - Robert G. Wilson v, Jan 06 2008
The first occurrence of a k-almost prime or 0 if not present: 1, 0, 4, 78, 16, 696, 5016, 95920, 46144, 10236034900, 600374208, 1613472, 21565696412400, 0, 24323590656, ..., . - Robert G. Wilson v, Jan 06 2008

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    a006508 n = a006508_list !! n
    a006508_list = iterate a002808 1  -- Reinhard Zumkeller, Oct 24 2011
  • Mathematica
    Composite[n_Integer] := FixedPoint[n + PrimePi@# + 1 &, n + PrimePi@n + 1]; NestList[ Composite@# &, 1, 45] (* Labos Elemer *)
    With[{c=Select[Range[42000],CompositeQ]},NestList[c[[#]]&,1,45]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Mar 03 2018 *)

Extensions

Edited by Robert G. Wilson v, Jan 06 2008

A324704 Lexicographically earliest sequence containing 1 and all numbers > 2 divisible by prime(m) for some m already in the sequence.

Original entry on oeis.org

1, 4, 6, 7, 8, 10, 12, 13, 14, 16, 17, 18, 19, 20, 21, 22, 24, 26, 28, 29, 30, 32, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 46, 48, 49, 50, 51, 52, 53, 54, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 70, 71, 72, 73, 74, 76, 77, 78, 79, 80, 82, 84
Offset: 1

Views

Author

Gus Wiseman, Mar 11 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:
   1: {}
   4: {1,1}
   6: {1,2}
   7: {4}
   8: {1,1,1}
  10: {1,3}
  12: {1,1,2}
  13: {6}
  14: {1,4}
  16: {1,1,1,1}
  17: {7}
  18: {1,2,2}
  19: {8}
  20: {1,1,3}
  21: {2,4}
  22: {1,5}
  24: {1,1,1,2}
  26: {1,6}
  28: {1,1,4}
		

Crossrefs

Programs

A245703 Permutation of natural numbers: a(1) = 1, a(p_n) = A014580(a(n)), a(c_n) = A091242(a(n)), where p_n = n-th prime, c_n = n-th composite number and A014580(n) and A091242(n) are binary codes for n-th irreducible and n-th reducible polynomials over GF(2), respectively.

Original entry on oeis.org

1, 2, 3, 4, 7, 5, 11, 6, 8, 12, 25, 9, 13, 17, 10, 14, 47, 18, 19, 34, 15, 20, 31, 24, 16, 21, 62, 26, 55, 27, 137, 45, 22, 28, 42, 33, 37, 23, 29, 79, 59, 35, 87, 71, 36, 166, 41, 58, 30, 38, 54, 44, 61, 49, 32, 39, 99, 76, 319, 46, 91, 108, 89, 48, 200, 53, 97, 75, 40, 50, 203, 70, 67, 57, 78, 64, 43, 51
Offset: 1

Views

Author

Antti Karttunen, Aug 02 2014

Keywords

Comments

All the permutations A091202, A091204, A106442, A106444, A106446, A235041 share the same property that primes (A000040) are mapped bijectively to the binary representations of irreducible GF(2) polynomials (A014580) but while they determine the mapping of composites (A002808) to the corresponding binary codes of reducible polynomials (A091242) by a simple multiplicative rule, this permutation employs index-recursion also in that case.

Crossrefs

Programs

  • PARI
    allocatemem(123456789);
    a014580 = vector(2^18);
    a091242 = vector(2^22);
    isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
    i=0; j=0; n=2; while((n < 2^22), if(isA014580(n), i++; a014580[i] = n, j++; a091242[j] = n); n++)
    A245703(n) = if(1==n, 1, if(isprime(n), a014580[A245703(primepi(n))], a091242[A245703(n-primepi(n)-1)]));
    for(n=1, 10001, write("b245703.txt", n, " ", A245703(n)));
    
  • Scheme
    ;; With memoization-macro definec.
    (definec (A245703 n) (cond ((= 1 n) n) ((= 1 (A010051 n)) (A014580 (A245703 (A000720 n)))) (else (A091242 (A245703 (A065855 n))))))

Formula

a(1) = 1, a(p_n) = A014580(a(n)) and a(c_n) = A091242(a(n)), where p_n is the n-th prime, A000040(n) and c_n is the n-th composite, A002808(n).
a(1) = 1, after which, if A010051(n) is 1 [i.e. n is prime], then a(n) = A014580(a(A000720(n))), otherwise a(n) = A091242(a(A065855(n))).
As a composition of related permutations:
a(n) = A245702(A135141(n)).
a(n) = A091204(A245821(n)).
Other identities. For all n >= 1, the following holds:
a(A007097(n)) = A091230(n). [Maps iterates of primes to the iterates of A014580. Permutation A091204 has the same property]
A091225(a(n)) = A010051(n). [Maps primes to binary representations of irreducible GF(2) polynomials, A014580, and nonprimes to union of {1} and the binary representations of corresponding reducible polynomials, A091242. The permutations A091202, A091204, A106442, A106444, A106446 and A235041 have the same property.]

A331683 One and all numbers of the form 2^k * prime(j) for k > 0 and j already in the sequence.

Original entry on oeis.org

1, 4, 8, 14, 16, 28, 32, 38, 56, 64, 76, 86, 106, 112, 128, 152, 172, 212, 214, 224, 256, 262, 304, 326, 344, 424, 428, 448, 512, 524, 526, 608, 622, 652, 688, 766, 848, 856, 886, 896, 1024, 1048, 1052, 1154, 1216, 1226, 1244, 1304, 1376, 1438, 1532, 1696
Offset: 1

Views

Author

Gus Wiseman, Jan 30 2020

Keywords

Comments

Also Matula-Goebel numbers of lone-child-avoiding rooted trees at with at most one non-leaf branch under any given vertex. A rooted tree is lone-child-avoiding if there are no unary branchings. The Matula-Goebel number of a rooted tree is the product of primes indexed by the Matula-Goebel numbers of the branches of the root, which gives a bijective correspondence between positive integers and unlabeled rooted trees.
Also Matula-Goebel numbers of lone-child-avoiding locally disjoint semi-identity trees. Locally disjoint means no branch of any vertex overlaps a different (unequal) branch of the same vertex. In a semi-identity tree, all non-leaf branches of any given vertex are distinct.

Examples

			The sequence of all lone-child-avoiding rooted trees with at most one non-leaf branch under any given vertex together with their Matula-Goebel numbers begins:
    1: o
    4: (oo)
    8: (ooo)
   14: (o(oo))
   16: (oooo)
   28: (oo(oo))
   32: (ooooo)
   38: (o(ooo))
   56: (ooo(oo))
   64: (oooooo)
   76: (oo(ooo))
   86: (o(o(oo)))
  106: (o(oooo))
  112: (oooo(oo))
  128: (ooooooo)
  152: (ooo(ooo))
  172: (oo(o(oo)))
  212: (oo(oooo))
  214: (o(oo(oo)))
  224: (ooooo(oo))
		

Crossrefs

These trees counted by number of vertices are A212804.
The semi-lone-child-avoiding version is A331681.
The non-semi-identity version is A331871.
Lone-child-avoiding rooted trees are counted by A001678.
Matula-Goebel numbers of lone-child-avoiding rooted trees are A291636.
Unlabeled semi-identity trees are counted by A306200, with Matula-Goebel numbers A306202.
Locally disjoint rooted trees are counted by A316473.
Matula-Goebel numbers of locally disjoint rooted trees are A316495.
Lone-child-avoiding locally disjoint rooted trees by leaves are A316697.

Programs

  • Maple
    N:= 10^4: # for terms <= N
    S:= {1}:
    with(queue):
    Q:= new(1):
    while not empty(Q) do
      r:= dequeue(Q);
      p:= ithprime(r);
      newS:= {seq(2^i*p,i=1..ilog2(N/p))} minus S;
      S:= S union newS;
      for s in newS do enqueue(Q,s) od:
    od:
    sort(convert(S,list)); # Robert Israel, Feb 05 2020
  • Mathematica
    uryQ[n_]:=n==1||MatchQ[FactorInteger[n],({{2,},{p,1}}/;uryQ[PrimePi[p]])|({{2,k_}}/;k>1)];
    Select[Range[100],uryQ]

Formula

Intersection of A291636, A316495, and A306202.

A358576 Matula-Goebel numbers of rooted trees whose node-height equals their number of internal (non-leaf) nodes.

Original entry on oeis.org

9, 15, 18, 21, 23, 30, 33, 35, 36, 39, 42, 46, 47, 49, 51, 57, 60, 61, 66, 70, 72, 73, 77, 78, 83, 84, 87, 91, 92, 93, 94, 95, 98, 102, 111, 113, 114, 119, 120, 122, 123, 129, 132, 133, 137, 140, 144, 146, 149, 151, 154, 156, 159, 166, 167, 168, 174, 177, 181
Offset: 1

Views

Author

Gus Wiseman, Nov 25 2022

Keywords

Comments

The Matula-Goebel number of a rooted tree is the product of primes indexed by the Matula-Goebel numbers of the branches of its root, which gives a bijective correspondence between positive integers and unlabeled rooted trees.
Node-height is the number of nodes in the longest path from root to leaf.

Examples

			The terms together with their corresponding rooted trees begin:
   9: ((o)(o))
  15: ((o)((o)))
  18: (o(o)(o))
  21: ((o)(oo))
  23: (((o)(o)))
  30: (o(o)((o)))
  33: ((o)(((o))))
  35: (((o))(oo))
  36: (oo(o)(o))
  39: ((o)(o(o)))
  42: (o(o)(oo))
  46: (o((o)(o)))
  47: (((o)((o))))
  49: ((oo)(oo))
  51: ((o)((oo)))
  57: ((o)(ooo))
  60: (oo(o)((o)))
  61: ((o(o)(o)))
		

Crossrefs

The version for edge-height is A209638.
Square trees are A358577, counted by A358589, ordered A358590.
The version for leaves instead of height is A358578, counted by A185650.
These trees are counted by A358587, ordered A358588.
A000081 counts rooted trees, ordered A000108.
A034781 counts rooted trees by nodes and height.
A055277 counts rooted trees by leaves, ordered A001263.

Programs

  • Mathematica
    MGTree[n_]:=If[n==1,{},MGTree/@Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Count[MGTree[#],[_],{0,Infinity}]==Depth[MGTree[#]]-1&]

Formula

A358552(a(n)) = A342507(a(n)).

A324698 Lexicographically earliest sequence containing 2 and all numbers > 1 whose prime indices already belong to the sequence.

Original entry on oeis.org

2, 3, 5, 9, 11, 15, 23, 25, 27, 31, 33, 45, 47, 55, 69, 75, 81, 83, 93, 97, 99, 103, 115, 121, 125, 127, 135, 137, 141, 155, 165, 197, 207, 211, 225, 235, 243, 249, 253, 257, 275, 279, 291, 297, 309, 341, 345, 347, 363, 375, 379, 381, 405, 411, 415, 419, 423
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}
   9: {2,2}
  11: {5}
  15: {2,3}
  23: {9}
  25: {3,3}
  27: {2,2,2}
  31: {11}
  33: {2,5}
  45: {2,2,3}
  47: {15}
  55: {3,5}
  69: {2,9}
  75: {2,3,3}
  81: {2,2,2,2}
  83: {23}
  93: {2,11}
  97: {25}
  99: {2,2,5}
		

Crossrefs

Programs

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

Original entry on oeis.org

1, 2, 3, 5, 8, 13, 19, 30, 54, 96, 156, 248, 440, 688, 1120, 1864, 3664, 5856, 11232, 16896, 31296, 53952, 91008, 137472, 270528, 516720, 863088, 1710816, 3173856, 4836672, 9329472, 14897376, 29788128, 52256448, 88429248, 166037184, 331648704, 497685888, 829449600
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(6) = 19 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}    {2,4}
                              {3,4}    {2,5}
                              {4,5}    {3,4}
                              {2,4,5}  {3,6}
                                       {4,5}
                                       {4,6}
                                       {5,6}
                                       {2,4,5}
                                       {3,4,6}
                                       {4,5,6}
An example for n = 20 is {5,6,7,9,10,12,14,15,16,19,20}, with prime indices:
   5: {3}
   6: {1,2}
   7: {4}
   9: {2,2}
  10: {1,3}
  12: {1,1,2}
  14: {1,4}
  15: {2,3}
  16: {1,1,1,1}
  19: {8}
  20: {1,1,3}
None of these prime indices {1,2,3,4,8} belong to the subset, as required.
		

Crossrefs

The maximal case is A324743. The strict integer partition version is A324751. The integer partition version is A324756. The Heinz number version is A324758. An infinite version is A304360.

Programs

  • Mathematica
    Table[Length[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]));
    ((k,b)->if(k>#p, 1, my(t=self()(k+1,b)); if(!bitand(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

A358580 Difference between the number of leaves and the number of internal (non-leaf) nodes in the rooted tree with Matula-Goebel number n.

Original entry on oeis.org

1, 0, -1, 1, -2, 0, 0, 2, -1, -1, -3, 1, -1, 1, -2, 3, -1, 0, 1, 0, 0, -2, -2, 2, -3, 0, -1, 2, -2, -1, -4, 4, -3, 0, -1, 1, 0, 2, -1, 1, -2, 1, 0, -1, -2, -1, -3, 3, 1, -2, -1, 1, 2, 0, -4, 3, 1, -1, -2, 0, -1, -3, 0, 5, -2, -2, 0, 1, -2, 0, -1, 2, -1, 1, -3
Offset: 1

Views

Author

Gus Wiseman, Nov 25 2022

Keywords

Comments

The Matula-Goebel number of a rooted tree is the product of primes indexed by the Matula-Goebel numbers of the branches of its root, which gives a bijective correspondence between positive integers and unlabeled rooted trees.

Examples

			The Matula-Goebel number of ((ooo(o))) is 89, and it has 4 leaves and 3 internal nodes, so a(89) = 1.
		

Crossrefs

Zeros are A358578, counted by A185650 (ordered A358579).
Positions of positive terms are counted by A358581, negative A358582.
Positions of nonnegative terms are counted by A358583, nonpositive A358584.
A000081 counts rooted trees, ordered A000108.
A034781 counts trees by nodes and height.
A055277 counts trees by nodes and leaves, ordered A001263.

Programs

  • Mathematica
    MGTree[n_]:=If[n==1,{},MGTree/@Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Count[MGTree[n],{},{0,Infinity}]-Count[MGTree[n],[_],{0,Infinity}],{n,100}]

Formula

a(n) = A109129(n) - A342507(n).
Previous Showing 41-50 of 284 results. Next