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 22 results. Next

A196050 Number of edges in the rooted tree with Matula-Goebel number n.

Original entry on oeis.org

0, 1, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 4, 4, 5, 4, 5, 5, 5, 5, 5, 6, 5, 6, 5, 5, 6, 5, 5, 6, 5, 6, 6, 5, 5, 6, 6, 5, 6, 5, 6, 7, 6, 6, 6, 6, 7, 6, 6, 5, 7, 7, 6, 6, 6, 5, 7, 6, 6, 7, 6, 7, 7, 5, 6, 7, 7, 6, 7, 6, 6, 8, 6, 7, 7, 6, 7, 8, 6, 6, 7, 7, 6, 7, 7, 6, 8, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 6, 7, 7, 7, 8, 6, 6, 8, 6, 8
Offset: 1

Views

Author

Emeric Deutsch, Sep 27 2011

Keywords

Comments

The Matula-Goebel number of a rooted tree is defined in the following recursive manner: to the one-vertex tree there corresponds the number 1; to a tree T with root degree 1 there corresponds the t-th prime number, where t is the Matula-Goebel number of the tree obtained from T by deleting the edge emanating from the root; to a tree T with root degree m>=2 there corresponds the product of the Matula-Goebel numbers of the m branches of T.
a(n) is, for n >= 2, the number of prime function prime(.) = A000040(.) operations in the complete reduction of n. See the W. Lang link with a list of the reductions for n = 2..100, where a curly bracket notation {.} is used for prime(.). - Wolfdieter Lang, Apr 03 2018
From Gus Wiseman, Mar 23 2019: (Start)
Every positive integer has a unique factorization (encoded by A324924) into factors q(i) = prime(i)/i, i > 0. For example:
11 = q(1) q(2) q(3) q(5)
50 = q(1)^3 q(2)^2 q(3)^2
360 = q(1)^6 q(2)^3 q(3)
In this factorization, a(n) is the number of factors counted with multiplicity. For example, a(11) = 4, a(50) = 7, a(360) = 10.
(End)
From Antti Karttunen, Oct 23 2023: (Start)
Totally additive with a(prime(n)) = 1 + a(n).
Number of iterations of A366385 (or equally, of A366387) needed to reach 1.
(End)

Examples

			a(7) = 3 because the rooted tree with Matula-Goebel number 7 is the rooted tree Y.
a(2^m) = m because the rooted tree with Matula-Goebel number 2^m is the star tree with m edges.
		

Crossrefs

Programs

  • Haskell
    import Data.List (genericIndex)
    a196050 n = genericIndex a196050_list (n - 1)
    a196050_list = 0 : g 2 where
       g x = y : g (x + 1) where
         y = if t > 0 then a196050 t + 1 else a196050 r + a196050 s
             where t = a049084 x; r = a020639 x; s = x `div` r
    -- Reinhard Zumkeller, Sep 03 2013
    
  • Maple
    with(numtheory): a := proc (n) local r, s: r := proc (n) options operator, arrow: op(1, factorset(n)) end proc: s := proc (n) options operator, arrow: n/r(n) end proc: if n = 1 then 0 elif bigomega(n) = 1 then 1+a(pi(n)) else a(r(n))+a(s(n)) end if end proc: seq(a(n), n = 1 .. 110);
  • Mathematica
    a[1] = 0; a[n_?PrimeQ] := a[n] = 1 + a[PrimePi[n]]; a[n_] := Total[#[[2]] * a[#[[1]] ]& /@ FactorInteger[n]];
    Array[a, 110] (* Jean-François Alcover, Nov 16 2017 *)
    difac[n_]:=If[n==1,{},With[{i=PrimePi[FactorInteger[n][[1,1]]]},Sort[Prepend[difac[n*i/Prime[i]],i]]]];
    Table[Length[difac[n]],{n,100}] (* Gus Wiseman, Mar 23 2019 *)
  • PARI
    a(n) = my(f=factor(n)); [self()(primepi(p))+1 |p<-f[,1]]*f[,2]; \\ Kevin Ryde, May 28 2021
    
  • Python
    from functools import lru_cache
    from sympy import isprime, primepi, factorint
    @lru_cache(maxsize=None)
    def A196050(n):
        if n == 1 : return 0
        if isprime(n): return 1+A196050(primepi(n))
        return sum(e*A196050(p) for p, e in factorint(n).items()) # Chai Wah Wu, Mar 19 2022

Formula

a(1)=0; if n = prime(t) (the t-th prime), then a(n)=1 + a(t); if n = r*s (r,s>=2), then a(n)=a(r)+a(s). The Maple program is based on this recursive formula.
a(n) = A061775(n) - 1.
a(n) = A109129(n) + A366388(n) = A109082(n) + A358729(n). - Antti Karttunen, Oct 23 2023

A109129 Width (i.e., number of non-root vertices having degree 1) of the rooted tree with Matula-Goebel number n.

Original entry on oeis.org

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

Views

Author

Keith Briggs, Aug 17 2005

Keywords

Comments

The Matula-Goebel number of a rooted tree is defined in the following recursive manner: to the one-vertex tree there corresponds the number 1; to a tree T with root degree 1 there corresponds the t-th prime number, where t is the Matula-Goebel number of the tree obtained from T by deleting the edge emanating from the root; to a tree T with root degree m >= 2 there corresponds the product of the Matula-Goebel numbers of the m branches of T.
A non-root vertex having degree 1 is called a leaf.
Every positive integer has a unique factorization (see A324924) into factors q(i) = prime(i)/i for i > 0. The number of ones in this factorization is a(n). For example, 30 = q(1)^3 q(2)^2 q(3), so a(30) = 3. - Gus Wiseman, Mar 23 2019

Examples

			a(7)=2 because the rooted tree with Matula-Goebel number 7 is the rooted tree Y.
a(2^m) = m because the rooted tree with Matula-Goebel number 2^m is a star with m edges.
		

Crossrefs

Programs

  • Haskell
    import Data.List (genericIndex)
    a109129 n = genericIndex a109129_list (n - 1)
    a109129_list = 0 : 1 : g 3 where
       g x = y : g (x + 1) where
         y = if t > 0 then a109129 t else a109129 r + a109129 s
             where t = a049084 x; r = a020639 x; s = x `div` r
    -- Reinhard Zumkeller, Sep 03 2013
    
  • Maple
    with(numtheory): a := proc (n) local r, s: r := proc (n) options operator, arrow: op(1, factorset(n)) end proc: s := proc (n) options operator, arrow: n/r(n) end proc: if n = 1 then 0 elif n = 2 then 1 elif bigomega(n) = 1 then a(pi(n)) else a(r(n))+a(s(n)) end if end proc: seq(a(n), n = 1 .. 110);
  • Mathematica
    Nest[Function[{a, n}, Append[a, If[PrimeQ@ n, a[[PrimePi@ n]], Total@ Map[#2 a[[#1]] & @@ # &, FactorInteger[n]] ]]] @@ {#, Length@ # + 1} &, {0, 1}, 105] (* Michael De Vlieger, Mar 24 2019 *)
  • PARI
    ML(n) = if( n==1, 1, my(f=factor(n)); sum(k=1,matsize(f)[1],ML(primepi(f[k,1]))*f[k,2])) ;
    A109129(n) = if( n==1, 0, ML(n) ); \\ François Marques, Mar 16 2021
    
  • Python
    from functools import lru_cache
    from sympy import primepi, isprime, factorint
    @lru_cache(maxsize=None)
    def A109129(n):
        if n <= 2: return n-1
        if isprime(n): return A109129(primepi(n))
        return sum(e*A109129(p) for p, e in factorint(n).items()) # Chai Wah Wu, Mar 19 2022

Formula

a(1)=0; a(2)=1; if n = p(t) (= the t-th prime) and t >= 2, then a(n) = a(t); if n = rs (r, s >= 2), then a(n) = a(r) + a(s). The Maple program is based on this recursive formula.
The Gutman et al. references contain a different recursive formula.

Extensions

Typo in formula fixed by Reinhard Zumkeller, Sep 03 2013

A324923 Number of distinct factors in the factorization of n into factors q(i) = prime(i)/i, i > 0.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 20 2019

Keywords

Comments

Also the number of distinct proper terminal subtrees of the rooted tree with Matula-Goebel number n. See illustrations in A061773.

Examples

			The factorization 22 = q(1)^2 q(2) q(3) q(5) has four distinct factors, so a(22) = 4.
		

Crossrefs

Programs

  • Mathematica
    difac[n_]:=If[n==1,{},With[{i=PrimePi[FactorInteger[n][[1,1]]]},Sort[Prepend[difac[n*i/Prime[i]],i]]]];
    Table[Length[Union[difac[n]]],{n,100}]
  • PARI
    A006530(n) = if(1==n, n, my(f=factor(n)); f[#f~, 1]);
    A324923(n) = { my(lista = List([]), gpf, i); while(n > 1, gpf=A006530(n); i = primepi(gpf); n /= gpf; n *= i; listput(lista,i)); #Set(lista); }; \\ Antti Karttunen, Oct 23 2023

Formula

a(n) = A317713(n) - 1.
a(n) = A196050(n) - A366386(n). - Antti Karttunen, Oct 23 2023

Extensions

Data section extended up to a(108) by Antti Karttunen, Oct 23 2023

A324924 Irregular triangle read by rows giving the factorization of n into factors q(i) = prime(i)/i, i > 0.

Original entry on oeis.org

1, 1, 2, 1, 1, 1, 2, 3, 1, 1, 2, 1, 1, 4, 1, 1, 1, 1, 1, 2, 2, 1, 1, 2, 3, 1, 2, 3, 5, 1, 1, 1, 2, 1, 1, 2, 6, 1, 1, 1, 4, 1, 1, 2, 2, 3, 1, 1, 1, 1, 1, 1, 4, 7, 1, 1, 1, 2, 2, 1, 1, 1, 8, 1, 1, 1, 2, 3, 1, 1, 1, 2, 4, 1, 1, 2, 3, 5, 1, 1, 2, 2, 9
Offset: 1

Views

Author

Gus Wiseman, Mar 20 2019

Keywords

Comments

Row n is the multiset of Matula-Goebel numbers of all proper terminal subtrees of the rooted tree with Matula-Goebel number n. For example, the rooted tree with Matula-Goebel number 1362 is (o(o)((oo)(oo))), with proper terminal subtrees {o,o,o,o,o,o,(o),(oo),(oo),((oo)(oo))}, which have Matula-Goebel numbers {1,1,1,1,1,1,2,4,4,49}, which is row 1362, as required.

Examples

			Triangle begins:
  {}
  1
  1  2
  1  1
  1  2  3
  1  1  2
  1  1  4
  1  1  1
  1  1  2  2
  1  1  2  3
  1  2  3  5
  1  1  1  2
  1  1  2  6
  1  1  1  4
  1  1  2  2  3
  1  1  1  1
  1  1  4  7
  1  1  1  2  2
  1  1  1  8
  1  1  1  2  3
  1  1  1  2  4
  1  1  2  3  5
  1  1  2  2  9
For example, row 65 is {1,1,1,2,2,3,6} because 65 = q(1)^3 q(2)^2 q(3) q(6).
		

Crossrefs

Programs

  • Mathematica
    difac[n_]:=If[n==1,{},With[{i=PrimePi[FactorInteger[n][[1,1]]]},Sort[Prepend[difac[n*i/Prime[i]],i]]]];
    Table[difac[n],{n,30}]

A324931 Integers in the list of quotients of positive integers by their product of prime indices.

Original entry on oeis.org

1, 2, 4, 3, 8, 6, 16, 12, 7, 5, 32, 9, 24, 14, 10, 64, 18, 48, 28, 20, 128, 36, 19, 13, 21, 15, 96, 27, 56, 40, 256, 72, 38, 26, 11, 42, 30, 192, 54, 112, 17, 80, 512, 144, 76, 52, 22, 84, 60, 384, 49, 23, 35, 53, 108, 37, 224, 25, 57, 39, 34, 160, 63, 1024
Offset: 1

Views

Author

Gus Wiseman, Mar 21 2019

Keywords

Comments

These quotients are given by A324932(n)/A324933(n).
This is a permutation of the positive integers, with inverse A324934.

Examples

			The sequence of quotients n/A003963(n) begins: 1, 2, 3/2, 4, 5/3, 3, 7/4, 8, 9/4, 10/3, 11/5, 6, 13/6, 7/2, 5/2, 16, ...
		

Crossrefs

Programs

  • Mathematica
    Select[Table[n/Times@@Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>PrimePi[p]^k],{n,100}],IntegerQ]

Formula

a(n) = A324850(n)/A003963(A324850(n)).

A324925 Number of integer partitions y of n such that Product_{i in y} prime(i)/i is an integer.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 5, 5, 5, 8, 9, 11, 17, 19, 21, 28, 32, 40, 51, 57, 67, 83, 96, 118, 142, 160, 189, 224, 260, 307, 363, 412, 479, 561, 649, 749, 874, 997, 1141, 1321, 1518, 1734, 1994, 2274, 2582, 2960, 3374, 3837, 4370, 4950, 5604, 6371, 7208, 8157, 9231, 10392
Offset: 0

Views

Author

Gus Wiseman, Mar 20 2019

Keywords

Comments

The Heinz numbers of these integer partitions are given by A324850.

Examples

			The a(1) = 1 through a(8) = 5 integer partitions:
  (1)  (11)  (21)   (211)   (2111)   (321)     (3211)     (32111)
             (111)  (1111)  (11111)  (411)     (4111)     (41111)
                                     (2211)    (22111)    (221111)
                                     (21111)   (211111)   (2111111)
                                     (111111)  (1111111)  (11111111)
For example, (3,2,1,1) is such a partition because (2/1) * (2/1) * (3/2) * (5/3) = 10 is an integer.
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],IntegerQ[Product[Prime[i]/i,{i,#}]]&]],{n,0,30}]

A324935 Matula-Goebel numbers of rooted trees whose non-leaf terminal subtrees are all different.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 14, 16, 17, 19, 20, 21, 22, 24, 26, 28, 29, 31, 32, 34, 35, 37, 38, 40, 41, 42, 43, 44, 48, 51, 52, 53, 56, 57, 58, 59, 62, 64, 67, 68, 70, 71, 73, 74, 76, 77, 79, 80, 82, 84, 85, 86, 88, 89, 91, 95, 96, 101, 102, 104
Offset: 1

Views

Author

Gus Wiseman, Mar 21 2019

Keywords

Comments

Every positive integer has a unique factorization into factors q(i) = prime(i)/i, i > 0. This sequence consists of all numbers where this factorization has all distinct factors, except possibly for any multiplicity of q(1). For example, 22 = q(1)^2 q(2) q(3) q(5) is in the sequence, while 50 = q(1)^3 q(2)^2 q(3)^2 is not.
The enumeration of these trees by number of vertices is A324936.

Examples

			The sequence of 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))
  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))))
  31: (((((o)))))
		

Crossrefs

Programs

  • Mathematica
    difac[n_]:=If[n==1,{},With[{i=PrimePi[FactorInteger[n][[1,1]]]},Sort[Prepend[difac[n*i/Prime[i]],i]]]];
    Select[Range[100],UnsameQ@@DeleteCases[difac[#],1]&]

A324934 Inverse permutation to A324931.

Original entry on oeis.org

1, 2, 4, 3, 10, 6, 9, 5, 12, 15, 35, 8, 24, 14, 26, 7, 41, 17, 23, 20, 25, 47, 52, 13, 58, 34, 28, 19, 79, 37, 184, 11, 87, 61, 53, 22, 56, 33, 60, 30, 145, 36, 92, 70, 65, 75, 164, 18, 51, 82, 98, 46, 54, 39, 178, 29, 59, 106, 293, 49, 122, 245, 63, 16, 125
Offset: 1

Views

Author

Gus Wiseman, Mar 21 2019

Keywords

Crossrefs

A324936 Number of unlabeled rooted trees with n vertices whose non-leaf terminal subtrees are all different.

Original entry on oeis.org

1, 1, 2, 4, 8, 17, 37, 83, 189, 436, 1014, 2373, 5578, 13156, 31104, 73665, 174665, 414427, 983606, 2334488
Offset: 1

Views

Author

Gus Wiseman, Mar 21 2019

Keywords

Comments

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

Examples

			The a(1) = 1 through a(6) = 17 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)))))
		

Crossrefs

Programs

  • Mathematica
    durt[n_]:=Join@@Table[Select[Union[Sort/@Tuples[durt/@ptn]],UnsameQ@@Cases[#,{},{0,Infinity}]&],{ptn,IntegerPartitions[n-1]}];
    Table[Length[durt[n]],{n,10}]

A325614 Unsorted q-signature of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 12 2019

Keywords

Comments

Every positive integer has a unique q-factorization (encoded by A324924) into factors q(i) = prime(i)/i, i > 0. For example:
11 = q(1) q(2) q(3) q(5)
50 = q(1)^3 q(2)^2 q(3)^2
360 = q(1)^6 q(2)^3 q(3)
Row n lists the nonzero multiplicities in the q-factorization of n, in order of q-index. For example, row 11 is (1,1,1,1) and row 360 is (6,3,1).

Examples

			Triangle begins:
  {}
  1
  1 1
  2
  1 1 1
  2 1
  2 1
  3
  2 2
  2 1 1
  1 1 1 1
  3 1
  2 1 1
  3 1
  2 2 1
  4
  2 1 1
  3 2
  3 1
  3 1 1
		

Crossrefs

Row lengths are A324923.
Row sums are A196050.
Row-maxima are A109129.

Programs

  • Mathematica
    difac[n_]:=If[n==1,{},With[{i=PrimePi[FactorInteger[n][[1,1]]]},Sort[Prepend[difac[n*i/Prime[i]],i]]]];
    Table[Length/@Split[difac[n]],{n,30}]
Showing 1-10 of 22 results. Next