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

A072411 LCM of exponents in prime factorization of n, a(1) = 1.

Original entry on oeis.org

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

Views

Author

Labos Elemer, Jun 17 2002

Keywords

Comments

The sums of the first 10^k terms, for k = 1, 2, ..., are 14, 168, 1779, 17959, 180665, 1808044, 18084622, 180856637, 1808585068, 18085891506, ... . Apparently, the asymptotic mean of this sequence is limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1.8085... . - Amiram Eldar, Sep 10 2022

Examples

			n = 288 = 2*2*2*2*2*3*3; lcm(5,2) = 10; Product(5,2) = 10, max(5,2) = 5;
n = 180 = 2*2*3*3*5; lcm(2,2,1) = 2; Product(2,2,1) = 4; max(2,2,1) = 2; it deviates both from maximum of exponents (A051903, for the first time at n=72), and product of exponents (A005361, for the first time at n=36).
For n = 36 = 2*2*3*3 = 2^2 * 3^2 we have a(36) = lcm(2,2) = 2.
For n = 72 = 2*2*2*3*3 = 2^3 * 3^2 we have a(72) = lcm(2,3) = 6.
For n = 144 = 2^4 * 3^2 we have a(144) = lcm(2,4) = 4.
For n = 360 = 2^3 * 3^2 * 5^1 we have a(360) = lcm(1,2,3) = 6.

Links

Crossrefs

Cf. A028234, A067029, A072412-A072414, A273058, A284569, A290103.
Similar sequences: A001222 (sum of exponents), A005361 (product), A051903 (maximal exponent), A051904 (minimal exponent), A052409 (gcd of exponents), A267115 (bitwise-and), A267116 (bitwise-or), A268387 (bitwise-xor).
Cf. also A055092, A060131.
Differs from A290107 for the first time at n=144.
After the initial term, differs from A157754 for the first time at n=360.

Programs

  • Mathematica
    Table[LCM @@ Last /@ FactorInteger[n], {n, 2, 100}] (* Ray Chandler, Jan 24 2006 *)
  • PARI
    a(n) = lcm(factor(n)[,2]); \\ Michel Marcus, Mar 25 2017
  • Python
    from sympy import lcm, factorint
def a(n):
    l=[]
    f=factorint(n)
    for i in f: l+=[f[i],]
    return lcm(l)
print([a(n) for n in range(1, 151)]) # Indranil Ghosh, Mar 25 2017

Formula

a(1) = 1; for n > 1, a(n) = lcm(A067029(n), a(A028234(n))). - Antti Karttunen, Aug 09 2016
From Antti Karttunen, Aug 22 2017: (Start)
a(n) = A284569(A156552(n)).
a(n) = A290103(A181819(n)).
a(A289625(n)) = A002322(n).
a(A290095(n)) = A055092(n).
a(A275725(n)) = A060131(n).
a(A260443(n)) = A277326(n).
a(A283477(n)) = A284002(n). (End)

Extensions

a(1) = 1 prepended and the data section filled up to 120 terms by Antti Karttunen, Aug 09 2016

A290107 a(1) = 1; for n > 1, a(n) = product of distinct exponents in the prime factorization of n.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Aug 13 2017

Keywords

Examples

			For n = 36 = 2^2 * 3^2, the only distinct exponent that occurs is 2, thus a(36) = 2.
For n = 144 = 2^4 * 3^2, the distinct exponents are 2 and 4, thus a(144) = 2*4 = 8.
For n = 4500 = 2^2 * 3^2 * 5^3, the distinct exponents are 2 and 3, thus a(4500) = 2*3 = 6.

Links

Crossrefs

Cf. A156061, A181819.
Differs from A005361 for the first time at n=36.
Differs from A072411 for the first time at n=144, and also from A157754 for the second time (after the initial term).

Programs

  • Mathematica
    Table[If[n == 1, 1, Apply[Times, Union[FactorInteger[n][[All, -1]] ]]], {n, 120}] (* Michael De Vlieger, Aug 14 2017 *)
  • PARI
    A290107(n) = factorback(vecsort((factor(n)[, 2]), ,8));
  • Scheme
    (define (A290107 n) (A156061 (A181819 n)))

Formula

a(n) = A156061(A181819(n)).

A158378 a(1) = 0, a(n) = gcd(A051904(n), A051903(n)) for n >= 2.

Original entry on oeis.org

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

Views

Author

Jaroslav Krizek, Mar 17 2009

Keywords

Comments

a(n) for n >= 2 equals GCD of minimum and maximum exponents in the prime factorization of n.
a(n) for n >= 2 it deviates from A052409(n), first different term is a(10800) = a(2^4*3^3*5^2), a(10800) = gcd(2,4) = 2, A052409(10800) = gcd(2,3,4) = 1.

Examples

			For n = 12 = 2^2 * 3^1 we have a(12) = gcd(2,1) = 1.

Links

Crossrefs

Cf. A157754, A051904, A051903, A052409.

Programs

  • Mathematica
    Table[GCD @@ {Min@ #, Max@ #} - Boole[n == 1] &@ FactorInteger[n][[All, -1]], {n, 100}] (* Michael De Vlieger, Jul 12 2017 *)
  • PARI
    A051903(n) = if((1==n),0,vecmax(factor(n)[, 2]));
A051904(n) = if((1==n),0,vecmin(factor(n)[, 2]));
A158378(n) = gcd(A051903(n),A051904(n)); \\ Antti Karttunen, Jul 12 2017
  • PARI
    a(n) = if(n == 1, 0, my(e = factor(n)[,2]); gcd(vecmin(e), vecmax(e))); \\ Amiram Eldar, Sep 11 2024

Formula

For n >= 2 holds: a(n)*A157754(n) = A051904(n)*A051903(n).
a(1) = 0, a(p) = 1, a(pq) = 1, a(pq...z) = 1, a(p^k) = k, for p = primes (A000040), pq = product of two distinct primes (A006881), pq...z = product of k (k > 2) distinct primes p, q, ..., z (A120944), p^k = prime powers (A000961(n) for n > 1), k = natural numbers (A000027).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1. - Amiram Eldar, Sep 11 2024
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.