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

A282446 Call d a recursive divisor of n iff the p-adic valuation of d is a recursive divisor of the p-adic valuation of n for any prime p dividing d; a(n) gives the number of recursive divisors of n.

Original entry on oeis.org

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

Views

Author

Rémy Sigrist, Feb 15 2017

Keywords

Comments

More informally, the prime tower factorization of a recursive divisor of n can be obtained by removing branches from the prime tower factorization of n (the prime tower factorization of a number is defined in A182318).
A recursive divisor of n is also a divisor of n, hence a(n)<=A000005(n) for any n, with equality iff n is cubefree (i.e. n belongs to A004709).
A recursive divisor of n is also a (1+e)-divisor of n, hence a(n)<=A049599(n) for any n, with equality iff the p-adic valuation of n is cubefree for any prime p dividing n.
This sequence first differs from A049599 at n=256: a(256)=4 whereas A049599(256)=5; note that 256=2^(2^3), and 2^3 is not cubefree.

Examples

			The recursive divisors of 40 are: 1, 2, 5, 8, 10 and 40, hence a(40)=6.
		

Crossrefs

Programs

  • Mathematica
    a[1] = 1; a[n_] := a[n] = Times @@ (1 + a/@ (Last /@ FactorInteger[n])); Array[a, 100] (* Amiram Eldar, Apr 12 2020 *)
  • PARI
    a(n) = my (f=factor(n)); return (prod(i=1, #f~, 1+a(f[i,2])))

Formula

Multiplicative, with a(p^k)=1+a(k) for any prime p and k>0.
a(A014221(n))=n+1 for any n>=0.

A279513 Multiplicative with a(p^k) = p*a(k) for any prime p and k>0.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 6, 6, 10, 11, 12, 13, 14, 15, 8, 17, 12, 19, 20, 21, 22, 23, 18, 10, 26, 9, 28, 29, 30, 31, 10, 33, 34, 35, 24, 37, 38, 39, 30, 41, 42, 43, 44, 30, 46, 47, 24, 14, 20, 51, 52, 53, 18, 55, 42, 57, 58, 59, 60, 61, 62, 42, 12, 65, 66, 67, 68
Offset: 1

Views

Author

Rémy Sigrist, Dec 13 2016

Keywords

Comments

To compute a(n): multiply (with multiplicity) each prime number appearing in the prime tower factorization of n (see A182318 for definition).
a(n) = n if n is squarefree.
a(n) <= A000026(n) for any n>0.
First differs from A000026 at n=256: a(256)=12 and A000026(256)=16.
If n = p_1 * p_2 * ... * p_k with p_1, p_2, ..., p_k primes, then a(p_1 ^ p_2 ^ ... ^ p_k) = n.

Examples

			a(6!) = a(2^(2^2)*3^2*5) = 2*2*2*3*2*5 = 240.
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Section 2.33 Hall-Montgomery Constant, p. 207.

Crossrefs

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=1, 1,
          mul(i[1]*a(i[2]), i=ifactors(n)[2]))
        end:
    seq(a(n), n=1..100);  # Alois P. Heinz, Aug 22 2020
  • Mathematica
    a[n_] := a[n] = If[n==1, 1, Times @@ (#[[1]] a[#[[2]]]& /@ FactorInteger[n] )]; Array[a, 256] (* Jean-François Alcover, Mar 31 2017 *)
  • PARI
    a(n) =  my (f=factor(n)); return (prod(i=1, #f~, f[i,1]*a(f[i,2])))

Formula

Sum_{k=1..n} a(k) ~ (1/2) * c * n^2, where c = Product_{p prime} (1 - 1/p^2 + (p-1)*Sum_{k>=2} a(k)/p^(2*k)) = 0.8351076361... (Gilman and Tschiersch, 1993; Finch, 2003; the constant was calculated by Kevin Ford). - Amiram Eldar, Nov 04 2022

A284456 Numbers such that there is no smaller number with the same factorization shape (see Comments for details).

Original entry on oeis.org

1, 2, 4, 6, 12, 16, 30, 36, 48, 60, 64, 144, 180, 192, 210, 240, 420, 576, 720, 900, 960, 1260, 1296, 1680, 2310, 2880, 3600, 4096, 4620, 5040, 5184, 6300, 6480, 6720, 12288, 13860, 14400, 18480, 20160, 25200, 25920, 30030, 32400, 36864, 44100, 45360, 46656
Offset: 1

Views

Author

Rémy Sigrist, Mar 27 2017

Keywords

Comments

We say that two numbers, say X and Y, have the same factorization shape iff X and Y have the same number of distinct prime factors, say x_1, ..., x_k and y_1, ..., y_k, and there is a permutation f on {1,..,k} such that, for any i between 1 and k, the x_i-adic valuation of X has the same factorization shape as the y_f(i)-adic valuation of Y.
This sequence is a subsequence of A279686 (two numbers with the same prime tower factorization class also have the same factorization shape).
This sequence is a subsequence of the products of primorial numbers (A025487).
This sequence is a supersequence of the primorial numbers (A002110).
The factorization shape of n can be identified with the rooted tree underlying the prime tower factorization of n (see A182318 for the definition of prime tower factorization); for example:
(2) o
| |
12 = 2^2*3 => (2) (3) => o o
\ / \ /
* O
Here are the sets corresponding to some factorization shapes:
- Shape "1": the number 1 (this is the only finite set),
- Shape "2": the prime numbers (A000040),
- Shape "4": the prime powers of prime numbers (A053810),
- Shape "6": the squarefree semiprimes (A006881),
- Shape "16": numbers of the form p^q^r, for p,q,r primes (A217709),
- Shape "30": the sphenic numbers (A007304).
If n belongs to this sequence, then 2^n belongs to this sequence.
If n_1 >= ... >= n_k belong to this sequence, then Product_{i=1..k} prime(i)^n_i belongs to this sequence.
This sequence is not a subsequence of A220219 (48 belongs to this sequence, hence 2^48 belongs to this sequence; but 48+1 is not prime, so 2^48 does not belong to A220219; in fact, a(9)=48 is the first term of the sequence not one less than a prime, and a(681)=2^48 is the first term of this sequence not in A220219).
All terms, except the initial term 1, are even.
If a(n) <= 2^a(m), then the p-adic valuation of a(n) is <= a(m) for any prime p; this property implies that, provided you know the first m terms, you can generate all terms up to 2^a(m) by enumerating the products of primorials <= 2^a(m) with exponents in {a(1), ..., a(m)}; hence, starting with the initial term a(1)=1, after n iterations, you have all terms <= A014221(n).

Crossrefs

A279686 Numbers that are the least integer of a prime tower factorization equivalence class (see Comments for details).

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 18, 30, 36, 40, 48, 60, 64, 72, 81, 90, 108, 144, 162, 180, 192, 200, 210, 225, 240, 256, 280, 320, 324, 360, 405, 420, 432, 450, 500, 512, 540, 576, 600, 630, 648, 720, 768, 810, 900, 960, 1260, 1280, 1296, 1350, 1400, 1536, 1575, 1600
Offset: 1

Views

Author

Rémy Sigrist, Dec 16 2016

Keywords

Comments

The prime tower factorization of a number is defined in A182318.
We say that two numbers, say n and m, belong to the same prime tower factorization equivalence class iff there is a permutation of the prime numbers, say f, such that replacing each prime p by f(p) in the prime tower factorization of n leads to m.
The notion of prime tower factorization equivalence class can be seen as a generalization of the notion of prime signature; thereby, this sequence can be seen as an equivalent of A025487.
This sequence contains all primorial numbers (A002110).
This sequence contains A260548.
This sequence contains the terms > 0 in A014221.
If n appears in the sequence, then 2^n appears in the sequence.
If n appears in the sequence and k>=0, then A002110(k)^n appears in the sequence.
With the exception of term 1, this sequence contains no term from A182318.
Odd numbers appearing in this sequence: 1, 81, 225, 405, 1575, 2025, 2835, 6125, 10125, 11025, 14175, 15625, 16875, 17325, 31185, 33075, 50625, 67375, 70875, 99225, ...
Here are some prime tower factorization equivalence classes:
- Class 1: the number one (the only finite equivalence class),
- Class p: the prime numbers (A000040),
- Class p*q: the squarefree semiprimes (A006881),
- Class p^p: the numbers of the form p^p with p prime (A051674),
- Class p^q: the numbers of the form p^q with p and q distinct primes,
- Class p*q*r: the sphenic numbers (A007304),
- Class p*q*r*s: the products of four distinct primes (A046386),
- Class p*q*r*s*t: the products of five distinct primes (A046387),
- Class p*q*r*s*t*u: the products of six distinct primes (A067885).

Examples

			2 is the least number of the form p with p prime, hence 2 appears in the sequence.
6 is the least number of the form p*q with p and q distinct primes, hence 6 appears in the sequence.
72 is the least number of the form p^q*q^p with p and q distinct primes, hence 72 appears in the sequence.
36000 is the least number of the form p^q*q^r*r^p with p, q and r distinct primes, hence 36000 appears in the sequence.
		

Crossrefs

A284695 Greatest prime number appearing in the prime tower factorization of n (when n > 1); a(1) = 1.

Original entry on oeis.org

1, 2, 3, 2, 5, 3, 7, 3, 3, 5, 11, 3, 13, 7, 5, 2, 17, 3, 19, 5, 7, 11, 23, 3, 5, 13, 3, 7, 29, 5, 31, 5, 11, 17, 7, 3, 37, 19, 13, 5, 41, 7, 43, 11, 5, 23, 47, 3, 7, 5, 17, 13, 53, 3, 11, 7, 19, 29, 59, 5, 61, 31, 7, 3, 13, 11, 67, 17, 23, 7, 71, 3, 73, 37, 5
Offset: 1

Views

Author

Rémy Sigrist, Apr 01 2017

Keywords

Comments

See A182318 for the definition of the prime tower factorization of a number.
a(n) >= A006530(n) for any n>0.

Examples

			8 = 2^3, hence a(8) = max(2, 3) = 3.
		

Crossrefs

Formula

a(n) = A006530(A279513(n)) for any n>0.

A300955 In the prime tower factorization of n, replace 2's with 3's and 3's with 2's.

Original entry on oeis.org

1, 3, 2, 27, 5, 6, 7, 9, 8, 15, 11, 54, 13, 21, 10, 7625597484987, 17, 24, 19, 135, 14, 33, 23, 18, 125, 39, 4, 189, 29, 30, 31, 243, 22, 51, 35, 216, 37, 57, 26, 45, 41, 42, 43, 297, 40, 69, 47, 15251194969974, 343, 375, 34, 351, 53, 12, 55, 63, 38, 87, 59
Offset: 1

Views

Author

Rémy Sigrist, Mar 17 2018

Keywords

Comments

The prime tower factorization of a number is defined in A182318.
This sequence is a self-inverse multiplicative permutation of the natural numbers.
This sequence has infinitely many fixed points (A300957); for any k > 0, at least one of k or 2^k * 3^a(k) is a fixed point.
This sequence is a recursive version of A182318.
This sequence has connections with A300948.

Examples

			a(6) = a(2 * 3) = 3 * 2 = 6.
a(16) = a(2 ^ 2 ^ 2) = 3 ^ 3 ^ 3 = 7625597484987.
		

Crossrefs

Cf. A064614, A182318, A300948, A300957 (fixed points).

Programs

  • Maple
    a:= n-> `if`(n=1, 1, mul(`if`(i[1]=2, 3, `if`(i[1]=3,
                 2, i[1]))^a(i[2]), i=ifactors(n)[2])):
    seq(a(n), n=1..60);  # Alois P. Heinz, Mar 17 2018
  • Mathematica
    a[n_] := If[n == 1, 1, Product[If[i[[1]] == 2, 3, If[i[[1]] == 3,
       2, i[[1]]]]^a[i[[2]]], {i, FactorInteger[n]}]];
    Table[a[n], {n, 1, 60}] (* Jean-François Alcover, Jan 30 2025, after Alois P. Heinz *)
  • PARI
    a(n) = my (f=factor(n)); prod(i=1, #f~, my (p=f[i,1]); if (p==2, 3, p==3, 2, p)^a(f[i,2]))

Formula

Multiplicative with a(p^k) = A064614(p)^a(k).
a(a(n)) = n.

A279510 Multiplicative with a(p(i)^j) = p(i+1)^a(j) for i-th prime p(i) and j>0.

Original entry on oeis.org

1, 3, 5, 27, 7, 15, 11, 243, 125, 21, 13, 135, 17, 33, 35, 7625597484987, 19, 375, 23, 189, 55, 39, 29, 1215, 343, 51, 3125, 297, 31, 105, 37, 2187, 65, 57, 77, 3375, 41, 69, 85, 1701, 43, 165, 47, 351, 875, 87, 53, 38127987424935, 1331, 1029, 95, 459, 59
Offset: 1

Views

Author

Rémy Sigrist, Dec 13 2016

Keywords

Comments

This sequence establishes a bijection between the natural numbers and A182318.
Records: 1, 3, 5, 27, 243, 7625597484987, 38127987424935, 53379182394909, 7450580596923828125, 22351741790771484375, ..., .
Records appear at: 1, 2, 3, 4, 8, 16, 48, 80, 81, 162, 256, 768, 1280, 1792, 2304, 6400, 6561, 13122, 26244, ..., . Robert G. Wilson v, Dec 14 2016

Examples

			a(6) = a(2*3) = 3*5 = 15
a(16) = a(2^(2^2)) = 3^(3^3) = 7625597484987
		

Crossrefs

Programs

  • Mathematica
    a[n_] := Block[{fi = FactorInteger@ n}, Times @@ (NextPrime[#[[1]]]^a[#[[2]]] & /@ fi)]; a[1] = 1; Array[a, 60] (* Robert G. Wilson v, Dec 14 2016 *)
  • PARI
    a(n) =  my (f=factor(n)); return (prod(i=1, #f~, nextprime(1+f[i,1])^a(f[i,2])))

A284694 Least prime number appearing in the prime tower factorization of n (when n > 1); a(1) = 1.

Original entry on oeis.org

1, 2, 3, 2, 5, 2, 7, 2, 2, 2, 11, 2, 13, 2, 3, 2, 17, 2, 19, 2, 3, 2, 23, 2, 2, 2, 3, 2, 29, 2, 31, 2, 3, 2, 5, 2, 37, 2, 3, 2, 41, 2, 43, 2, 2, 2, 47, 2, 2, 2, 3, 2, 53, 2, 5, 2, 3, 2, 59, 2, 61, 2, 2, 2, 5, 2, 67, 2, 3, 2, 71, 2, 73, 2, 2, 2, 7, 2, 79, 2, 2
Offset: 1

Views

Author

Rémy Sigrist, Apr 01 2017

Keywords

Comments

See A182318 for the definition of the prime tower factorization of a number.
a(n) <= A020639(n) for any n>0.

Examples

			9 = 3^2, hence a(9) = min(3, 2) = 2.
		

Crossrefs

Formula

a(n) = A020639(A279513(n)) for any n>0.

A284761 a(n) = gcd(A279513(n), A279513(n+1)).

Original entry on oeis.org

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

Views

Author

Rémy Sigrist, Apr 02 2017

Keywords

Comments

Two consecutive numbers, say n and n+1, cannot share a prime factor (gcd(n, n+1)=1). However, their prime tower factorizations can share some prime numbers; this is the case iff a(n)>1 (see A182318 for the definition of the prime tower factorization of a number).
If p is prime, then a(p-1) = a(p) = 1.
If p is an odd prime, then a(p^2) = 2.
This sequence contains a multiple of p for any prime p:
- let m = A074792(p)^p-1,
- m is a multiple of p, hence p divides A279513(m),
- m+1 = A074792(p)^p, hence p divides A279513(m+1),
- hence p divides gcd(A279513(m), A279513(m+1)) = a(m).
This sequence contains infinitely many distinct values; see A284821 for these distinct values in order of appearance, and A284822 for the corresponding indexes.

Examples

			a(8) = gcd(A279513(8), A279513(9)) = gcd(A279513(2^3), A279513(3^2)) = gcd(2*3, 3*2) = 6.
		

Crossrefs

A287957 Table read by antidiagonals: T(n, k) = greatest common recursive divisor of n and k; n > 0 and k > 0.

Original entry on oeis.org

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

Views

Author

Rémy Sigrist, Jun 03 2017

Keywords

Comments

We use the definition of recursive divisor given in A282446.
More informally, the prime tower factorization of T(n, k) is the intersection of the prime tower factorizations of n and k (the prime tower factorization of a number is defined in A182318).
This sequence has connections with the classical GCD (A003989).
For any i > 0, j > 0 and k > 0:
- T(i, j) = 1 iff gcd(i, j) = 1,
- A007947(T(i, j)) = A007947(gcd(i, j)),
- T(i, j) >= 1,
- T(i, j) <= min(i, j),
- T(i, j) <= gcd(i, j),
- T(i, 1) = 1,
- T(i, i) = i,
- T(i, j) = T(j, i) (the sequence is commutative),
- T(i, T(j, k)) = T(T(i, j), k) (the sequence is associative),
- T(i, i*j) <= i,
- if gcd(i, j) = 1 then T(i*j, k) = T(i, k) * T(j, k) (the sequence is multiplicative),
- T(i, 2*i) = A259445(i).
See also A287958 for the LCM equivalent.

Examples

			Table starts:
n\k|    1   2   3   4   5   6   7   8   9   10
---+-----------------------------------------------
1  |    1   1   1   1   1   1   1   1   1    1  ...
2  |    1   2   1   2   1   2   1   2   1    2  ...
3  |    1   1   3   1   1   3   1   1   3    1  ...
4  |    1   2   1   4   1   2   1   2   1    2  ...
5  |    1   1   1   1   5   1   1   1   1    5  ...
6  |    1   2   3   2   1   6   1   2   3    2  ...
7  |    1   1   1   1   1   1   7   1   1    1  ...
8  |    1   2   1   2   1   2   1   8   1    2  ...
9  |    1   1   3   1   1   3   1   1   9    1  ...
10 |    1   2   1   2   5   2   1   2   1   10  ...
...
T(4, 8) = T(2^2, 2^3) = 2.
		

Crossrefs

Programs

  • PARI
    T(n,k) = my (g=factor(gcd(n,k))); return (prod(i=1, #g~, g[i,1]^T(valuation(n, g[i,1]), valuation(k, g[i,1]))))
Showing 1-10 of 31 results. Next